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NOMENCLATURE  OF  IMPORTANT  PARAMETERS 


A  Cross-sectional  area 

b  Width 

CT  Cryogenic  temperature 

D  Tank  diameter;  flexural  rigidity  of  a  plate 

DLF  Dynamic  load  factor 

dof  Degree  of  freedom 

E  Elastic  modulus  of  a  material 

El  Flexural  rigidity  of  a  plate-stiffener  combination 
F  Dynamic  force 

f  Frequency 

G  Shear  modulus  of  a  material 

g  Acceleration  of  gravity 

CA  Shear  rigidity  of  a  material 

H  Tank  height 

h  liquid  filling  in  a  tank;  plate  thickness 

I  Second  moment  of  area  of  stiffener  and  associated  plate;  beam 

section  moment  of  inertia 

K  Axial  spring  constant  of  a  stiffener 

Kp  Nondimensional  force  coefficient 

Kp  Dynamic  load  factor 

Kjn  Dynamic  mass  factor 

Kp  Nondimensional  dynamic  pressure  coefficient 

K'p  Nondimensional  dynamic  pressure  coefficient  corrected  for  oscilla¬ 

tion  amplitude  effects 

Kp  Nondimensional  time  coefficient 
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%  Tank  length;  plywood  specimen  length 

LNG  Liquified  Natural  Gas  (Liquid  Methane) 

LN2  Liquid  nitrogen 

LPG  Liquified  Petroleum  Gas 

m  Mass 

N  Load  applied  over  a  surface 

P  Dynamic  pressure  created  by  liquid  sloshing 

R  Resultant  dynamic  force  created  by  liquid  sloshing  in  a  spherical 

tank 

RF  Reduction  factor 

RT  Room  temperature 

S  Sectional  area  of  stiffener  and  associated  plate 

T  Time ;  period 

Tr  Resonant  period 

tr  Rise  time 

U  Internal  strain  energy 

V  Spherical  tank  volume;  shear  reaction 

W  Static  weight  of  liquid  in  a  tank;  external  work 

w  Deflection  of  a  beam  under  loading 

x  Translational  (horizontal  or  surge)  excitation  amplitude  of  a 

prismatic  or  rectangular  ship  tank 

Y  Lateral  dynamic  force  created  by  liquid  sloshing  in  a  spherical  tank 

y  Translational  (vertical  or  heave)  excitation  amplitude  of  a  prismatic 

or  rectangular  ship  tank 

Z  Vertical  dynamic  force  created  by  liquid  sloshing  in  a  spherical  tank 

?.  Distance  from  a  ship  tank  bottom 

a  Angle  (relative  to  vertical)  of  the  resultant  dynamic  force  (R) 

created  by  liquid  sloshing  in  a  ship  tank 
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NOMENCLATURE  (Cpnt,) 


Y  Specific  weight  of  a  liquid 

<J>  Pitch  or  roll  angle  of  a  ship  tank 

p  Mass  density 

r)  Horizontal  excitation  amplitude  of  a  spherical  ship  tank 
i|)  See  Equation  V-9  ,  Pge.  131, 

a  Stress 

v  Poisson’s  ratio 

a)  Circular  frequency 


I.  INTRODUCTION 


This  report  presents  all  results  of  Project  SR-1251,  "Evaluation  of 
Liquid  Dynamic  Loads  in  Slack  LNG  Cargo  Tanks."  The  study  was  conducted 
under  the  direction  of  the  Ship  Research  Committee  of  the  National  Academy 
of  Sciences. 

The  potential  for  significant  loads  resulting  from  sloshing  liquids 
in  slack  cargo  tanks  has  been  realized  by  the  marine  industry  for  years. 

In  the  past  decade,  with  the  advent  of  the  supertanker  and  large  LNG  ships, 
sloshing  loads  have  been  of  even  greater  concern.  In  the  case  of  LNG 
ships,  certain  operational  constraints  call  for  the  transport  of  liquid 
in  slack  tanks  and,  in  addition,  the  absence  of  tank  internals  results  in 
no  damping  of  the  liquid  motions.  Also,  as  ship  tanks  have  grown  in  size, 
the  probability  of  resonant  sloshing  has  increased  since  resonant  sloshing 
periods  and  ship  motions  more  closely  match.  In  view  of  this,  a  signifi¬ 
cant  activity  has  been  undertaken  by  various  agencies  throughout  the  world 
to  establish  sloshing  loads  in  LNG  cargo  tanks.  Experimental  programs  have 
been  conducted  using  instrumented  scale  model  tanks.  These  test  programs 
have  covered  many  different  ship  tank  geometries,  excitation  amplitudes 
and  frequencies,  and  liquid  fill  depths.  Experiments  were  necessary  be¬ 
cause  large  amplitude  sloshing  is  not  amenable  to  theoretical  analysis. 

To  develop  usable  liquid  dynamic  load  criteria  for  sloshing  requires  the  com¬ 
bining  of  current  experimental  information  with  new  experiments  and  analysis. 

Four  project  tasks  were  undertaken  to  develop  dynamic  load  criteria 
for  slack  tanks.  In  Task  1,  all  test  reports  on  LNG  slosh  loads  were  as¬ 
sembled,  and  the  data  reviewed  and  presented  on  a  uniform  basis.  The  data 
were  summarized  with  regard  to  source,  test  conditions,  and  types  of  mea¬ 
surements  recorded.  To  present  the  data  on  a  uniform  basis,  nondimens ional 
coefficients  were  utilized  where  data  from  the  various  reports  were  pre¬ 
sented  in  the  same  format.  Concurrent  with  this  effort,  the  structural 
details  of  typical  LNG  tank  designs  were  established,  and  methods  and  back¬ 
ground  information  necessary  to  develop  LNG  slosh  load  design  methodology 
were  identified. 

In  Task  2,  scale  model  tanks  were  utilized  to  conduct  additional 
sloshing  experiments  whereby  important  loads  and  test  conditions  not  cur¬ 
rently  covered  in  the  literature  were  investigated.  Scale  model  LNG  ship 
tanks  were  utilized  whereby  dynamic  sloshing  pressures  were  measured  at 
resonant  sloshing  conditions.  Both  single  and  combined  degree  of  freedom 
excitation  experiments  were  conducted  to  evaluate  the  ability  of  single 
degree  of  freedom  tests  to  produce  maximum  impact  loads.  Significant  em¬ 
phasis  was  placed  on  establishing  pressure-time  histories  for  resonant 
sloshing  impact  loads  for  use  in  dynamic  tank  wall  response  analysis. 
Experimental  studies  were  also  conducted  to  reproduce  full-scale  impulsive 
impact  pressures  on  representative  segments  of  membrane  tank  structures 
and  to  measure  the  structural  response. 

In  Task  3,  analytical  efforts  were  undertaken  to  develop  methods  for 
predicting  LNG  tank  wall  and  support  structural  response  to  typical  slosh 


loads.  Analytical  methods  were  used  to  examine  the  stresses  and  deforma¬ 
tions  in  typical  LNG  tank  structures  and  their  supports  when  subjected  to 
dynamic  sloshing  loads.  Design  guidelines  and  design  methods  were  formu¬ 
lated  utilizing  the  analytical  studies  conducted  in  this  task;  the  experi¬ 
mental  studies  from  Task  2;  and  information  on  past  research  and  analyses, 
collected  and  analyzed  in  Task  1. 

Finally,  in  Task  4,  the  composite  of  information  generated  in  the 
project  was  utilized  to  generate  simplified  design  procedures  to  account 
for  dynamic  loads  in  slack  cargo  tanks. 


II.  BACKGROUND 


II. 1  History  of  Slosh  Problem 

Space  Related  Activities 

The  basic  problem  of  determining  the  dynamic  loads  which  result  from 
the  motions  ("sloshing")  of  liquids  in  partially  filled  moving  containers 
was  studied  extensively  in  the  1950's  and  60' s  in  connection  with  the  de¬ 
velopment  of  large  rocket  vehicles  for  the  space  program.  Liquid  sloshing 
in  spherical  and  cylindrical  containers  has  been  studied  for  space  appli¬ 
cations,  both  analytically  and  experimentally. (1)  The  nature  of  slosh 
loading  in  these  types  of  tanks  and  its  prediction  are  probably  better 
understood  than  for  prismatic  tanks,  but  analytical  techniques  for  predict¬ 
ing  large  amplitude  sloshing  are  still  not  fully  developed,  and  such  loads 
are  extremely  important  in  designing  the  support  structure  and  internal 
components  of  ship  tanks.  In  addition,  much  of  the  sloshing  technology 
developed  for  space  applications  is  not  applicable  because  emphasis  was 
placed  on  frequencies  and  total  forces  as  they  related  to  control  system 
requirements,  and,  therefore,  the  effects  of  local  peak  impact  pressure  on 
structural  requirements  were  not  studied  to  any  extent.  Further,  the  ex¬ 
citation  amplitudes  considered  in  space  applications  are  too  small  for 
ship  motion  simulation. 

LNG  Ship  Related  Slosh  Problems 

Current  activity  in  the  design  of  super  tankers  as  well  as  ships  for 
liquified  gas  transport  has  resulted  in  renewed  consideration  of  the  influ¬ 
ence  of  the  contained  liquids  on  cargo  tank  design, (2-75 )  especially  since 
the  probability  of  exciting  a  resonant  slosh  mode  is  increased  in  the  lar¬ 
ger  tanks.  In  many  cases,  the  transport  of  liquid  cargos  in  partially 
filled  tanks  is  prohibited.  However,  several  factors  make  partial  filling 
either  unavoidable  or  highly  attractive.  For  example,  in  the  case  of  liq¬ 
uid  natural  gas  (LNG)  ships,  partially  filled  conditions  are  needed  because 
(1)  chilled-down  liquid  is  required  to  maintain  cold  tanks  on  return  trips, 
(?.)  higher  specific  gravity  liquids  than  LNG  are  transported  in  tanks  de¬ 
signed  for  LNG,  (3)  partial  unloading  is  desirable  when  multi-port  stops 
are  made,  and  (4)  loading  or  unloading  at  sea  creates  significant  time 
periods  at  undesirable  fill  depths.  For  all  liquid  cargo  ships,  partial 
filling  in  ballast  tanks  and  fuel  tanks  occurs,  and  conditions  (3)  and  (4) 
above  for  LNG  ships  also  apply.  Therefore,  the  designer  of  all  types  of 
liquid  carriers  must  be  aware  of  the  consequences  of  liquid  sloshing  and 
be  able  to  predict  the  resulting  loads. 

As  the  ~u'p  designs  have  evolved  (Figure  II-l) ,  several 

imports-  cypes  of  unique  ship  loads  have  been  considered  by  the  designers. 
The  slosh-generated  loads  are  one  of  these  and  have  a  considerable  influ¬ 
ence  on  the  tank  and  support  structure  design.  Several  factors  make  slosh 
loads  more  important  with  regard  to  LNG  ship  design.  A  tank  failure  in  an 
LNG  ship  merits  special  consideration  because  of  (1)  the  risk  of  brittle 
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fracture  of  the  primary  structure  (low-temperature  shock),  (2)  the  expensive 
repair  cost  of  the  complicated  tank  designs  (3)  the  high  out-of-service 
costs,  and  (4)  the  potential  for  large  volume  vapor  release.  Also,  the  com¬ 
plexity  of  the  tank  design  in  LNG  carriers  is  such  that  at  least  some  LNG 
tanks  are  more  susceptible  to  damage  from  slosh  loading  than  tanks  for  trans¬ 
porting  oil  or  other  petroleum  products. 

LNG  Ship  Tank  Designs 

There  are  presently  over  tee.  designs  of  LNG  tanks  that  are  either 
currently  in  use  or  under  major  consideration.  These  tanks  generally  fall 
into  two  categories:  namely,  freestanding  (independent)  and  non-freestand¬ 
ing  (membrane)  tanks  (Figure  II— 2) .  Integral  tanks  used  for  LPG  transport 
are  not  acceptable  for  LNG  since  their  use  is  restricted  to  temperatures 
greater  than  -10°C.  The  freestanding  or  independent  tank  is  usually  of  the 
spherical  or  cylindrical  design,  and  because  of  its  geometry  it  is  amenable 
to  stress  analysis  and  other  conventional  analytical  techniques.  Because 
the  stresses  can  be  calculated,  a  secondary  barrier  system  is  not  required 
as  is  the  case  of  the  non-freestanding  tanks.  An  exception  is  the  pris¬ 
matic  freestanding  tank,  which  does  require  a  secondary  barrier.  Free¬ 
standing  tanks  are  also  easier  to  fabricate,  and  the  insulation  is  easier 
to  install  than  on  other  systems.  One  drawback  to  the  freestanding  design 
is  the  disadvantage  of  requiring  a  larger  ship  per  given  cargo  volume. 

Since  freestanding  tank  walls  can  be  designed  to  withstand  large  impact 
pressures,  the  primary  problem  associated  with  LNG  sloshing  in  freestand¬ 
ing  tanks  results  from  the  slosh  loads  on  internal  components  and  on  the 
tank  support  structure. 

The  second  general  tank  type,  the  non-freestanding  or  membrane  tank, 
is  essentially  built  into  the  ship’s  hold,  making  use  of  the  ship’s  struc¬ 
ture  for  support.  The  membrane  tanks  use  a  thin  internal  layer  of  metal 
to  act  as  a  liquid  barrier  and  are  directly  supported  by  insulation  mater¬ 
ial.  The  insulation  is  applied  directly  to  the  hull  with  no  access  space, 
which  makes  this  type  of  tank  difficult  to  repair  after  material  fracture 
or  other  damage.  Because  of  the  complex  structure,  membrane  tanks  are  not 
amenable  to  analysis.  In  addition,  because  of  this  and  the  thinness  of  the 
membrane,  a  complete  secondary  liquid  barrier  is  required.  The  primary 
problem  associated  with  sloshing  in  membrane  tanks  is  the  potential  damage 
to  the  tank  walls  from  impulsive  slosh  pressures.  Severe  impulsive  slosh 
loads  in  the  membrane  tank  can  occur  at  small  fill  depths  as  a  result  of 
large-amplitude  traveling  wave  impact,  which  is  not  amenable  to  analysis. 

Also,  severe  slosh  loads  can  occur  on  or  near  the  tank  top  as  a  consequence 

of  standing  slosh  waves  in  partially  filled  tanks.  Since  this  type  of  tank 

cannot  be  analyzed  to  determine  its  failure  strength,  special  load  tests 

must  be  performed  on  representative  segments  of  '  le  structure  to  determine 

its  load  bearing  strength.  An  estimate  of  the  "equivalent"  static  slosh 

loads  that  occur  in  these  types  of  tanks  is  then  utilized  to  determine  if 

the  structure  has  the  required  strength.  i 

i 
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FIGURE  I 1-2.  EXAMPLE  LNG  TANK  DESIGNS 


Recorded  LNG  Tank  Damage  from  Sloshing 


As  of  1979,  over  80  years  of  operating  experience  have  been  gained 
with  numerous  ships  of  various  tank  designs.  During  these  years,  several 
structural  problems  have  been  recorded  which  have  resulted  from  slosh  loads 
on  LNG  cargo  tanks. (6)  Slosh  related  loads  causing  tank  damage  have  oc¬ 
curred  on  two  ships  with  membrane  tanks,  the  "Polar  AJ „ska"  and  the  "Arctic 
Tokyo/'  On  the  "Polar  Alaska,"  supports  of  the  electric  cables  supplying 
the  cargo  pumps  were  broken  by  liquid  sloshing  loads.  This  occurred  when 
the  tank  was  approximately  15-20%  full.  The  broken  cable  supports  resulted 
in  damage  to  the  bottom  of  the  membrane  tank.  On  the  "Arctic  Tokyo,"  a 
leak  in  the  number  1  tank  was  caused  from  liquid  sloshing  when  the  tank 
was  about  20%  full.  Inspection  revealed  that  the  leak  was  located,  along 
with  four  deformed  points  in  the  membranes,  in  the  aft  corners  of  the 
transverse  and  longitudinal  bulkheads  at  about  the  liquid  surface  level. 
Subsequent  model  tests(15,43)  performed  on  scale  models  of  the  damaged 
"Arctic  Tokyo"  tank  revealed  that  a  15-30%  fill  depth  with  respect  to  tank 
length  resulted  in  appreciable  impact  loads  from  sloshing.  The  model  tests 
were  not  successful,  however,  in  establishing  peak  impact  pressures(7)  that 
could  have  caused  the  damage.  As  a  result,  additional  work(10)  was  under¬ 
taken  to  investigate  more  thoroughly  all  aspects  of  modeling  LNG  sloshing 
and  to  provide  a  greater  understanding  of  the  slosh  generated  loads  and 
their  implications  to  tank  and  ship  design.  As  a  result  of  these  studies, 
operations  with  partially  filled  tanks  other  than  nearly  fi.ll  or  nearly 
empty  have  been  prohibited. 

In  the  spring  of  1978,  the  first  of  the  125,000  membrane  tank  LNG 
ships  was  put  into  service.  On  one  of  the  early  cargo-laden  voyages,  the 
ship  experienced  heavy  seas  and  the  crew  heard  loud  sloshing  impact  noises 
in  the  cargo  tanks.  After  cargo  discharge  and  subsequent  inspection  of  the 
tanks,  damage  to  the  tank  structure  was  noted.  Subsequent  studies(74 , 75) 
concluded  that  the  damage  was  slosh  induced,  even  though  the  tanks  were 
approximately  95%  full  during  the  voyage.  Thus,  the  designer  must  also  be 
concerned  with  slosh  loads  at  near  full  conditions. 

II. 2  Nature  of  Liquid  Sloshing 

General  Conditions 


In  general,  sloshing  is  affected  by  liquid  fill  depth,  tank  geometry, 
and  tank  motion  (amplitude  and  frequency).  The  liquid  motion  inside  a  tank 
has  an  infinite  number  of  natural  periods,  but  it  is  the  lowest  mode  that 
is  most  likely  to  be  excited  by  the  motions  of  a  ship.  Most  studies  have, 
therefore,  concentrated  on  investigating  forced  harmonic  oscillations  in  the 
vicinity  of  the  lowest  natural  period,  which  is  defined  as  that  predicted 
by  linear  theory.  Nonlinear  effects  result  in  the  frequency  of  maximum  re¬ 
sponse  being  slightly  different  from  the  natural  frequency  and  dependent  on 
amplitude.  The  most  significant  type  of  ship  tank  slosh  loads  occur  with 
large  excitation  amplitudes  where  nonlinear  effects  are  present. 

The  sloshing  phenomena  in  cargo  tanks  that  are  basically  rectangular 
in  shape  can  usually  be  described  by  considering  only  two-dimensional  fluid 
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flow.  Sloshing  in  spherical  or  cylindrical  tanks,  however,  must  usually 
consider  three-dimensional  flow  effects. 

Two-Dimensional  Flow 


Tanks  with  two-dimensional  flow  are  divided  into  two  classes:  low 
and  high  liquid  fill  depths.  The  low  fill  depth  case  is  represented  by 
h/£  <  0.2,  where  h  is  the  still  liquid  depth  and  £  is  the  tank  length  in 
the  direction  of  motion.  The  low  fill  depth  case  is  characterized  by  the 
formation  of  hydraulic  jumps  and  traveling  waves  for  excitation  periods 
around  resonance.  At  higher  fill  depths,  large  standing  waves  are  usually 
formed  in  the  resonant  frequency  range.  When  hydraulic  jumps  or  traveling 
waves  are  present,  extremely  high  impact  pressures  can  occur  on  the  tank 
walls.  Figure  II-3a  shows  typical  pressure  traces  recorded  under  this 
sloshing  condition.  Impact  pressures  typical  of  those  shown  in  Figure 
II-3a  can  also  occur  on  the  tank  top  when  tanks  are  filled  to  the  higher 
fill  depths.  The  pressure  pulses  are  similar  to  those  experienced  in  ship 
slamming,  and  the  pressure  variation  is  neither  harmonic  nor  periodic  since 
the  magnitude  and  duration  of  the  pressure  peaks  vary  from  cycle  to  cycle 
even  though  the  tank  is  experiencing  a  harmonic  oscillation.  Figure  II-3b 
shows  typical  pressure  traces  that  result  when  small  amplitude  sloshing  is 
occurring  away  from  resonance  at  any  fill  depth. 

Three-Dimensional  Flow 

Three-dimensional  flow  occurs  in  spherical  tanks,  usually  in  the 
form  of  a  swirl  mode.(l)  Similar  three-dimensional  effects  can  be  present 
in  cylindrical  or  rectangular  tanks  under  certain  excitation  conditions. 

The  prediction  of  sloshing  forces  in  the  neighborhood  of  resonance  with 
swirling  is  extremely  difficult,  and  experimental  data  obtained  with  scale 
model  tanks  are  usually  needed  to  establish  pressures  and  forces  with  this 
type  of  sloshing. 

Design  Implications 

The  design  of  a  liquid  cargo  tank  to  withstand  the  dynamic  slosh- 
induced  loads  requires  that  the  designer  be  able  to  predict  the  resonant 
slosh  periods  at  different  fill  depths  for  the  required  tank  geometry. 

These  periods  can  then  be  compared  with  the  expected  ship  periods  to  deter¬ 
mine  the  probability  of  resonant  sloshing.  An  estimate  of  the  maximum  dy¬ 
namic  loads  to  be  expected  is  then  made  to  determine  a  proper  design.  Most 
theoretical  analyses  are  not  able  to  predict  slosh  pressures  and  forces  in 
the  neighborhood  of  resonance,  especially  for  the  large  amplitude  excita¬ 
tions  typical  of  a  ship  cargo  tank.  However,  several  theories  are  avail¬ 
able  to  predict  tank  loadings  at  off-resonant,  low  arapltiuae  sloshing  con¬ 
ditions.  Depending  on  the  likelihood  of  resonance  and  the  expected  exci¬ 
tation  amplitudes,  the  designer  can  either  use  theory  or  turn  to  experimen¬ 
tal  model  data  to  provide  the  required  design  information. 

Available  linear  and  nonlinear  theories  are  discussed  in  detail  in 
References  2  and  6.  A  review  of  the  theoretical  efforts  reveals  that  most 
studies  are  limited  to  linear  sloshing,  which  is  valid  only  for  small 
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(a)  Impact  pressure  traces  with  large 
amplitude  resonant  sloshing 


(b)  Pressure  trace  for  non-resonant  sloshing 

FIGURE  I 1-3.  TYPICAL  PRESSURE  WAVEFORMS  ON  TANK  WALLS  WITH  SLOSHING 
LIQUIDS 
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amplitudes  and  frequencies  and  predicts  an  infinite  response  at  resonance. 
Some  nonlinear  theories  are  available  for  specific  tank  geometries,  and 
these  theories  allow  a  prediction  of  slosh-induced  dynamic  forces  and 
pressures  on  tank  structures  at  resonance.  However,  the  nonlinear  theories 
are  also  limited  to  small  excitation  amplitudes  and  cannot  be  used  for  a 
general  tank  shape  or  when  certain  real  sloshing  effects  are  present  such 
as  liquid  impacting  on  the  tank  top.  Therefore,  emphasis  has  been  placed 
on  utilizing  model  test  data  for  predicting  full  scale  tank  loads  for  de¬ 
sign  purposes. 


II. 3  Previous  Studies 


A  significant  number  of  scale  model  studies  have  been  conducted  to 
investigate  sloshing  in  LNG  cargo  tanks.  These  efforts  have  been  under¬ 
taken  primarily  by  three  worldwide  laboratories: 


o  Southwest  Research  Institute  (SwRI) 

o  Det  norske  Veritas  (DnV) 

o  Bureau  Veritas  (BV) 


These  studies  are  reviewed  in  detail  in  Task  1.  Nearly  all  model  tests  to 
date  have  considered  the  six  degrees  of  ship  motion  individually  and  inves¬ 
tigated  sloshing  by  varying  amplitude  and  frequency  harmonically,  usually 
in  heave,  surge,  pitch,  or  roll.  Also,  water  has  been  used  almost  exclu¬ 
sively  as  the  model  liquid.  In  most  studies,  the  scaling  of  impact  load 
data  to  full  scale  for  use  in  tank  design  has  considered  only  Froude  scal¬ 
ing  and  thus  eliminated  any  possible  effects  of  fluid  properties  such  as 
viscosity,  compressibility,  or  vapor  pressure  (cavitation).  Under  these 
assumptions,  preseures  scale  by 


(iff),  •  (sO.  <IW> 

where  the  subscripts  m  and  p  are  for  the  model  and  prototype,  respectively. 
The  periods  between  prototype  and  model  are  given  by 

(T/g7T)  =  (T/i7I)m  (H-2) 

In  scaling  pressure  data,  a  pressure  coefficient  is  defined  as 


K 


P 


P 

Pg*<J> 


(II— 3) 


where  is  the  pitch,  roll,  or  yaw  angle.  For  translation,  is  usually 
replaced  by  the  translational  amplitude,  x.  (See  Figures  IV-1  and  IV-4  for 
definitions  of  <p  and  x,  respectively.) 
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Scaling  Effects 


The  scaling  criteria  that  should  be  used  in  predicting  full-scale 
slosh  loads  from  model  data  are  discussed  in  References  2,  10,  44,  and  74. 

Most  model  studies  have  utilized  Froude  scaling  to  predict  full-scale  loads, 
and  no  allowance  for  fluid  effects  was  considered.  Depending  on  the  cargo 
to  be  carried,  some  of  these  fluid  properties  would  appear  important.  For 
example  LNG  is  transported  at  a  tank  pressure  slightly  above  its  vapor 
pressure,  and,  therefore,  cavitation  and  thermodynamic  (vapor  condensation) 
effects  could  be  important.  Also,  LNG  has  an  extremely  low  viscosity  com¬ 
pared  to  water,  and  therefore  model  tests  using  water  could  produce  noncon¬ 
servative  predictions  of  full-scale  loads  if  the  model  tests  were  over¬ 
damped.  Also,  compressibility  of  the  impacting  liquid/vapor  could  be  im¬ 
portant  in  scaling  the  slosh  loads. 

SwRI  and  DnV  experimental  programs  have  been  conducted  to  determine 
the  effects  of  these  fluid  properties  on  scaling  sloshing  loads.  The  test 
results  indicate  that  fluid  properties  will  have  a  minor  effect  on  scaling 
impact  pressures  when  large  amplitude  sloshing,  typical  of  a  ship  cargo 
tank,  is  present.  To  determine  the  validity  of  Froude  scaling  for  large 
amplitude  sloshing,  full-scale  pressure  measurements  were  recorded  in  a 
partially  water-filled  0B0  tank  under  rolling  motion  (Reference  2).  Subse¬ 
quent  model  tests  in  l/30th  scale  were  conducted  with  the  full-scale  roll 
motions  reproduced  on  the  model .  Model  pressures  converted  to  full-scale 
using  Froude  scaling  showed  excellent  agreement  for  both  the  magnitudes  and 
distributions  of  pressures.  Since  water  was  used  in  both  model  and  full-scale 
tests,  the  effects  of  liquid  viscosity  and  vapor  condensation  were  not  in¬ 
cluded.  However,  an  evaluation  of  these  effects  in  References  10,  44,  and  74 
indicates  they  are  of  small  importance  to  large  amplitude  slosh  scaling. 

As  a  result,  Froude  scaling  is  appropriate,  and  Equations  II-2  and  II-3  are 
used  for  scaling  periods  and  loads,  respectively. 


'  III .  TASK  1  -  DATA  REVIEW  AND  EVALUATION 


The  objective  of  Task  1  of  the  project  was  to  review  and  present 
currently  available  slosh  test  data  on  a  uniform  basis,  to  identify  data 
required  to  develop  design  methods,  and  to  outline  required  experimental 
and  analytical  studies.  The  majority  of  the  research  work  in  LNG  sloshing 
has  been  done  by  SwRI,  Det  norske  Veritas,  and  Bureau  Veritas.  Numerous 
model  tests  have  been  conducted  by  these  groups  and  others  to  study  slosh¬ 
ing  loads  on  LNG  ship  tanks  and  the  effects  of  these  loads  on  the  tanks  and 
on  the  ship  structures.  A  principal  part  of  Task  1  has  been  to  compile  all 
of  this  information  and  present  it  on  a  uniform  basis. 

The  initial  step  in  Task  1  was  to  conduct  a  thorough  literature 
search  to  find  all  information  that  is  presently  available  on  LNG  sloshing 
in  ship  tanks  (including  pressures,  forces,  and  tank  response)  and  to  ob¬ 
tain  information  on  tank  structural  details.  This  search  was  broken  down 
into  three  segments:  (1)  a  manual  search  of  appropriate  journals  and  peri¬ 
odicals,  (2)  a  computer  search  of  the  pertinent  data  bases,  and  (3)  written 
inquiries.  Written  inquiries  were  sent  to 

General  Dynamics /Quincy  Shipbuilding  Division 

Avondale  Shipyards,  Inc. 

ABS/Research  and  Development  Division 

seeking  principally  tank  structural  details  and  the  results  of  tank  analy¬ 
ses  for  sloshing  loads.  In  search  of  sloshing  data  which  might  not  yet 
appear  in  the  open  literature,  inquiries  were  sent  to 

Mitsui  Shipbuilding  and  Engineering  Co.,  Ltd. 

Bureau  Veritas 

Det  norske  Veritas 

In  addition  to  these  written  inquiries,  personnel  at  Newport  News  Shipbuild¬ 
ing,  El  Paso  Marine,  and  Kaverner-Moss  (U.S.  office)  were  contacted  for  in¬ 
formation. 

The  manual  literature  search  consisted  of  a  survey  of  all  journals 
and  periodicals  that  might  contain  information  on  LNG  sloshing  in  ship 
tanks  or  on  the  structural  analysis  of  ship  tanks  for  slosh-induced  loads. 
The  following  list  consists  of  all  the  sources  that  were  reviewed  in  the 
search.  The  majority  of  the  information  was  listed  in  the  Marine  Research 
Information  Service. 

Sources  Reviewed  by  the  Manual  Literature  Search 

1.  Marine  Research  Information  Service 

2.  Marine  Technology 

3.  Shipping  World  and  Shipbuilder 

4.  Royal  Institute  of  Naval  Architecture 

5.  The  Society  of  Naval  Architects  and  Marine  Engineers 
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6.  Tanker  and  Bulk  Carrier 

7.  European  Shipbuilding 

8.  Northeast  Coast  Institute  of  Engineers  and  Shipbuilders 

9.  Ship  Structure  Committee 

10.  International  Ship  Structures  Committee 

11.  United  States  Coast  Guard  Report 

12.  International  Shipbuilding  Progress 

13.  Norwegian  Maritime  Research 

14.  Marine  Engineering  Log 

15.  Naval  Engineers  Journal 

16.  The  Naval  Architect 


The  computer  literature  search  included  searching  several  data  bases 
to  locate  articles  that  contained  "key  words"  that  were  related  to  LNG 
sloshing  in  ship  tanks  and  structural  analysis  of  ship  tanks.  Most  of  the 
articles  that  were  located  by  the  computer  search  were  also  located  by  the 
manual  search.  This  fact  gives  confidence  that  the  search  thoroughly  in¬ 
vestigated  the  literature  that  is  currently  available  on  this  topic.  The 
following  list  consists  of  all  data  bases  that  were  searched  by  computer. 

Sources  Reviewed  by  the  Computer  Literature  Search 

1.  Computerized  Engineering  Index 

2.  Mechanical  Engineering  Information  Service 

3.  National  Technical  Information  Service 

4 .  Oceanic  Abstracts 

5 .  Energy  Line 

A  reference  list  of  all  the  literature  uncovered  during  this  Task  1 
activity  is  included  in  Section  VIII. 


III.l  Scale  Model  Sloshing 


The  reports  identified  by  the  literature  search  that  contain  infor¬ 
mation  on  model  tests  of  LNG  sloshing  in  ship  tanks  are  summarized  in  Table 
III-l.  The  table  includes  information  on  sources  ->f  the  reports,  tank  ge¬ 
ometries,  test  conditions,  formats  of  presented  results,  and  general  com¬ 
ments  and  observations.  The  data  contained  in  these  reports  were  revic  .ed 
and  analyzed.  An  attempt  was  made  to  reduce  all  pertinent  information  from 
the  various  reports  to  a  common  form  for  presentation  in  this  report. 


A  wide  range  of  test  conditions  is  covered  by  the  experimental  stud¬ 
ies  performed  to  date.  Various  parameters  such  as  impact  pressures,  forces, 
and  moments  have  been  measured  during  these  studies.  A  number  of  tank  ge¬ 
ometries,  tank  motions,  and  test  liquids  have  also  been  investigated. 


Due  to  the  complexity  of  the  liquid  sloshing  phenomena  and  differ¬ 
ences  in  the  methods  of  data  acquisition  by  the  various  investigators,  the 
data  from  the  model  studies  contain  a  significant  amount  of  scatter.  For 
instance,  most  of  the  early  model  experiments  recorded  data  for  relatively 
few  sloshing  cycles.  Later,  it  was  determined  (Reference  44)  that  the 
small  sample  sizes  of  the  early  experiments  (N  <  200)  were  not  sufficient 
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TABLE  HI-1.  SUMMARY  OF  MODEL  TANK  LNG  SLOSHING  EXPERIMENTS 

(Concl’d) 


for  accurately  predicting  worst  case  sloshing  pressures.  Therefore,  the 
more  recent  studies  (containing  200  or  more  cycles  of  sloshing  data)  pt?- 
dict  worst  case  sloshing  loads  that  are  considerably  higher  than  earlier 
estimates. 

Because  of  the  inappropriateness  of  some  of  the  model  test  results, 
only  a  select  number  of  references  were  chosen  for  the  data  presented  in 
this  report.  The  predicted  worst,  case  sloshing  loads  from  the  selected 
studies  provide  a  more  realistic  definition  of  loads  than  the  early  model 
test  data.  The  majority  of  the  summarized  data  have  been  extracted  from 
References  9,  10,  44,  64,  and  74. 

Nondimensional  Coefficients 


Of  principal  importance  in  many  tank  analyses  is  the  slosh-induced 
dynamic  impact  pressure  acting  on  the  tank  walls.  Figure  III-l  shows  a 
typical  pressure  trace  and  distribution  of  pressures  taken  from  laboratory 
test  data.  There  are  several  formats  that  can  be  used  to  present  pressure 
information.  First,  there  is  the  average  pressure,  PavG*  This  value  is 
defined  in  Figure  III-l  as  the  average  pressure  for  a  given  number  of 
sloshing  cycles.  Next,  there  is  the  maximum  (worst  case)  pressure,  Pj.jaX* 
which  is  the  highest  pressure  occurring  during  any  one  of  a  given  number 
of  sloshing  cycles.  Finally,  there  is  the  "q%  exceedance  level"  pressure, 

Pq.  This  is  Che  pressure  level  which  is  exceeded  by  q%  of  the  measured 
pressures  in  a  given  sample  size.  For  example,  the  median  pressure  would 
be  the  50%  exceedance  level.  Much  of  the  data  presented  by  Det  norske 
Veritas  is  presented  in  this  form. 

For  simplicity  and  uniformity,  the  pressure  results  presented  in  this 
report  are  in  nondimensional  form.  The  nondimensional  pressure  is  denoted 
as  Kp.  The  pressure  value  used  in  Kp  is  either  the  average  pressure  (AVG), 
maximum  pressure  (MAX),  or  the  "q%  exceedance  level"  pressure  (q)  as  indicated. 
The  pressure  has  been  nondimensional ized  using  Froude  and  Euler  scaling  as 
discussed  in  References  44  and  74.  The  resulting  nondimensional  pressure 
is  defined  by  the  following  relation: 


K 


P 


P 

Pg£<}> 


or 


P 

Pg£(f> 


(III-l) 


vihere  P 
P 
g 
l 

A 

X 


dynamic  wall  impact  pressure 
liquid  density 
acceleration  of  gravity 

tank  length  in  the  direction  of  the  tank  motion 
angular  excitation  amplitude  (in  radians)* 
translational  excitation  amplitude  (in  units  of  length) 

^Single  amplitude  value 


Another  important  parameter  to  the  design  analysis  of  a  tan'.,  is  the 
inertial  force  created  by  the  sloshing  of  liquid  in  the  tank.  A  nondimen¬ 
sional  force,  Kp,  similar  to  the  nondimensional  pressure  has  been  developed 
for  the  presentation  of  experimental  force  data.  Measured  peak  sloshing 
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ill 
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forces  do  not  exhibit  the  cyclic  variations  typical  of  localized  impact 
pressures.  As  a  result,  the  average  peak  force  (F)  measured  over  a  given 
number  of  sloshing  cycles  is  representative  of  the  peak  force  occurring  on 
any  given  cycle.  For  a  prismatic  or  rectangular  tank  geometry,  this  non- 
dimensional  force  is  defined  as 


K 


F 


P8thb*  °r  pgthbCf) 


(III— 2) 


where  F  =  dynamic  force  (only  inertial  forces  caused  by  the  liquid 
are  included) 

h  =  liquid  filling  height 

b  =  tank  width  transverse  to  the  direction  of  the  tank  motion 

For  a  spherical  tank,  dynamic  impact  pressures  are  not  of  critical 
importance  to  tank  design. (24, 44)  However,  the  tank  forces  created  by  the 
inertia  of  the  sloshing  liquid  are  significant  to  tank  design.  The  nondi- 
mensional  force  coefficient  used  for  spherical  tanks  is  defined  as: 


K 


F 


F 

PgV 


(III-3) 


where  F  =>  dynamic  force  (see  Figure  III -9  for  the  definitions  of  the 
various  forces) 

V  =  tank  volume 

From  Froude  scaling,  the  nondimensional  time  is  defined  as: 


Ky  -  T  Sg/I 


(III-4) 


where  T  =  time. 

In  addition  to  resonant  sloshing  periods,  slosh-induced  dynamic 
pressure  rise  times  are  normalized  by  Equation  (III-4) , 

Table  III-2  provides  equations  for  calculating  the  natural  period, 
Tr,  of  liquid  in  rectangular,  spherical,  and  cylindrical  tanks  as  a  func¬ 
tion  of  fill  depth.  These  equations  are  derived  in  Reference  1. 

Statistical  Analysis  of  Slosh-Induced  Pressures 

Because  of  the  random  nature  of  resonant  sloshing  in  an  LNG  tank,  a 
sufficient  sample  size  must  be  determined  that  will  give  a  representative 
distribution  of  peak  sloshing  pressures.  References  44  and  64  have  both 
addressed  this  problem.  These  two  references  used  somewhat  different  ap¬ 
proaches  to  obtain  similar  results. 


TABLE  II 1-2.  EQUATIONS  FOR  CALCULATING  LIQUID  NATURAL  PERIODS 
IN  RECTANGULAR,  SPHERICAL,  AND  CYLINDRICAL  TANKS 


Rectangular: 


V?  tanh  (f  ) 

Spherical  Tank: 

tr  -  ~^(t)  -  (f  )2  £or  (°- 

with  the  values  of  C^  shown  on  the  figure  below. 


h  =  liquid  filling  height 
R  =  tank  radius 
D  =  tank  diameter 

Cylindrical  Tank:  (Vertical  orientation) 


2tt 


~yf 7.664  tanh  7.664  (^j 
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Reference  44  recorded  1000  cycles  of  resonant  sloshing  pressures. 
Then,  the  largest  number  of  consecutive  peak  pressures  less  than  85%  of  the 
maximum  (worst  case)  pressure  was  determined.  This  number  varied  from  140 
to  160  cycles.  Based  on  this  information,  a  sample  size  of  200  cycles  was 
considered  representative.  Histograms  of  the  200-cycle  and  1000-cycle  test 
data  are  pictured  in  Figure  III-2.  The  distributions  are  quite  similar  in 
shape,  indicating  good  repeatability  between  a  200-cycle  sample  and  a  1000- 
cycle  sample. 

In  Reference  64,  various  sample  sizes  up  to  400  cycles  were  evalu¬ 
ated.  It  was  determined  that  200  cycles  were  required  to  give  good  repeat¬ 
ability  of  the  average  pressure  for  a  given  sample  size.  The  results  of 
References  44  and  64  show  that:  (1)  200  cycles  are  adequate  to  define  the 
average  sloshing  pressure  (Pavg)  anc*  (2)  that  the  worst  case  pressure 
(PjjAx)  will  be  at  least  85%  (or  greater)  of  Pj^ax  f°r  1000  cycles  of  slosh¬ 
ing. 


Reference  44  developed  a  theoretical  pressure  distribution  for  pre¬ 
dicting  extreme  pressure  values.  This  distribution  was  based  on  a  three- 
parameter  Weibull  curve  fit.  The  resulting  equation  is: 


where  Q(P)  =  the  probability  that  the  response  variable  will  be 
greater  than  P 
P  =  dynamic  pressure 

A0,A,M  =  experimentally  determined  parameters  of  the  distribution 

Utilizing  this  distribution.  Reference  44  presents  a  procedure  for  estimat¬ 
ing  long-term  distributions  of  slosh-induced  dynamic  impact  pressures  oc¬ 
curring  during  the  life  of  a  ship.  This  procedure  is  based  on  the  statisti¬ 
cal  analyses  of  the  slosh  pressures  and  the  tank  excitation  motions. 


III. 1.1  Prismatic  Tanks 

Model  scale  data  for  slosh- induced  dynamic  wall  pressures  vary  con¬ 
siderably  with  test  conditions.  While  the  general  nature  of  resonant 
sloshing  in  prismatic  tanks  of  various  geometries  is  similar,  the  impact 
pressures  for  a  given  tank  can  vary  significantly  with  wall  location,  ex¬ 
citation  motion  and  amplitude,  liquid  density,  and  fill  level.  Since  it 
is  impossible  to  measure  all  possible  combinations  of  these  parameters,  the 
composite  of  the  sloshing  pressure  data  from  all  sources  is  presented  so 
that  the  range  of  load  coefficients  can  be  established  for  design  purposes. 
Figure  III-3  shows  the  highest  (for  any  measurement  location)  KpAVG  vs  fill 

level  (h/£)  for  a  range  of  excitation  motions  and  amplitudes,  tank  geome¬ 
tries  (all  are  prismatic  or  rectangular),  pressure  transducer  locations, 
and  test  liquids.  All  points  on  the  graph  are  for  resonant  sloshing  con¬ 
ditions.  The  ranges  between  1  and  25.  No  distinct  relation  between 

Kp^G  and  the  liquid  filling  level  is  noted. 
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»  nondiraensional  pressure  coefficient  (averaged  over  nany 
«VG  cycles) 

h  »  liquid  filling  height 

i  »  tank  length 

NOTE: 


Pressures  are  highest  measured  values  for  any  test  conditions 
regardless  of  tank  motion,  tank  geometry,  transducer  location, 
and  test  liquid. 


FIGURE  1 1 1 -3 .  HIGHEST  AVERAGE  NOND I MENS I ONAL  PRESSURE  VS 
NONDIMENS I ONAL  TANK  FILLING  LEVEL  FOR  ALL  MODEL  TESTS 
RUN  TO  DATE 
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Ranges  of  values  compiled  from  the  model- scale  data  are  plot¬ 

ted  in  Figure  III-4.  These  values  are  for  nondimensional  excitation  ampli¬ 
tudes  from  0.05  to  0.30.  It  should  be  noted  that  the  data  points  on  this 
graph  are  worst  case  values  for  at  least  200  cycles  of  resonant  sloshing. 

As  the  filling  level  increases,  the  location  of  the  highest  pressure  moves 
from  the  vertical  walls  to  the  tank  top.  Otherwise,  no  clearly  defined 
trend  of  Kp^^  with  h /l  is  noted.  The  location  of  highest  pressures  on  the 

vertical  walls  is  usually  near  the  static  liquid  filling  level,  while  on 
the  tank  top  the  highest  pressure  occurs  near  the  corners. 

Except  for  one  data  point,  all  Kp^^  values  on  Figure  III-4  range 

between  10  and  125.  The  Kp  value  of  325  is  considered  unusual  since  the 
other  data  points  represent  a  wide  range  of  sloshing  conditions  covering 
several  sources. 

The  form  of  nondimensional  pressure  coefficient  equation,  Equation 
III-l,  reveals  that  the  pressure  P  is  assumed  to  be  a  linear  function  of 
excitation  amplitude,  <J>  (i.e.,  P  =  Kppg£<{>).  In  actuality,  pressures  are 
only  a  linear  function  of  <}>  for  amplitudes  less  than  -0.1.  Figure  III-5 
from  Reference  44,  exhibits  pressure  vs  amplitude  data  showing  this  trend. 

It  is  evident  from  Figure  III-5  that  as  <f>  approaches  0.3,  the  slopes  of  the 
Kpio<f>  vs  curves  approach  horizontal.  This  means  the  influence  of  the  ex¬ 
citation  amplitude  on  the  impact  pressure  diminishes  as  the  amplitude  in¬ 
creases.  The  relation  between  Kp4>  and  <j>  exhibited  by  Figure  III-5  is  typi¬ 
cal  for  all  areas  of  the  tank  that  experience  slosh- induced  impact  pres¬ 
sures.  Thus,  if  pressure  coefficients  obtained  at  experimental  amplitudes 
less  than  0.1  are  used  to  predict  pressures  at  larger  amplitudes,  the  re¬ 
sults  will  be  conservative. 

Figure  III-6  presents  sloshing  force  coefficient  data  summarized 
from  References  9,  10,  and  44  for  a  range  of  amplitudes.  Figure  III-7  showt. 
the  same  data  for  a  given  amplitude  (±8°).  The  data  indicate  that  the  force 
coefficients  decrease  as  the  filling  level  increases.  However,  the  actual 
force  will  not  necessarily  decrease  with  h/2,  as  a  result  of  the  definition 
of  KF  (i.e.,  F  =  Kppg£2  (h/£.)4>) . 

Table  III-2  presents  an  equation  that  predicts  the  resonant  period 
of  the  liquid  in  a  rectangular  (or  prismatic)  tank.  Nondimens ionalizing  the 
equation  from  Table  III-2,  it  becomes: 

T-\/f~)  =  “V  4tt  ftanh  (^l  (HI-6) 

*  /RESONANCE  L  '  ' J 

A  comparison  (Figure  III-8)  of  the  resonant  periods  calculated  from  this 
equation  to  measured  periods  from  the  model  tests  for  a  variety  of  pris¬ 
matic  tanks  indicates  that  the  theoretical  values  correlate  well  with  the 
experimental  results. 
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FIGURE  I I 1-4.  MAXIMUM  PRESSURE  FIGURE  I I 1-5.  PRESSURE  COEFFICIENTS  VS  PITCH 

COEFFICIENT  VS  TANK  FILLING  LEVEL  AMPLITUDE  (FROM  REFERENCE  44) 


1-7.  AVERAGE  FORCE 
FILLING  LEVEL  (Pit 
=  +0.14  Rad.  or  ±8° 


The  information  presented  in  this  section  covers  the  important  con¬ 
siderations  in  the  design  of  a  rectangular  or  prismatic  tank.  Topics  such 
as  scaling  criteria  for  model  tests,  random  excitation  motions,  statistical 
analysis  of  cargo  sloshing,  and  several  other  problems  related  to  cargo 
sloshing  in  LNG  tanks  have  also  been  addressed  by  various  investigators. 
Their  works  are  listed  in  this  report  in  the  Reference  List  of  currently 
available  literature. 

Ill .1.2  Spherical  Tanks 

As  previously  stated,  forces  and  not  impact  pressures  are  of  primary 
concern  in  the  design  of  a  spherical  LNG  ship  tank.  References  24  and  44 
present  a  thorough  set  of  data  on  sloshing  forces  acting  on  a  spherical 
tank.  The  data  presented  in  this  section  are  excerpted  from  these  reports. 

The  pertinent  forces  acting  on  a  spherical  tank  are  defined  in  Fig¬ 
ure  III-9,  Figures  III-10  through  III— 17  show  typical  force  magnitudes  for 
various  tank  filling  levels  and  excitation  amplitudes.  All  data  compiled 
were  obtained  using  a  horizontal  tank  excitation  motion.  Reaction  forces 
produced  by  pitch  and  roll  motions  about  the  tank  centerline  were  not  re¬ 
ported.  The  authors  stated  that  the  reaction  forces  produced  by  pitch  and 
roll  are  small  relative  to  those  produced  by  horizontal  motions. 

The  following  general  trends  were  discovered  during  the  analysis  of 
the  test  data.  The  maximum  lateral  force  usually  occurs  under  the  same 
conditions  as  the  maximum  resultant  force.  The  direction  angle,  a,  of  the 
resultant  force  is  greatest  for  h/D  less  than  0.2.  The  angle  decreases  as 
h/D  increases  from  0.4  to  0.6  and  remains  about  constant  for  h/D  values 
above  0.6. 

The  lateral  force  attains  its  maximum  value  when  h/D  is  greater  than 
0.8  for  large  amplitude  excitation  motion  (large  amplitude  being  d/D  greater 
than  0.05).  The  maximum  lateral  force  is  about  equivalent  to  the  inertial 
force  of  a  full  sphere.  For  small  amplitude  excitations  (r|/D  values  near 
0.01),  the  maximum  lateral  force  occurs  when  h/D  is  about  0.5. 

The  maximum  resultant  force  is  nearly  equal  to  the  weight  of  liquid 
contained  in  a  full  sphere.  This  maximum  occurs  during  the  same  conditions 
that  create  the  maximum  lateral  force.  These  maximum  resulting  and  lateral 
forces  occur  at  relatively  high  excitation  frequencies. 

The  maximum  dynamic  vertical  force  occurs  when  the  filling  level  of 
the  tank  is  about  0.4.  However,  the  maximum  static  plus  dynamic  vertical 
force  is  created  when  the  tank  filling  level  is  greater  than  0.8.  The  maxi¬ 
mum  static  plus  dynamic  force  is  about  equal  in  magnitude  to  the  maximum 
resultant  force. 

The  resonant  slosh  frequencies  for  the  liquid  in  a  spherical  tank 
vary  substantially  depending  on  the  filling  level.  High  filling  levels 
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FIGURE  I I 1-9 .  SLOSH  FORCES  ON  SPHERE  FIGURE  III-10.  RESULTANT  FORCE  ON  SPHERE  VS 

FILL  DEPTH  (Reference  44) 


RE  III— 15 .  VERTICAL  STATIC  +  DYNAMIC 
IRCE  ON  SPHERE  VS  FILL  DEPTH 
eference  44) 


have  resonant  frequencies  high  enough  that  the  ship  motion  energy  in  this 
frequency  range  is  so  lew  that  resonant  sloshing  is  unlikely  to  occur.  For 
lower  filling  levels,  the  liquid  response  often  corresponds  with  the  ship's 
motion,  causing  resonant  sloshing  to  occur  in  the  tank. 

The  resonant  sloshing  periods  predicted  by  the  equation  for  spheri¬ 
cal  tanks  from  Table  III-2  are  shown  in  Figures  III-ll,  III-13,  and  III-16. 

111. 2  Full-' Scale  Sloshing  Data 

Full-scale  sloshing  data  are  very  scarce.  The  LNG  carrier,  the  Ben 
Franklin .  has  been  instrumented  for  recording  slosh-related  information^ 
but  the  ship  has  yet  to  go  into  service.  Two  El  Paso  Marine  Company  car¬ 
riers,  the  Sonatrach  and  the  Arzew.  have  also  been  instrumented,  but  no 
data  are  available  at  this  time.  Full-scale  data  on  the  sloshing  of  fuel 
oil  in  a  tank  with  no  internal  structures  are  presented  in  Reference  42; 
however,  this  information  is  only  of  very  limited  value  to  the  problem  of 
LNG  sloshing  in  ship  tanks.  Full-scale  impact  pressures  measured  with  wa¬ 
ter  sloshing  in  an  0B0  tank  are  compared  to  model  scale  pressures  in  Refer¬ 
ence  2.  The  model  data  were  obtained  in  a  geometrically  similar  1/30  scale 
tank  with  the  recorded  full-scale  roll  motions  reproduced  on  the  model. 
Pressures  in  model  scale  were  converted  to  full  scale  using  Froude  scaling. 
A  comparison  of  predicted  to  actual  pressures  showed  excellent  agreement 
for  both  the  magnitudes  and  distributions  of  pressures. 

111. 3  Review  of  Tank  Structural  Detail 

A  general  description  of  LNG  tanks  was  given  in  Section  II.  The 
importance  of  sloshing-induced  forces  in  the  design  of  the  tank  depends 
largely  upon  the  tank  type.  For  this  discussion  of  tank  structural  de¬ 
tails,  it  is  convenient  to  divide  the  tanks  into  two  general  categories: 
pressure  tanks  and  nonpressure  tanks.  Nonpressure  tanks  include  membrane, 
semi-membrane,  and  independent  prismatic  tanks.  These  tanks  are  con¬ 
structed  primarily  of  plane  surfaces  and  are  designed  for  low  vapor  pres¬ 
sures.  Design  for  low  internal  pressures,  particularly  near  the  tank  top 
where  the  liquid  head  is  a  minimum,  makes  these  tanks  susceptible  to  local¬ 
ized  sloshing  pressure  in  this  region.  Also,  prismatic  tanks  with  plane 
surfaces  are  more  likely  to  experience  high  liquid  impact  forces  during 
LNG  sloshing  than,  for  example,  pressure  tanks  with  curved  boundaries. 

In  contrast  to  nonpressure  tanks,  pressure  tanks  are  designed  for 
higher  internal  pressures,  approximately  3  to  10  times  higher  than  for  non¬ 
pressure  tanks.  In  addition,  as  reported  by  DnV,(24)  impulsive  sloshing 
pressures  are  unlikely  to  occur  in  pressure  tanks  because  of  their  curved 
boundaries,  and  the  limited  measurements  which  have  been  made  indicate  that' 
the  pressures  are  low.  Another  factor  is  that  pressure  tanks  which  are 
spherical  or  cylindrical  in  geometry  will  react  to  the  sloshing  pressures, 
even  localized  ones,  primarily  with  membrane  action  (as  opposed  to  bending 
action  in  tanks  with  plane  surfaces),  and  thus  pressure  tanks  are  less  sus¬ 
ceptible  to  damage  from  local  sloshing  pressure  if  they  do  occur. 
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Pressure  tanks*  structural  arrangements  are  much  simpler  than  tanks 
constructed  with  plane  surfaces.  The  structure  consists  principally  of  a 
smooth  wall  sphere  or  cylinder  of  relatively  heavy  gauge  (1/2  in.).  Attach¬ 
ments  to  the  tank  are  made  only  at  its  support  points  and  at  filling  loca¬ 
tions.  A  nonstructural  insulation  is  bonded  to  the  outer  surface  of  the 
tank  wall.  The  most  widely  used  LNG  pressure  tank  is  the  Moss-Rosenburg 
system.  This  system  uses  an  aluminum  spherical  tank  which  is  attached  to 
the  ship  by  a  cylindrical  skirt  at  the  equatorial  ring  (see  Figure  II-2b) . 

A  discussion  of  design  loads  for  this  tank  by  Glasfeld(23)  indicates  that 
sloshing  forces  influence  only  the  design  of  the  tank  support. 

As  indicated  above,  nonpressure  tanks  are  more  likely  than  pressure 
tanks  to  experience  high  sloshing  pressures  and  to  be  damaged  by  them. 

Since  the  structural  details  of  pressure  tanks  are  not  unique,  only  the 
more  complex  and  unique  aspects  of  I,NG  prismatic  tank  structures  will  be 
reviewed.  The  structures  of  two  types  of  membrane  tanks  and  one  independent 
prismatic  tank  are  described  in  the  following  paragraphs.  These  tanks  are 
representative  of  the  range  of  designs  of  nonpressure  type  tanks  which  are 
being  built  today. 

III. 3.1  Membrane  Tanks 


Gaz  Transport  Design.  In  the  Gaz  Transport  membrane  system,  invar 
membranes  are  supported  from  plywood  boxes  which  are  filled  with  Perlite 
to  provide  thermal  insulation.  Two  layers  of  boxes  are  used,  primary  and 
secondary.  Primary  boxes  are  on  the  inside  of  the  tank  adjacent  to  the  LNG 
and  support  the  primary  membrane.  Primary  boxes  are  separated  from  the  sec¬ 
ondary  boxes  by  the  secondary  membrane.  Both  secondary  and  primary  boxes 
are  directly  attached  to  the  ship's  hull,  but  compression  forces  on  the  pri¬ 
mary  boxes  are  transferred  to  the  hull  through  the  secondary  boxes. 

The  geometry  of  a  12-cell  box  is  given  in  Figure  III-18,  and  the 
arrangement  of  the  insulation  boxes  in  the  tank  is  shown  in  Figure  III-19. 
Attachments  of  the  boxes  to  the  hull  are  not  shown.  Primary  boxes  have 
both  a  9-cell  and  12-cell  construction.  Twelve-cell  boxes  are  used  in 
areas  where  maximum  overpressures  are  expected.  All  secondary  boxes  have 
nine  cells.  Secondary  boxes  bear  directly  against  the  inner  hull  of  the 
ship  on  the  bottom  and  sides  of  the  tank,  against  transverse  bulkheads 
(cofferdams)  at  the  tank  ends,  and  against  the  bottom  of  the  trunk  deck  at 
the  tank  top.  An  example  of  this  type  of  structure  is  shown  in  Figure 
III-20.  It  is  a  section  of  the  cofferdam  near  the  top  of  the  tanks.  Local 
structure  between  vertical  webs  is  plating  stiffened  by  bulb  angles. 

Sloshing  pressures  on  the  membrane  apply  loads  directly  to  the  covers 
of  the  primary  boxes.  These  forces  are  transferred  principally  through  the 
longitudinal  (parallel  to  the  box  axis)  stiffeners  in  the  primary  box  to 
the  transverse  stiffeners  in  the  secondary  boxes  at  the  stiffener  crossing 
points.  This  can  be  seen  from  Figure  III-19.  Transverse  stiffeners  in  the 
secondary  boxes  transfer  the  sloshing  forces  to  the  ship's  hull.  As  shown  in 
Figure  III-20,  each  stiffener  in  the  cofferdam  bulkhead  reacts  the  loading 
from  approximately  two  box  widths.  Because  the  insulation  boxes  are  very 
stiff  (loads  are  transferred  through  bending  of  the  covers  in  the  12-cell 
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FIGURE  1 1 1-20.  SECTION  OF  COFFERDAM  BULKHEAD  FOR  GAZ  TRANSPORT  TANKS 


box  and  through  compression  of  the  stiffeners) ,  the  sloshing  forces  will 
be  attenuated  very  little  at  the  bulkhead.  Structures  which  must  be  ana¬ 
lyzed  for  sloshing  pressures  include  the  primary  box  covers,  stiffeners 
in  the  boxes,  and  local  stiffened  plates  in  the  hull.  Because  sloshing 
pressures  are  localized,  occurring  only  near  the  static  liquid  level,  over¬ 
all  bending  of  the  ship's  bulkheads  is  not  a  problem. 

Technigaz  Membrane  System.  As  in  the  Gaz  Transport  system,  a  thin 
metal  membrane  is  attached  to  an  insulation  system,  which  in  turn  is  at¬ 
tached  directly  to  the  inner  hull  of  the  ship.  In  the  Technigaz  system 
the  secondary  membrane  is  a  continuous  plywood  barrier. 

A  section  through  the  membrane  and  insulation  of  the  Technigaz  Sys-' 
tem  is  shown  in  Figure  III-21.  The  insulation  and  membrane  support  in  this 
design  consists  of: 

(1)  An  end  grain  balsa  pad  which  is  adjacent  to  the  membrane, 

(2)  A  thick  sandwich  panel  built  with  plywood  surfaces  and  balsa 
core, 

(3)  Grounds  which  are  attached  to  the  ship's  hull  and  which 
support  compression  from  the  balsa-plywood  panels, 

(4)  Glass  wool  which  fills  the  voids  between  the  sandwich  panel 
and  the  ship's  hull  and  between  the  grounds. 

Structure  in  the  ship's  hull  which  supports  the  membrane  system  is  similar 
to  that  for  the  Gaz  Transport  design. 

Sloshing  pressures  on  the  membrane  bear  directly  on  the  balsa  pad. 
These  forces  are  transferred  into  the  balsa-plywood  panels,  producing  bend¬ 
ing  in  the  panels  between  the  grounds  and  shear  reactions  at  the  grounds. 
The  grounds  transfer  the  shear  reactions  in  the  panels  to  the  ship's  hull. 
Compression  in  the  balsa  pad,  bending  and  shear  in  the  panels,  and  local 
deformation  of  the  ship's  hull  must  be  checked  for  adequate  strength  to 
react  to  the  sloshing  pressures. 

III. 3. 2  CONCH  Independent  Tank 

The  CONCH  LNG  tank  system  is  designed  to  contain  the  static  and  dy¬ 
namic  pressures  produced  by  the  LNG  independent  from  the  ship's  hull.  It 
is  keyed  to  the  hull  to  prevent  translation  and  rotation  relative  to  the 
ship.  Glass  wool  fills  the  void  space  between  the  tank  and  the  ship  inner 
hull  (between  the  supporting  blocks) .  Tanks  are  subdivided  into  approxi¬ 
mately  four  equal  volumes  by  centerline  and  swash  bulkheads.  Construction 
of  the  tanks  is  similar  to  that  of  a  ship.  Bulkheads  and  tank  walls  are 
constructed  of  flat  plates  which  are  reinforced  by  stringers  and  girders. 

A  typical  horizontal  girder  (horizontal  section  through  one  of  the  tanks) 
is  shown  in  Figure  III-22.  The  weakest  parts  of  the  tank  appear  to  be  the 
longitudinal  and  swash  bulkheads,  the  sides  of  the  tank  near  the  top,  and 
the  tank  top.  These  parts  of  the  tank  have  0. 50- in. -wall  thickness  and  the 


FIGURE  III -22 .  TYPICAL  HORIZONTAL  GIRDER  -  CONCH  TANK 


widest  frame  spacing.  Figures  III-23  and  III-24  show  typical  vertical  sec¬ 
tions  taken  in  the  longitudinal  and  transverse  directions,  respectively. 

Both  the  frames  and  stiffened  panels  between  the  frames  of  the  tank 
must  be  checked  for  damage  to  sloshing  forces.  This  is  true  for  the  tank 
walls,  bulkheads,  and  tank  top.  In  addition,  tank  supports  must  be  analyzed 
for  the  reactions  produced  by  the  resultant  of  the  sloshing  forces. 
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23.  TYPICAL  LONGITUDINAL  GIRDER 
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IV.  TASK  2  -  EXPERIMENTAL  STUDIES 


IV. 1  Experimental  Study  Objectives 


As  a  result  of  the  Task  1  review  of  previous  experimental  work,  the 
following  experiments  were  performed. 

o  Combined  degree  of  freedom  excitation  (pitch  and  heave, 
surge  and  heave) 

o  Impulsive  pressure-time  history  measurements 


Dynamic  loads  test  on  full"scale  LNG  tank  structural 
segments 


Previous  studies  provide  loads  data  for  a  wide  range  of  tank  geome¬ 
tries,  fill  depths,  and  excitation  conditions.  Previous  studies  were  con¬ 
ducted  by  exciting  the  model  tank  in  only  one  degree  of  freedom  (usually 
pitch,  roll,  or  surge)  and  recording  loads  at  resonant  sloshing  conditions! 
To  determine  if  combined,  simultaneous  excitation  in  more  than  one  degree 
of  freedom  would  increase  the  sloshing  pressures,  additional  two  degree  of 
freedom  experiments  were  performed. 

The  impact  pressure-time  history  is  of  importance  to  the  prismatic 
tank  wall  structural  design.  Previous  studies  determined  peak  pressure 
magnitudes  and  distributions,  but  little  data  were  available  on  pressure¬ 
time  histories.  Additional  experiments  were  performed  to  establish 
P(T)  vs  T  characteristics. 

To  substantiate  the  analytical  procedures  developed  in  this  study 
and  to  evaluate  the  response  of  LNG  tank  wall  structures  to  dynamic  slosh 
pressures,  typical  full-scale  dynamic  pressures  (predicted  from  the  model 
tests)  were  applied  to  a  representative  section  of  a  full-scale  LNG  mem¬ 
brane  tank  structure  in  a  set  of  laboratory  tests. 


Experimental  Facilities 


IV. 2.1  Model  Tank  Sloshins 


jaratus 


The  laboratory  work  on  this  project  involved  sloshing  liquid  in 
scale  model  ship  tanks  and  measuring  slosh-induced  dynamic  pressures.  The 
laboratory  facilities  at  SwRI  include  two  slosh  rigs  designed  for  this  pur¬ 
pose.  One  machine  is  capable  of  producing  an  angular  motion  (to  simulate 
pitching  or  rolling) ,  and  the  other  can  provide  independent  and  simultane¬ 
ous  horizontal  and  vertical  motions  (to  simulate  surge  or  sway,  and  heave). 
Appendix  A  provides  a  detailed  description  of  both  these  sloshing  systems. 


IV. 2.2  Sloshing  Pressure  Measurement  System 


Sloshing  dynamic  pressures  were  recorded  with  Columbia  Research 
Labs'  Model  100-P  piezoelectric  pressure  transducers  in  combination  with 
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Model  4101  amplifiers.  These  transducer/araplifier  combinations  have  a  range 
of  0.005  to  4000  pounds  per  square  inch,  with  a  frequency  response  from  0.1 
to  10,000  Hz.  Transducers  were  isolated  from  the  walls  of  the  tank  by  mount¬ 
ing  them  in  a  plug  of  silastic  material.  The  silastic  eliminated  the  possi¬ 
bility  of  the  transducer  casing  being  loaded  by  flexing  of  the  walls  of  the 
model  tank. 

The  pressure  signals  were  input  to  a  light-beam  galvanometer-type 
oscillograph.  The  oscillograph  produced  a  permanent  record  of  the  pressure 
signals  on  oscillograph  paper.  The  oscillograph  is  a  Honeywell  Model  1858- 
0790  G  Visicorder  oscillograph  with  a  frequency  response  from  0  to  3000  Hz 
with  +0%  to  -5%  error.  A  Hewlett-Packard  Model  9864A  digitizer  in  combi¬ 
nation  with  a  Hewlett-Packard  Model  9830A  computer  was  used  to  digitize  the 
pressure  signals  from  the  oscillograph  traces.  The  digitized  pressure  val¬ 
ues  were  stored  on  magnetic  disk  files  for  final  data  reduction. 

IV. 2. 3  Model  Ship  Tank 

The  model  tank  used  for  all  laboratory  tests  is  a  1/50  scale  replica 
of  a  prismatic  tank  from  a  125,000  m^  LNG  carrier.  The  tank  geometry  is 
typical  of  most  prismatic  tanks.  The  model  tank  walls  are  made  of  plexiglass 
with  sufficient  rigidity  such  that  tank  wall  response  was  not  a  factor  in  the 
pressure  measurements.  Figure  IV-1  shows  the  tank  and  all  pressure  trans¬ 
ducer  locations.  Table  IV-1  lists  the  model  scale  and  corresponding  full” 
scale  tank  dimensions.  The  pitch  axes  for  the  25%  and  75%  tank  filling 
levels  investigated  are  denoted  in  Figure  IV-1  as  zp^  and  zP2»  respectively. 

For  all  tests,  the  model  test  liquid  was  water. 

IV. 3  Combined  Degree  of  Freedom  Model  Tests 

As  part  of  Task  2,  the  effects  of  combined  tank  motions  on  slosh 
pressures  were  investigated.  The  purpose  of  this  work  was  to  determine  if 
slosh-induced  pressures  were  augmented  by  combined  motions. 

Experimental  tank  motions  included  pitching  plus  heaving  and  surging 
plus  heaving.  The  results  for  these  motions  were  compared  to  the  results 
for  the  individual  motions  of  pitching,  surging,  and  heaving.  Low  (25%) 
and  high  (75%)  tank  filling  levels  were  used  for  each  case.  All  laboratory 
test  motions  are  summarized  in  Table  IV-2.  The  detailed  test  results  are 
given  in  Tables  TV-3  through  IV-7,  with  the  data  summarized  in  Tables  IV-8 
through  IV-11.  The  excitation  amplitudes  were  chosen  to  represent  typical 
extreme  values  that  a  ship  tank  might  experience  while  in  service.  Each 
test  was  run  for  200  cycles  at  resonant  sloshing  conditions,  while  the  tank 
was  undergoing  a  harmonic  excitation  motion  at  resonance. 

IV. 3. I  Pitching  Plus  Heaving  Test 


The  angular  motion  sloshing  facility  was  used  to  provide  the  pitch¬ 
ing  plus  heaving  motion.  The  iscale  model  ship  tank  was  placed  approximately 
two  tank  lengths  out  from  the  pitch  axis  on  a  cantilevered  beam.  This  lo¬ 
cation  represents  the  typical  placement  of  the  fore  or  aft  tank  in  an  LNG 
carrier  with  its  pitching  axis  at  the  midship.  Figure  IV-2  shows  the  loca¬ 
tion  of  the  model  tank  with  respect  to  the  sloshing  rig  axis. 
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O  ~  Denotes  pressure  transducer  location 


FIGURE  IV-1.  SCALE  MODEL  TANK  DIMENSIONS  AND  PRESSURE  TRANSDUCER 
LOCATIONS  FOR  1/50  SCALE  PRISMATIC  TANK  FROM  A  125,000  m3  SHIP 


TABLE  IV-1.  TANK  DIMENSIONS  AND  PRESSURE 
TRANSDUCER  LOCATIONS 


Dimension 

Model  (ft) 

Full  Scale  (ft) 

l 

2.23 

111.55 

b 

2.44 

122.04 

H 

1.57 

j 

78.51 

X1 

0.219 

10.94 

X2 

1.22 

60.94 

X3 

2.22 

110.94 

x4 

0.333 

16.67 

X5 

0.885 

44.27 

*1 

0.104 

5.21 

21 

0.146 

7.29 

Z2 

0.417 

20.83 

Z3 

0.792 

39.58 

z4 

1.23 

61.46 

Z5 

1.44 

71.88 

ZP1 

0.196 

9.81 

ZP2 

0.589 

29.44 

Conversion  Factor:  1  foot  =  0.3048  meters. 
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TABLE  I V-2 .  TEST  CONDITIONS  FOR  COMBINED 
MOTION  LABORATORY  TESTS 


Tank  Motion 

Tank  Amplitude 

Tank 

Filling 

Level 

Nondimensional 
Resonant  Slosh  Period 

(T  /g7z) 

Pitching 

4  -  +4° 

25% 

5.170 

Pitching 

$  *  +4° 

75% 

3.566 

Surging 

f  -  0.019 

25% 

4.869 

Surging 

f  -  0.019 

75% 

3.663 

Surging 

|  -  0.038 

25% 

4.969 

Surging 

|  -  0.038 

75% 

3.663 

Pitching  +  Heaving 

$  -+40  and  ^*0.417 

25% 

<4.748 

Pitching  +  Heaving 

$  «»+40  and  0.417 

75% 

7.235 

(2nd  Harmonic) 

Surging  +  Heaving 

0.038  and  0.053 

25% 

4.869 

Surging  +  Heaving 

0.038  and  0.053 

75% 

3.663 

See  Figures  IV- 1  and  IV-4  for  definitions  of  x,  y,  i,  H,  and  <p 
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NOTE:  See  Equation  III-l  for  definition  of  KP 
Conversion  Factor:  1  psi  -  6. 895 (10) “3  N/nna2. 
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NOTE:  See  Equation  III-l  for  definition  of  KP 
Conversion  Factor:  1  psi  »  6.895(10)“^  N/nm^. 
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NOTE:  See  Equation  III-l  for  definition  of  KP 
Conversion  Factor:  1  psi  «  6. 895 (10) N/am2. 


TABLE  IV— 7.  TEST  PRESSURE  DATA  FOR  SURGING 
PLUS  HEAVING  TANK  MOTION 
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NOTE:  See  Equation  III-l  for  definition  of  KP 
Conversion  Factor:  1  psi  *  6. 895 (10 )~3  N/xnm2. 


TABLE  IV-8.  AVERAGE  KP  VALUES  FOR  A  1/50  SCALE 
PRISMATIC  TANK  FROM  A  125,000  m3  LNG  CARRIER 
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TABLE  IV-11.  MAXIMUM  FULL  SCALE  PRESSURE  VALUES  FOR  A 
1/50  SCALE  PRISMATIC  TANK  FROM  A  125,000  m3  LNG  CARRIER 


Due  to  the  design  of  the  test  rig,  the  pitching  and  heaving  motions 
were  simultaneous,  but  not  independent  of  each  other.  The  two  motions  were 
at  the  same  frequency  and  in  phase  for  these  tests.  During  testing,  it  was 
noticed  that  the  fundamental  excitation  frequency  for  a  75%  full  tank  under¬ 
going  pitching  plus  heaving  did  not  produce  substantial  dynamic  pressures, 
so  the  frequency  was  increased  to  the  first  harmonic.  At  this  frequency, 
substantial  sloshing  pressures  occurred  and  were  recorded  for  inclusion  in 
the  test  results.  The  theoretical  (Equation  III-6)  and  experimental  exci¬ 
tation  frequencies  are  shown  in  Figure  IV-3. 

The  recorded  pressure  values  for  pitching  plus  heaving  and  only 
pitching  are  contained  in  Tables  IV-3  and  IV-4.  No  values  are  found  in 
the  tables  for  heaving  only  since  the  resonant  frequency  for  a  heaving  mo¬ 
tion  is  twice  the  natural  frequency  of  the  liquid  in  the  tank.  The  heave 
resonant  frequency  is  outside  the  range  of  interest  for  ship  motions. 
Therefore,  no  attempt  was  made  to  excite  the  tank  into  resonance  in  heave 
and  record  pressures. 

The  pressures  presented  in  Tables  IV-3  and  IV-4  for  pitching  and 
pitching  plus  heaving  indicate  there  was  no  augmenting  of  the  pressures 
using  the  combined  motion.  In  addition,  many  test  pressures  were  somewhat 
higher  for  pitching  only  motion. 

IV. 3. 2  Surging  Plus  Heaving  Test 

The  translational  motion  sloshing  facility  was  utilized  for  the 
surging  plus  heaving  tests.  As  with  the  pitching  plus  heaving  test,  the 
surging  plus  heaving  test  results  are  compared  to  the  results  of  the  indi¬ 
vidual  excitation  motion  (surging).  Figure  IV-4  illustrates  the  surging 
and  heaving  motions. 

Two  surging  amplitudes  were  examined.  These  are  typical  extreme 
valves  that  a  full-scale  tank  might  experience  during  operation.  To  deter¬ 
mine  the  worst  case  combined  surging  plus  heaving  condition,  the  tank  was 
oscillated  with  the  surging  frequency  set  at  resonance,  and  the  heaving 
frequency  swept  through  a  range  of  values  from  0  to  5  Hz  (model  scale)  or 
0  to  0.707  Hz  (full  scale).  Next,  the  tank  was  oscillated  with  the  heav¬ 
ing  motion  set  at  its  theoretical  resonant  frequency  and  the  surging  exci¬ 
tation  frequency  swept  through  a  range  of  values  [0  to  5  Hz  (model  scale) 
or  0  to  0.707  Hz  (full  scale)].  The  worst  case  was  chosen  as  the  one  that 
created  the  most  violent  sloshing  (by  visual  observation)  and  the  highest 
recorded  pressures. 

The  results  of  these  tests  indicate  that  the  worst  case  sloshing 
occurred  when  the  surging  and  heaving  excitation  frequencies  were  equal  to 
the  resonant  frequency  of  the  liquid  for  a  surging  only  motion.  The  test 
results  were  recorded  with  both  the  surging  and  heaving  motions  at  the  same 
frequency  and  in  phase.  The  tank  displacements  in  the  surging  and  heaving 
directions  were  also  equal  in  magnitude  (x  =  y  =  1  in.®  2.54  cm),  with 
(x/2)  =  0.038  and  (y/H)  =  0.053. 
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10 


Tank  Filling  Level 
(Percent) 


\ 


Tr  ■  resonant  sloshing  period 

g  »  acceleration  of  gravity 

i  ■  tank  length  in  direction  of  motion 

FIGURE  IV-3.  EXPERIMENTAL  AND  THEORETICAL  NONDIMENSIONAL  RESONANT 
SLOSH  PERIOD  VERSUS  TANK  FILLING  LEVEL 
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x  *  surge  motion  displacement 
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i.  ■  tank  length 
H  ■  tank  height 

FIGURE  IV-4.  MOTION  DEFINITIONS  FOR  SURGE  AND  HEAVE. 
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The  test  results  for  surging  plus  heaving  and  surging  only  are  found 
in  Tables  IV-5,  IV-6,  and  IV-7.  The  results  do  not  indicate  any  signifi¬ 
cant  augmenting  of  the  sloshing  pressures  caused  by  the  combined  tank  mo¬ 
tions  of  surging  plus  heaving,  Visual  observations  also  suggested  that  the 
sloshing  for  the  surging  plus  heaving  was  comparable  in  magnitude  to  that 
for  only  surging. 

Based  on  the  results  of  these  laboratory  tests,  there  does  not  ap¬ 
pear  to  be  an  increase  in  the  slosh-induced  dynamic  impact  pressures  caused 
by  multi-degree  of  freedom  tank  motions.  However,  it  is  realized  that 
these  tests  are  limited  in  that  simultaneous  excitation  of  all  six  degrees 
of  freedom  could  not  be  investigated.  In  spite  of  this  limitation,  the 
results  do  provide  insight  into  the  effects  of  combined  degrees  of  freedom 
excitation  and  indicate  that  worst  case  loads  can  be  obtained  by  resonant 
sloshing  in  single  degree  of  freedom  tests.  The  excellent  full -scale/model- 
scale  pressure  agreement  for  the  0B0  tank  (discussed  in  Section  III-2)  where 
model  tests  were  included  in  roll  only  further  support  this  conclusion. 

IV. A  Dynamic  Pressure-Time  Histories 

Pressure-time  histories  were  measured  in  model  scale  to  establish 
information  important  to  the  analysis  of  the  tank's  structural  response  to 
dynamic  loading  caused  by  cargo  sloshing. 

A  typical  pressure-time  history  (nondimensional)  is  shown  in  Figure 
IV-5.  Fer  analytical  purposes,  this  pressure  trace  will  be  approximated 
by  straight  lines  connecting  points  1,  2,  3,  A,  and  5.  The  most  important 
elements  of  the  pressure  trace  are  (1)  the  time  required  to  reach  the  peak 
pressure  (point  2) ,  (2)  the  magnitude  of  the  peak  pressure  at  point  2,  (3) 
the  duration  of  the  pressure  spike  (point  1  to  point  3),  (A)  the  magnitude 
of  the  tailing  pressure  between  points  3  and  A,  (5)  the  duration  of  the 
tailing  pressure  betweei.  points  3  and  A,  and  (6)  the  integrated  pressure 
value  over  the  entire  pressure  cycle.  By  knowing  these  characteristics  of 
the  pressure-time  history,  one  can  evaluate  the  structural  response  of  the 
tank  wall  to  an  impact  pressure. 

Wall  impact  pres.,ure-time  histories  for  both  a  25%  and  a  75%  full 
tank  were  recorded.  The  tank  excitation  motion  was  harmonic  pitching. 

Two  hundred  resonant  sloshing  cycles  were  recorded  for  each  test.  For  the 
25%  full  test,  pressures  were  measured  at  transducer  locations  1  through  9 
(Figure  IV-1) .  Pressures  at  locations  7  through  15  were  recorded  for  the 
75%  full  case. 

The  test  results  include  two  forms  of  the  pressure-time  history. 
First,  fer  each  transducer  location,  there  is  a  composite  worst  case  pres¬ 
sure-time  history.  This  profile  includes  the  shortest  pressure  rise  time 
(T2) ,  the  longest  time  durat.  ..s  fer  T3,  T4,  and  T5,  and  the  highest  pres¬ 
sures  for  all  points  (P2,  P3,  P4>  and  P5)  recorded  during  any  of  tire  200 
sloshing  cycles.  This  trace  represents  the  worst  case  values  (with  respect 
to  the  loading  on  a  tanv  wall  structure)  tuat  the  particular  transducer 
location  experienced  during  the  200  cycles  of  resonant  sloshing.  In  addi¬ 
tion,  a  pressure-time  history  showing  the  average  values  (for  200  resonant 
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The  Trapezoidal  Rule  using  the  five  numbered  points  Is  used 
to  determine  the  area  under  the  curve  (which  Is  the  value  of 
the  pressure  integrated  over  the  entire  period) . 


FIGURE  IV-5.  TYPICAL  PRESSURE-TIME  HISTORY  FOR  SLOSH-INDUCED  IMPACT 
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sloshing  cycles)  for  the  points  on  the  curve  are  plotted.  Also,  the  inte¬ 
grated  pressure  values  versus  the  pressure  duration  (T5)  a.e  presented. 

An  effort  was  made  to  determine  if  any  discernible  correlations 
exist  between  the  peak  pressure  (P2) ,  the  pressure  rise  time  (T2) ,  and  the 
pressure  spike  duration  (T3) .  Plots  showing  peak  pressure  versus  pressure 
rise  time,  peak  pressure  versus  pressure  spike  duration,  and  pressure  spike 
duration  versus  pressure  rise  time  are  included  in  the  results. 

The  nondimensional  pressure-time  history  test  results  for  represen¬ 
tative  transducer  locations  1  and  1A  are  found  in  Tables  IV-12  and  IV-13, 
and  Figures  IV-6  through  IV-15.  Similar  data  for  the  other  transducer  lo¬ 
cations  are  presented  in  Appendix  B.  Maximum,  average,  and  minimum  pressure 
and  time  values  are  included  in  the  tables.  It  should  be  noted  that  loca¬ 
tions  7,  8,  and  9  experienced  virtually  no  pressures  when  the  tank  was  25% 
full.  Thus,  there  are  no  tables  for  these  locations.  The  same  is  true  of 
location  15  for  the  75%  full  tank. 

The  pressure  data  obtained  at  transducer  locations  1,  A,  and  7  are 
representative  results  for  the  test  case  with  the  tank  25%  full.  Dynamic 
pressures  were  measured  at  locations  1  through  6.  The  pressure-time  his¬ 
tories  at  transducer  1  were  similar  in  magnitude  and  shape  to  those  at 
locations  2  and  3.  The  pressure  traces  at  transducer  location  A  were  simi¬ 
lar  to  those  at  locations  5  and  6. 

The  pressure-time  history  for  location  1  for  the  25%  full  tank  is 
found  in  Figure  IV-6.  The  integrated  pressure  value  for  this  transducer 
is  shown  in  Figure  IV-7.  Plots  exemplifying  the  relationships  between  the 
pressure  rise  time  (T2),  pressure  spike  duration  (T3),  and  the  peak  pres¬ 
sure  (P2)  for  this  location  are  shown  in  Figures  IV-8  through  IV-10.  It 
should  be  noted  that  several  of  the  plots  have  data  points  clustered  in 
small  areas.  This  was  caused  by  round-off  limitations  in  the  computation 
methodology  that  was  used  for  data  reduction. 

A  group  of  data  similar  to  that  just  described  for  the  test  case 
with  a  25%  full  tank  was  also  compiled  for  the  case  with  the  tank  75%  full. 
Data  from  pressure  transducer  locations  7,  11,  13,  and  1A  were  chosen  as 
representative  for  the  75%  full  tank.  The  pressure-time  histories  at  lo¬ 
cation  7  are  similar  to  those  at  locations  8  and  9.  Likewise,  the  pressure¬ 
time  histories  at  transducer  11  are  comparable  to  those  at  10  and  12  (with 
the  exception  that  the  magnitudes  at  location  10  are  somewhat  smaller). 

Graphs  showing  pressure-time  histories,  integrated  pressures,  and 
relationships  between  the  pressure  rise  time,  the  pressure  spike  duration, 
and  the  peak  pressure  for  transducer  1A  ate  found  in  Figures  IV-11  through 
IV-15.  Locations  7,  8,  and  9  were  below  the  static  liquid  level,  and,  as 
a  result,  these  locations  did  not  experience  significant  dynamic  loadings. 
The  pressure  rise  times  were  relatively  slow  and  the  peak  pressure  small. 
Pressures  of  a  more  dynamic  nature  were  measured  at  locations  10  through  1A. 
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FIGURE  IV-6.  NONDIMENS I ONAL  PRESSURE-TIME  FIGURE  IV-7.  INTEGRATED  NONDIMENSIONAL  PRESSURE 
HISTORY  VALUES  FOR  200  RESONANT  SLOSHING  VALUES  FOR  200  RESONANT  SLOSHING  CYCLES  AT 

CYCLES  AT  TRANSDUCER  LOCATION  1  FOR  A  25%  TRANSDUCER  LOCATION  1  FOR  A  25%  FULL  TANK 

FULL  TANK 


ICICKT/KP 


RESOWW  SLOSHING  CONDITIONS 


RESOKXIT  SLCSHIIS  COI©ITIClS 


IMPOLSC  RISC  TIKE/T2 
(KOI©  IN.) 

FIGURE  IV-8.  NONDIMENSIONAL 
PRESSURE  VS  IMPULSE  RISE  TIME 
FOR  200  RESONANT  SLOSHING  CYCLES 
AT  TRANSDUCER  LOCATION  1  FOR  A 
25%  FULL  TANK 


IHPUL5E  MJR3T1CS/T3 
(KOI©  IK.) 

FIGURE  IV-9.  NONDIMENSIONAL  PRESSURE 
VS  IMPULSE  DURATION  FOR  200 
RESONANT  SLOSHING  CYCLES  AT  TRANS¬ 
DUCER  LOCATION  1  FOR  A  25%.  FULL 
TANK 


KsaiRiT  straws  conditiris 

PITCH 

TffiNS.LOC.IIO.  I 
10.  OF  CTO.ES:  230 
FILL  Ol?TH;  2B 


KSBIfiNT  SLOSHIKS  CONDITIONS 
PITCH 
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FILL  DEPTH'-  TK 


i'll  8.20  8.30  B.HB  l.i 

IMPULSE  RISE  T1RE/T2 
(NONDIH.) 

FIGURE  IV-10.  NONDIMENSIONAL  IM¬ 
PULSE  DURATION  VS  IMPULSE  RISE 
TIME  FOR  200  RESONANT  SLOSHING 
CYCLES  AT  TRANSDUCER  LOCATION  1 
FOR  A  25%  FULL  TANK 
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FIGURE  IV- 11.  NONDIMENSIONAL  PRESSURE¬ 
TIME  HISTORY  VALUES  FOR  200  RESONANT 
SLOSHING  CYCLES  AT  TRANSDUCER  LOCATION 
14  FOR  A  75%  FULL  TANK 
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FIGURE  IV-12.  INTEGRATED  NON- 
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200  RESONANT  SLOSHING  CYCLES  AT 
TRANSDUCER  LOCATION  14  FOR  A  75* 
FULL  TANK 
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FIGURE  IV-13.  NONDIMENSIONAL 
PRESSURE  VS  IMPULSE  RISE  TIME 
FOR  200  RESONANT  SLOSHING  CYCLES 
AT  TRANSDUCER  LOCATION  14  FOR  A 
75*  FULL  TANK 
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FIGURE  IV-14.  NONDIMENSIONAL 
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AT  TRANSDUCER  LOCATION  14  FOR  A 
75*  FULL  TANK 
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FIGURE  IV-15.  NONDIMENSIONAL  IMPULSE 
DURATION  VS  IMPULSE  RISE  TIME  FOR 
200  RESONANT  SLOSHING  CYCLES  AT 
TRANSDUCER  LOCATION  14  FOR  A  75* 
FULL  TANK 
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IV. 5  Dynamic  Load  Simulator  for  Plywood  Insulation  Box  Tests 


A  series  of  strength  tests  were  performed  on  representative  plywood 
insulation  boxes  (Figure  III-18)  from  a  typical  membrane  tank  insulation 
system.  The  purpose  was  to  (1)  provide  structural  response  data  for  com¬ 
parison  to  analytical  predictions,  (2)  to  evaluate  the  box  response  to  var¬ 
ious  projected  full  scale  loading  conditions,  and  (3)  to  determine  pressure 
levels  a  box  could  withstand  without  sustaining  damage.  The  loading  con¬ 
ditions  were  designed  to  simulate  what  a  structure  might  encounter  in  actual 
operation.  Tests  included  both  static  and  dynamic  loadings. 

Because  of  the  unusual  test  conditions,  a  special  test  fixture  was 
fabricated.  It  was  necessary  to  uniformly  load  the  entire  cover  of  the 
test  box.  In  addition,  a  typical  full  -scale  pressure-time  history  (Figure 
IV-5)  had  to  be  produced  on  the  test  structure. 

The  test  system  consisted  of  a  chamber  of  water  located  on  top  of 
the  box  cover  with  a  rubber  diaphragm  between  the  two.  The  box  cover  was 
loaded  by  pressurizing  the  water  chamber.  This  created  a  relatively  uni¬ 
form  load  over  the  cover  since  the  water  deformed  to  the  shape  of  the  cover 
as  it  deflected  under  loading.  For  the  dynamic  pressure  tests,  three  fast¬ 
opening  valves  located  between  the  pressure  source  and  the  water  chamber 
were  used  to  shape  the  pressure-time  history.  A  schematic  of  the  test 
apparatus  is  presented  in  Figure  IV-16. 

The  following  paragraphs  describe  the  function  of  each  component 
shown  in  Figure  IV-16.  The  pressure  supply  was  a  tank  filled  with  nitrogen 
gas.  At  the  exit  port  to  the  pressure  tank  was  a  hydraulically-operated 
cock  valve.  Between  the  cock  valve  and  the  water  chamber  was  a  fast-open¬ 
ing,  hydraulically  operated  plunger  valve.  When  this  valve  and  the  cock 
valve  were  both  opened,  the  high  pressure  nitrogen  gas  pressurized  the  wa¬ 
ter  volume.  A  solenoid  valve  located  between  the  cock  valve  and  the  plunger 
valve  was  used  as  a  relief  valve.  The  pressure-time  history  for  a  dynamic 
pressure  test  was  controlled  by  the  timing  of  the  opening  and  closing  of 
these  three  valves.  The  apparatus  was  capable  of  creating  dynamic  pressure 
rise  times  as  short  as  20  milliseconds.  Based  on  the  pressure-time  history 
experiments,  this  is  a  typical  rise  time  that  a  plywood  insulation  box  might 
experience  while  in  service.  Peak  dynamic  pressures  of  about  200  psi  were 
attainable.  Static  pressures  on  the  order  of  250  psi  were  achieved  with 
this  particular  test  rig.  These  pressures  represent  worst  case  pressures 
that  could  be  encountered  in  full  scale -resonant  sloshing. 

The  supporting  structure  for  the  plywood  test  box  (called  the  pri¬ 
mary  box)  consisted  of  1 /2-in. -thick,  marine-grade  fir  plywood  boxes  (called 
secondary  boxes) .  These  boxes  simulated  the  supporting  structure  of  the 
primary  box  in  the  ship  tank.  Figure  III-19  shows  the  supporting  configur¬ 
ation  as  it  is  in  the  tank.  The  secondary  boxes  used  in  the  loading  tests 
had  stiffener  locations  similar  to  those  shown  in  Figure  III-19.  This  al¬ 
lowed  the  load  transfer  through  the  stiffeners  of  the  primary  test  box  to 
closely  approximate  the  actual  full-scale  conditions.  The  secondary  boxes 
in  the  cest  apparatus  rested  on  a  rigid  surface. 
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FIGURE  IV-16. 


SCHEMATIC  DIAGRAM  OF  PLYWOOD  INSULATION  BOX  DYNAMIC  LOADER 


Instrumentation  for  the  loading  tests  included  pressure  measuring 
devices  to  monitor  the  pressure  on  the  box  cover  and  displacement  trans¬ 
ducers  to  measure  the  deflection  of  the  cover  during  loading.  For  the  dy¬ 
namic  loading  tests,  PCB  piezoelectric  pressure  transducers  (Model  102A15) 
were  used  uo  measure  the  pressures  in  the  water  chamber.  These  transducers 
have  an  operating  range  of  0.01  to  500  psi  with  a  frequency  range  of  0.5  to 
50,000  Hz.  Three  transducers  were  used  to  monitor  the  uniformity  of  the 
pressure  in  the  chamber.  Transducers  were  located  at  the  midpoint  of  both 
end  walls  of  the  chamber  and  at  the  midpoint  of  one  side  wall.  All  three 
transducers  were  located  approximately  1-1/2  in.  above  the  surface  of  the 
test  box  cover.  The  general  locations  of  the  three  pressure  transducers 
and  the  pressure  gage  are  pictured  in  Figure  IV-16.  Figure  IV-17  shows 
typical  pressure  traces  for  the  three  transducers  for  a  dynamic  pressure 
test.  The  top  photograph  shows  the  pressure-time  histories  for  a  test  that 
did  not  experience  a  box  failure.  The  bottom  picture  illustrates  the  pres¬ 
sure  waveform  when  a  failure  occurred.  These  pressure  signals  were  re¬ 
corded  on  a  storage  oscilloscope,  as  were  all  of  the  dynamic  pressure  mea¬ 
surements. 

One  important  aspect  of  these  plywood  box  tests  was  the  determina¬ 
tion  of  the  box  cover  response  during  loading.  The  stiffener  configuration 
divided  the  boxes  into  12  cells.  The  cover  deflections  were  measured  at 
the  centers  of  the  unsupported  cell  areas.  Figure  IV-18  shows  the  cell 
divisions  for  a  typical  12-celled  box  and  the  locations  of  the  displacement 
measurements.  For  static  loading  tests,  spring-loaded,  linear  displacement 
transducers  were  utilized  for  the  measurements.  For  dynamic  loadings,  a 
Kaman  (Model  KD2300-25)  inductance- type  displacement  transducer  was  used. 

The  performance  of  the  two  types  of  transducers  correlated  well  as  evidenced 
by  the  results  presented  in  Figure  IV-19.  'phis  figure  shows  that  the 
cover  response  was  similar  for  static  and  dynamic  loading  conditions. 

The  cover  deflection  measurements  are  presented  in  Tables  IV-14 
through  IV-21.  In  addition  to  these  results.  Table  IV-22  summarizes  the 
tests  and  results  for  all  boxes.  Some  boxes  were  modified  to  study  the  in¬ 
fluence  of  various  parameters  on  box  strength. 

Three  types  of  failures  were  observed  during  the  box  loading  tests. 
First,  some  covers  experienced  shear  failures.  Usually  the  shear  failure 
would  occur  at  a  stiffener  support.  Second,  other  covers  experienced  bend¬ 
ing  failures  along  the  centerline  of  the  unsupported  area  of  an  individual 
cell.  The  third  type  of  failure  was  a  buckling  failure  of  a  stiffener. 

This  was  observed  in  only  one  test.  Figures  IV-20,  IV-21,  and  IV-22  show 
typical  shear,  bending,  and  stiffener  failures,  respectively. 

The  test  results  indicated  that  the  original  12-celled  boxes  could 
withstand  a  pressure  of  about  120  psi  without  suffering  catastrophic  fail¬ 
ure.  The  range  of  failure  pressures  was  120  to  158.6  psi  for  both  static 
and  dynamic  loading.  Also,  test  results  showed  no  appreciable  difference 
in  failure  pressures  for  static  or  dynamic  pressures.  A  typical  dynamic 
loading  is  shown  in  Figure  IV-17.  The  pressure  rise  time  is  about  20  ms,  the 
minimum  rise  time  obtainable  with  the  loading  device.  These  observations  are 
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Time  Scale:  50  msec/div 


Time  Scale:  50  msec/div 


FIGURE  IV-17.  TYPICAL  PRESSURE-TIME 
PLYWOOD  INSULATION  BOXES 


Pressure  Scale:  50  psi/div 
(0.345  N/mmVdiv) 

(Typical  pressure  trace  when 
box  failure  does  not  occur.) 


I  Pressure  Scale:  100  psi/div 
(0.690  N/n*n2/div) 

(Typical  pressure  trace 
showing  box  failure.) 


HISTORIES  FOR  DYNAMIC  LOADING  TESTS  ON 
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FIGURE  IV-21.  BENDING  FAILURE  OF  PLYWOC 
(TEST  BOX  NO.  8) 


INSULATION  BOX  COVER 
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TABLE  IV-14,  PLYWOOD  COVER  DEFLECTION  MEASUREMENTS 
FOR  TEST  EOX  NUMBER  1  (Static  Loading  Test) 


Pressure 

(psi) 

Cover  Deflection  (Inches) 

Cell  No.  1 

Cell  No.  3 

Cell  No.  10 

Cell  No.  12 

0 

0.00 

0.00 

0.00 

0.00 

35 

0.00 

0.00 

0.01 

0.00 

50 

0.02 

0.00 

0.04 

0.01 

60 

0.04 

0.00 

0.05 

0.03 

65 

0.04 

0.00 

0.05 

0.03 

70 

0.04 

0.00 

0.06 

0.03 

80 

0.05 

0.00 

0.06 

0.04 

90 

0.07 

0.05 

0.07 

0.04 

100 

0.08 

0.08 

0.09 

0.06 

110 

0.10 

0.08 

0.13 

0.08 

120 

0.20 

0.16 

0.17 

0.18 

Cell  number  11  failed  in  shear  at  130  psi  (0.896  N/mra2). 

Conversion  Factors;  1  psi  •  6.895(10)"^  N/wt2 
1  in.  «  25.4  mra 

TABLE  IV-15.  PLYWOOD  COVER  DEFLECTION  MEASUREMENTS 
FOR  TEST  t'OX  NUMBER  2  (Static  Loading  Test) 


Pressure 

(psi) 

Cover  reflection  (inches)  | 

Cell 
No.  1 

Cell 
No.  2 

Cell 
No.  3 

Cel) 
No.  6 

Cell 
No.  9 

Cell 
No.  10 

Coll 
No.  11 

Cell 

No.  12 

0 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

30 

0.00 

0.00 

0.00 

0.01 

0.01 

0.03 

0.01 

0.03 

40 

0.02 

0.00 

0.00 

0.01 

0.01 

0.03 

0.03 

0.05 

50 

0.04 

0.03 

0.03 

0.02 

0.01 

0.07 

0.04 

0.06 

60 

0.06 

0.04 

0.05 

0.03 

0.01 

0.08 

0.04 

0.08 

70 

0.06 

0.05 

0.06 

0.03 

0,01 

0.09 

0.05 

0.08 

80 

0.07 

0.06 

0.07 

0.03 

0.02 

0.10 

0.08 

0.10 

90 

0.08 

0.03 

0.08 

0.03 

0.03 

0.10 

0.08 

0.11 

100 

0.10 

0.05 

0.09 

0.03 

0.03 

0.12 

0.10 

0.16 

110 

0.11 

0.07 

0.10 

0.04 

0.05 

0.16 

0.10 

0.20 

Cell  number  12  failed  in  shear  at  120  psi  (0.827  N/mm^), 

Conversion  Factors:  1  psi  B  6.895(10)"^  N/mm2 
1  in.  =  25.4  mm 
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TABLE  IV-16.  PLYWOOD  COVER  DEFLECTION  MEASUREMENTS 
FOR  TEST  BOX  NUMBER  5 (Static  Loading  Test) 


Pressure 

(psi) 

Cover  Deflection  (inches) 

Cell 
No.  1 

Cell 

No.  2 

Cell 
No.  3 

Cell 
No.  6 

Cell 
No.  10 

Cell 
No.  11 

Cell 

No.  12 

0 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

— 

30 

0.01 

0.02 

0.02 

0.00 

I 

0.02 

40 

0.02 

0.03 

0.03 

0.00 

0.03 

50 

0.03 

0.04 

0.08 

0.01 

0.08 

0.05 

0.07 

60 

0.04 

0.05 

0.09 

0.02 

0.09 

0.06 

0.08 

70 

0.05 

0.06 

0.11 

0.03 

0.10 

0.08 

0.09 

80 

0.07 

0.08 

0.11 

0.03 

0.10 

0.09 

0.10 

90 

0.08 

0.09 

0.13 

0.04 

0.12 

0.10 

0.10 

100 

0.10 

0.10 

0.17 

0.05 

0.15 

0.10 

0.11 

110 

0.11 

0.11 

0.20 

0.05 

0.19 

0.12 

0.11 

no* 

0.12 

0.11 

0.22 

0.08 

0.20 

0.14 

0.16 

liO* 

0.13 

0.11 

0.23 

0.07 

0.20 

0.14 

0.18 

no* 

0.13 

0.11 

0.23 

0.08 

0.20 

0.13 

0.17 

no* 

0.13 

0.11 

0.23 

0.08 

0.20 

0.14 

0.17 

90 

0.10 

0.10 

0.20 

0.06 

O.lfi 

0.12 

0.12 

no 

0.15 

0.10 

0.23 

0.08 

0.  .0 

0.16 

0.7  7 

30 

0.05 

0.04 

0.10 

0.03 

C.iO 

0 .08 

0.08 

0 

0.03 

0.02 

0.00 

0.04 

0.02 

0.03 

no 

0.18 

0.11 

0.28 

0.07 

0.20 

0.13 

0.17 

0 

0.03 

0.01 

0.03 

0.01 

O.C7 

0.03 

0.04 

120 

0.18 

0.12 

0.30 

0.07 

0.21 

... 

MM 

0.19 

Cell  number  3  failed  in  shear  at  130  psi  (1.03  N/mm2). 

*Pressure  was  reduced  to  90  psi  (0.621  N/ram2)  and  then  increased  back  to 
110  psi  (0.758  N/mm2). 

Conversion  Factors:  1  psi  =  6.895(10)“2  N/mm2 

1  in.  =  25.4  mm 
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TABLE  IV-17.  PLYWOOD  COVER  DEFLECTION  MEASUREMENTS 
FOR  TEST  BOX  NUMBER  6  (Static  Loading  Test) 


Pressure 

(psi) 

Cover  Deflection  (Inches) 

Cell  No.  1 

Cell  No.  3 

Cell  No.  10 

Cell  No.  12 

0 

0.00 

0.00 

0.00 

0.00 

30 

0.02 

0.02 

0.01 

0.00 

AO 

0.03 

0.03 

0.02 

0.00 

50 

0.03 

0.03 

0.03 

0.01 

60 

0.0A 

0.03 

0.0A 

0.03 

70 

0.05 

0.05 

0.0A 

0.03 

80 

0.06 

0.05 

0.05 

0.0  A 

90 

0.06 

0.06 

0.06 

0.0A 

100 

0.09 

0.08 

0.07 

0.05 

110 

0.10 

0.09 

0.1A 

0.09 

0 

0.05 

0.00 

0.05 

0.01 

Cell  nunber  1  failed  In  (hear  at  160  pal  (1.10  N/an2). 

Conversion  Factora:  1  psi  ■  6.895(10)"3  N/n»2 
1  In.  *  25. A  mm 


TABLE  IV-18.-  PLYWOOD  COVER  DEFLECTION  MEASUREMENTS 
FOR  TEST  BOX  NUMBER  7  (Static  Loading  Test) 


Pressure 

(psi) 

Cover  Deflection  (inches) 

Cell  No.  1 

Cell  No.  3 

Cell  No.  10 

Cell  No.  12 

0 

0.00 

0.00 

0.00 

0.00 

6C 

0.02 

0.02 

0.03 

0.02 

70 

0.03 

0.03 

0.03 

0.02 

80 

0.03 

0.03 

O.OA 

0.03 

90 

O.OA 

O.OA 

0.05 

0.03 

100 

0.05 

0.05 

0.07 

O.OA 

110 

0.06 

0.07 

0.98 

0.05 

120 

0.08 

0.08 

0.09 

0.06 

Cell  number  1  failed  in  shear  at 


2A0  psi  (1.65  N/ma2). 


1  psi  =  6. 895 (10) ~3  N/ma2 
1  in.  »  25. A  mm 
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Conversion  Factors: 


TABLE  IV- 19 4  PLYWOOD  COVEP.  DEFLECTION  MEASUREMENTS 
FOR  TEST  BOX  NUMBER  8  (Dynamic  Loading  Test) 


Pressure 

(psl) 

Cover  Deflection  (inches) 

Cell  No.  6 

0 

0.000 

45.5 

0.033 

68.1 

0.038 

81.8 

0.042 

93.2 

0.047 

107.8 

0.081 

119.3 

0.091 

130.7 

0.095 

TABLE  IV-20.  PLYWOOD  COVER  DEFLECTION 
MEASUREMENTS  FOR' TEST  BOX  NUMBER  10 
(Static  Loading  Teat) 


tell  nuabor  5  failed  In  bending  at  142.1  psl  <0.980  N/m2). 
Conversion  Factors:  1  pal  «  6.895(10)-3  N/tm2 


1  In.  *  25.4  u: 


Pressure 

(psl) 

Cover  Deflection  (Inches) 

Cell  No.  1 

Cell  No,  3 

Cell  No.  10 

0 

0.00 

0.00 

o.oc 

10 

0.00 

0.00 

0.00 

40 

0.00 

0.00 

0.00 

50 

0.02 

0.00 

0.00 

to 

0.02 

0.01 

0.00 

70 

0.01 

0.02 

0.01 

80 

0.04 

0.02 

0.01 

90 

0.05 

0.01 

0.02 

100 

0.05 

0.03 

0.01 

110 

0.06 

0.01 

0.01 

120 

0.06 

0.04 

0.05 

TABLE  IV-21.  PLYWOOD  COVER  DEFLECTION 
MEASUREMENTS  FOR  TEST  BOX  NUMBER  11 
(Static  Loading  Test) 


The  center  longitudinal  stiffener  buckled  between  cells 
4  and  7  and  6  and  9  at  a  prcaaure  between  ISO  end  180 
pal  (1.01  and  1.24  K/am‘). 


Conversion  Factors i 


1  psl  •  6.89S(10>'S  H/m2 
1  In.  »  2S.4  an 


Pressure 

(psl) 

Cover  Deflection  (inches) 

1 

Location 
No.  1 

Location 
No.  2 

Location 

No.  3 

0 

0.00 

0.00 

0.00 

5.0 

0.10 

0.15 

0.10 

7.5 

0.14 

0.23 

0.13 

10.0 

0.22 

0.32 

0.20 

12.5 

0.28 

0.42 

0.23 

15.0 

0.30 

— 

0.28 

17.5 

0.37 

0.33 

The  cover  cracked  (bending  failure)  at  the  20  psl  (0.13s)  N/nm2) 
pressure  level. 


Conversion  Factors:  1  psl  =  6.895(10)~3  N/mm2 
1  In.  =  25.4  sen 
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idlnally, 


TABLE  IV-22.  SUMMARY  TABLE  OF  LOADING  TEST  RESULTS 
FOR  THE  PLYWOOD  INSULATION  BOXES  (Cond'd) 


consistent  with  analytical  calculations  of  dynamic  response  presented  in 
Section  V,  which  show  that  rise  times  must  be  about  1  -  2  ms  to  produce  dynamic 
amplification  of  the  loading. 


By  modifying  the  test  boxes,  loading  pressures  up  to  240  psi  could  be 
obtained  before  box  failure  occurred.  It  was  observed  that  the  box  covers 
could  withstand  approximately  0.2-in.  deflections  at  the  centers  of  the  un¬ 
supported  cell  areas  before  failure  of  the  cover  was  observed.  Permanent 
deformations  of  the  cover  (not  total  failure)  were  witnessed  when  pressure 
levels  reached  80-90  psi;  These  deformations  could  be  described  as 
"scalloping"  of  the  unsupported  .cell  areas  with-  the 'magnitudes  of  these 
deformations  in  the  range  of  0.05-  to  0.10-in.  maximum. 


Material  Properties  Tests 


A  series  of  tests  on  plywood  specimens  taken  from  the  Gaz  Transport 
boxes  were  conducted.  The  purpose  of  these  tests  was  to  provide  material 
property  data  for  the  analysis  of  the  plywood  box  and  to  aid  in  the  inter¬ 
pretation  of  tests  to  measure  box  strengths.  These  tests  were  necessary 
because  data  from  tests  by  Mitsui(?6)  and  by  Gaz  Transport (7 7, 78)  gave 
widely  varying  results  for  plywood  strength  and  stiffness. 


All  plywood  tests  were  conducted  with  simple  beams  in  3-point  bend¬ 
ing.  As  shown  in  Figure  IV-23,  the  specimens  were  approximately  50mm  wide, 
and  support  spacing  was  either  121mm,  88mm,  or  50mm.  Spacings  of  121mm  and 
88mm  match  the  widths  of  the  cells  in  the  9-cell  and  12-cell  boxes,  respec¬ 
tively.  The  spacing  of  50mm  was  used  to  check  the  effect  of  very  short 
spans.  All  specimens  were  cut  transverse  (perpendicular)  to  the  grain. * 
This  is  the  weakest  direction  and  also  the  direction  across  the  narrow  di¬ 
mension  of  the  cells  in  the  plywood  boxes.  Principal  bending  in  the  actual 
boxes  will  occur  in  the  cross-grain  direction. 


Specimens  were  cut  from  the  covers  of  boxes  which  had  been  in  ser¬ 
vice  ("old"  boxes)  and  from  new  boxes  which  had  never  been  exposed  to  cryo¬ 
genic  temperatures.  Tests  on  the  specimens  from  the  old  boxes  were  per¬ 
formed  at  two  loading  rates.  A  slow  rate  of  2.54  mm/min  represents  a  static 
test,  and  the  fast  loading  rate  of  508  mm/min  is  the  limit  of  the  test  ma¬ 
chine.  Tests  with  specimens  from  the  new  boxes  were  conducted  at  both  room 
temperature  (RT)  and  at  cryogenic  temperature  (CT).  All  tests  were  per¬ 
formed  at  the  slow  loading  rate. 

Results  for  the  old  plywood  are  presented  in  Table  IV-23.  Consider¬ 
able  variation  in  the  stiffness  and  stress  values  are  evident.  Stress  values 
range  from  a  low  of  20.9  N/mm2  to  a  high  of  65.5  N/ram2,  a  factor  greater 
than  3.  Average  values  of  0  and  E  for  the  two  different  loading  rates  and 
overall  are: 


E,  N/mm2  (psi) 


Fast  Loading 
Slow  Loading 
Overall 


'40.1  (5,816) 

49.4  (7,165) 

44.2  (6,411) 


2144  (310,987) 
3797  (550,755) 
2941  (426,592) 


*Grain  direction  of  the  exterior  ply. 
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121  mm 
88  mm 
50  mm 


Although  0  and  E  for  the  fast  loading  rate  are  lower  than  for  the  slow  rate, 
this  difference  was  disregarded  for  three  reasons.  First,  the  maximum  and 
minimum  values  for  stress  occurred  at  the  fast  loading  rate,  which  indicates 
that  there  is  no  consistent  trend.  Second,  even  the  fast  loading  rate  used 
in  the  tests*  is  slow  by  several  orders  of  magnitude  compared  to  the  ex¬ 
pected  response  rate  of  the  plywood  covers  in  the  insulation  boxes.  -Third, 
drop- tests  conducted  by  Mitsui, (76)  which  give  faster  loading  rates,  showed 
no  difference  in  results  from  static  tests.  For  these  reasons,  the  average 
values  from  all  tests  were  taken  as  representative  of  plywood  properties 
from  the  old  boxes.  The  spread  in  the  results  for  different  specimens  im¬ 
plies  that  there  is  considerable  variation  in  the  plywood  even  for  a  single 
box  cover.** 

Results  for  the  new  plywood  boxes  are  given  in  Table  IV-24.  Tests 
were  performed  at  both  cryogenic  and  room  temperatures.  Plywood  strength 
at  cryogenic  temperatures  is  consistently  higher  than  at  room  temperature. 
Also,  scatter  is  scmcw^at  less  in  the  strength  results  for  the  new  boxes. 

In  the  room  temperature  tests,  for  example,  the  values  differ  by  a  factor 
of  2.45  as  opposed  to  3.13  for  the  old  boxes,  and  this  is  based  on  about 
twice  as  many  tests  for  the  new  plywood. 

Plywood  data  for  all  tests,  conducted  by  Mitsui,  Gaz  Transport  and 
SwRI,  are  summarized  in  Table  IV-25.  These  data  include  one  test  on  the 
new  box  material  which  is  compressive  parallel  to  the  grain  of  the  exterior 
plies.  This  test  was  to  give  some  indication  of  the  basic  strength  of  the 
vertical  stiffeners  in  the  boxes. 


*Fast  loading  rate  used  in  these  tests  was  the  highest  rate  obtainable  with 
the  test  device. 

**A11  specimens  from  the  old  boxes  were  cut  from  a  single  box  cover. 
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Specimens  are  Identified  by  Box  ^Submerged  in  Liquid  Nitrogen  <*  -32l#F 

Conversion  Factors:  1  psi  “  6.S95(10)~3  N/«*2 


As  noted,  considerable  scatter  was  found  in  tests  on  the  specimens 
from  the  old  box  cover.  This  cover  had  sustained  localized  visible  damage, 
but  specimens  were  cut  from  areas  which  appeared  to  be  undamaged.  However, 
the  large  scatter  may  indicate  that  some  unnoticed  damage  was  present  in 
some  of  the  specimens  tested.  Thus,  plywood  properties  for  in-service 
boxes  were  established  from  tests  on  both  the  old  and  new  boxes.  Tests  on 
the  new  boxes  were  used  to  establish  th3  minimum,  average,  and  maximum 
values  for  both  cryogenic  temperature  (CT)  and  room  temperature  (RT). 

These  values,  taken  from  Table  IV-24,  are  listed  in  Table  IV-26. 

TABLE  IV-26.  RANGE  C/F  PLYWOOD  PROPERTIES 


Birch  Plywood,  New 

Mipimua  1 

Properties  Perpendicular 
to  Exterior  Plywood 

BT 

CT 

o,  N/mm2  (pel) 

31.17 

(4,538) 

59.23 

(8,587) 

B,  N/««2  (pel) 

1213 

(175,946) 

3801 

(522,325) 

77.80  103.70 

(11.28$)  (15,042) 


A  reduction  factor,  RF,  for  in-service  conditions  was  established  by 
the  ratio  of  average  RT  values  for  the  old  and  new  plywood.  From  Tables 
IV-23  and  J.V-24,  these  average  RT  values  are: 


Old  plywood:  O^vg  *  44.20  N/mm2  (4,411  psi) 

eAVG  "  2941  N/mm2  (426,592  psi) 

New  Plywood:  a^y q  ■  55.96  N/mm2  (8,117  psi) 

Eavg  ■  3070  N/mm2  (445,304  ps 

For  "in-service’1  conditions,  the  reduction  factors  become: 


55.96  N/mm2  (8,117  psi) 
3070  N/mm2  (445,304  psi) 


For  a: 


44.20 

55.96 


For  E: 


RF  - 


0.96 


V.  TASK  3  -  ANALYTICAL  STUDIES 


V.l  Response  Prediction  Method 

Analytical  studies  were  performed  to  develop  wall- response  prediction 
methods  for  prismatic  tanks.  The  methods  were  developed  for  and  applied  to 
the  two  specific  tank  designs  (membrane  and  prismatic)  for  which  structural 
details  were  available.  Example  calculations  for  these  two  tank  types  are 
given  in  Appendix  D.  Also  in  this  appendix,  predicted  wall  response  was 
compared  to  the  experimental  response  data  presented  in  Section  IV. 5,  to 
evaluate  the  validity  of  the  analytical  methodology.  Pressure  type  tanks 
were  not  analyzed  since  impulsive  sloshing  pressures  on  the  walls  of  spheri¬ 
cal  and  cylindrical  tanks  are  low  and  are  not  significant  in  the  wall  design. 
Sloshing  forces  are  created  on  cylindrical  and  spherical  tanks,  but  only  the 
resultant  forces  are  of  interest.  Methods  of  predicting  these  forces  are 
presented  in  Section  III. 2,  and  design  for  these  forces  is  discussed  in  Sec¬ 
tion  VI. 

The  response  o'  structures  to  dynamic  loads  can  be  quite  complicated; 
however,  simplifying  assumptions  can  be  made  which  make  the  response  calcu¬ 
lation  tractable  and  suitable  for  design  purposes.  One  approach  is  to  as¬ 
sume  or  prescribe  a  deformation  pattern  for  the  structure.  When  the  defor¬ 
mation  pattern  is  prescribed,  then  the  displacements  at  all  points  in  the 
structure  are  known  when  the  displacement  at  only  one  point  in  the  structure 
has  been  determined.  This  is  also  true  of  velocity  and  acceleration;  thus, 
the  transient  response  of  the  structure  can  be  described  by  one  degree  of 
freedom  (dof) .  To  apply  this  approach,  an  equivalent  one-dof  system  must 
be  developed.  Properties  of  the  equivalent  one-dof  system  are  computed  from 
energy  principles  as  explained  in  Appendix  C.  This  explanation  is  centered 
around  the  work  by  Biggs, (79)  but  the  basic  principles  are  the  same  for  all 
similar  approaches. 

This  method  of  calculating  the  response  of  structures  to  dynamic  load¬ 
ing  is  used  widely  for  design  purposes  (see  References  80  through  83)  be¬ 
cause  of  its  simplicity  and  ease  of  use.  Comparisons  between  calculated 
and  measured  behavior  have  been  made  by  Westine  and  Baker, (84)  Westine  and 
Cox, (85)  and  Westine,  Cox  and  Esparza; (86)  and  the  one  degree  of  freedom 
approximation  has  been  shown  to  give  good  results.  Although  much  work  of 
this  nature  has  been  for  elastic-plastic  behavior,  for  the  design  of  ship 
tanks  only  elastic  behavior  is  required.  Our  approach  in  this  work  is  to 
assume  that  the  structure  will  be  designed  for  elastic  behavior,  even  for 
the  worst  case  loadings.  This  is  certainly  a  necessity  for  structures  made 
of  plywood  which  exhibit  brittle  behavior  at  cryogenic  temperatures.  In 
other  materials  some  plastic  behavior  could  perhaps  be  tolerated,  but  it  is 
not  treated  here. 

As  explained  in  Appendix  C,  methods  for  computing  equivalent  one-dof 
systems  usually  prescribe  the  static  deformed  shape  as  the  deformation  pat¬ 
tern;  however,  it  was  also  shown  that  for  symmetrical  structures  with  uni¬ 
form  loading,  the  fundamental  mode  shape  gives  almost  identical  results. 

Three  other  important  points  are  made  in  Appendix  C. 
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(1)  Solutions  for  elastic  behavior  can  be  based  upon  a  dynamic 
load  factor  (DLF)  which  is  a  function  of  the  load-time  history 
and  of  the  fundamental  period  of  the  structure. 

(2)  The  dynamic  load  factor,  when  multiplied  by  the  peak  force 
which  occurs  in  F(t),  gives' an  equivalent  static  force  for 
design  purposes. 

(3)  For  elastic  behavior  of  structures  which  are  subjected  to 
loads  of  long  duration,  relative  to  the  structure's  funda¬ 
mental  period,  the  equivalent  static  force  gives  good  esti¬ 
mates  of  maximum  shear  reactions  in  the  structures. 

Locally,  sloshing  pressures  are  uniform,  so  that  the  assumption  of  the  fun¬ 
damental  mode  of  vibration  as  the  deformation  pattern  is  valid  for  struc¬ 
tural  components  in  the  tank.  Thus,  the  approach  used  for  the  analysis  of 
ship  tanks  will  be  to  compute  a  dynamic  load  factor  (DLF)  for  structural 
elements,  which  is  based  upon  the  fundamental  period  of  the  structure  and 
the  worst  case  pressure-time  history.  A  static  analysis  of  the  structure 
is  then  performed  using  the  peak  load  multiplied  by  DLF. 

To  apply  this  approach ,one  must  be  able  to  compute  the  fundamental 
frequency  of  the  component  being  analyzed,  establish  the-  maximum  DLF  to  be 
expected  from  the  sloshing  pressures  (using  the  fundamental  period),  deter¬ 
mine  the  peak  value  of  the  sloshing  pressure  or  load,  and  perform  a  static 
analysis  of  the  structural  component  for  the  design  loading,  l.e.,  (Peak 
sloshing  pressure)  x  DLF.  Procedures  for  computing  structural  frequencies 
and  establishing  the  DLF  are  given  in  the  next  section.  Guidelines  are  also 
given  for  the  analysis  of  the  structure. 

V.2  Design  Procedures 


V.2.1 


laraic  Load  Factors 


Measurements  of  rise  and  decay  times  for  sloshing  pressure,  as  re¬ 
ported  in  Section  IV. 4,  showed  no  consistent  relationship  between  the  rise 
time,  decay  time,  and  peak  pressure.  Thus,  it  was  not  possible  to  estab¬ 
lish  correlations  which  might  show  that  maximum  sloshing  pressures  are 
associated  with  the  minimum  rise  times  or  that  the  opposite  is  true.  Like¬ 
wise,  no  relationship  could  be  established  between  peak  pressure  and  decay 
time  or  between  rise  time  and  decay  time.  Thus,  a  complete  range  of  values 
must  be  considered  if  the  worst  case  loading  for  any  particular  structure 
is  to  be  obtained. 

The  wave  form  for  typical  sloshing  pressures  is  shown  in  Figure  IV-5, 
Section  IV. 4,  and  the  ranges  of  values  for  pressure  and  time  are  given  in 
Tables  IV-12  and  IV-13  and  in  Appendix  B.  These  values  are  in  nondimen- 
sional  form.  To  better  understand  the  range  of  pressures  and  times,  maxi¬ 
mum,  minimum,  and  average  values  were  extracted  from  Tables  IV-12  and  IV-13 
and  B-13  through  B-24.  These  values  are  given  in  Table  V-l,  both  in  non- 
dimensional  form  and  for  an  LNG  tank  which  is  36  m  (118  ft)  long.  This 
length  is  typical  of  prismatic  tanks  in  today's  LNG  ships. 
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TABLE  V-l.  RANGE  OF  PARAMETERS  WHICH  DESCRIBE 
SLOSHING  PRESSURES  IN  LNG  TANKS 


r—  ■  . . 

Slosh  Wave 
Parameters 
(See  Fig.  IV-5) 

Nondimensional 

36  m  (118-ft)  Ship  Tank 

Min. 

Max. 

Min. 

1 

Max. 

2 

0.0005 

0.1677 

0.7583 

0.0010 

0.3213 

1.4523 

P** 

r2 

0.0297 

7.0008 

33.7277 

0.0700 

16.4660 

79.3280 

T3 

0.0350 

0.3139 

0.9489 

0.0670 

0.6012 

1.8173 

>3 

0.0000 

0.5989 

2.2346 

0 

1.4090 

5.2560 

*4 

0.0548 

0.8491 

1.4716 

0.1050 

1.6262 

2.8184 

P4 

0.0000 

0.2719 

1.5348 

0 

0.6400 

3.6100 

T5 

0.1642 

1.3582 

2 .0900 

0.3145 

2.6012 

4.0028 

P5 

0.0000 

0.0006 

0.3310 

0 

0.0010 

0.7790 

Nondimensional  times  are  T  /g/j,  ;  dimensional  times  are  in  seconds . 

Nondimensional  pressures  are  P/pgll$.  For  the  36  m  (118-ft)  tank,  pres¬ 
sures  are  in  psi  per  0.1  radian  of  excitation  amplitude  <j».  Multiply  by 
0.006895  to  obtain  pressure  in  MPa  (N/ram^) , 


In  the  discussions  which  follow,  scaled  times  for  the  "typical"  36  m  (118-ft) 
LNG  tank  will  be  used.  For  different  tank  sizes,  the  times  will  scale  by 
the  square  root  of  the  tank  length  ratio,  increasing  for  larger  tanks  and 
decreasing  for  smaller  tanks.  These  discussions  lead  to  a  method  for  deter¬ 
mining  the  dynamic  load  factors  which  is  applicable  to  tanks  of  any  size. 

Note  from  Table  V-l  that  the  rise  times,  T2,  vary  from  a  minimum  of 
0.001  sec  to  a  maximum  of  1.4523  sec  and  that  the  decay  times,  T3,  for  the 
pressure  spike  vary  from  0.601  to  1.817  sec.  Also,  the  "trailing  pressure" 
(P3  and  P4)  is  usually  much  lower  than  the  pressure  spike  and  lasts  for  sev¬ 
eral  seconds.  Typical  structures  in  LNG  tanks  have  fundamental  periods  of 
less  than  50  ms;  thus,  they  will  be  excited  by  the  pressure  spike.  The 
trailing  pressure  will  act  as  a  static  load  and,  having  a  much  lower  magni¬ 
tude  than  the  pressure  spike,  will  be  a  much  less  severe  load.  Thus,  the 
trailing  pressure  is  neglected  in  the  analysis. 

Generally,  the  worst  case  sloshing  load  for  a  structural  component  is 
one  which  has  the  following  properties: 


(1)  highest  pressure 

(2)  shortest  rise  time 

(3)  longest  duration 

For  example,  if  the  rise  time  is  zero  and  the  duration  is  long  relative  to 
the  fundamental  period  of  the  structure,  then  the  loading  will  closely  ap¬ 
proximate  a  step  function  which  gives  a  dynamic  load  factor  of  2.0.  This 
is  seen  in  Figure  C-2  (Appendix  C)  for  the  constant  load  with  finite  rise 
time  (when  the  rise  time  goes  to  zero)  and  for  the  square  and  triangular 
force  pulses  (with  zero  rise  time)  when  their  duration  is  long. 

A  good  approximation  for  the  worst  case  loading  is  then  to  combine 
the  shortest  rise  time  with  the  longest  decay  time.  Figure  V-l  shows  such 
an  approximation  for  the  pressure  spike  at  resonant  sloshing.  The  rise 
time  of  1.0  ns  is  based  on  the  shortest  nondimensional  value  from  labora¬ 
tory  tests,  and  the  time  to  decay  of  1000  ms  is  a  reasonable  worst  case 
value.  Longer  durations  were  recorded  but  most  were  shorter,  and,  for 
structural  components  with  fundamental  periods  of  50  ms  and  less,  the  dif¬ 
ferences  in  structural  responses  for  decay  times  of  250,  500,  or  1000  ms 
are  insignificant.  This  will  be  demonstrated. 


1  ms  1000  ms 

TIME 


FIGURE  V-l.  SLOSHING  PRESSURE  WITH  MINIMUM  RISE  TIME  AND  LONG  DURATION 
FOR  36  m  (118- ft)  LONG  TANK 
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For  the  pressure-time  history  of  Figure  V-l,  the  DLF  computed  for  a 
one-dof  system  is  given  in  Figure  V-2.  This  dynamic  load  factor  wvs  found 
by  computing  the  maximum  transient  response  of  a  one-dof  system  to  a  force 
produced  by  the  P(T)  in  Figure  V-l.  The  maximum  transient  deformation, 
divided  by  the  deformation  produced  by  the  same  force  applied  statically, 
gives  the  DLF.  Stiffness  of  the  spring-mass  system  was  varied  to  vary  the 
periods.  It  covers  fundamental  periods  of  vibration  from  0.5  ms  to  28  ms. 
Note  that  the  DLF  has  a  maximum  value  of  about  2.0  (1.97)  and  drops  to  1.0 
(no  dynamic  effect)  at  1.0  ms  and  0.50  ms.  These  periods  are  whole  multi¬ 
ples  of  the  rise  time,  and  the  structure  responds  as  though  the  load  were 
applied  statically. 

Because  the  rise  times  vary  in  every  sloshing  cycle,  multiple  rise 
times  must  be  considered  to  obtain  the  worst  case  for  each  frequency  (each 
structural  component);  that  is,  a  structure  with  a  fundamental  period  of 
1.0  ms  would  not  be  designed  for  a  DLF  of  1.0  because  a  slightly  longer  or 
shorter  rise  time  would  produce  greater  response.  If  the  DLF  is  computed 
for  several  different  rise  times,  the  envelope  of  Figure  V-3  is  obtained. 
Curves  for  two  different  rise  times  are  shown  for  illustration.  The 
straight  lines  connecting  the  tops  of  the  curves  for  all  different  rise 
times  greater  than  1.0  ms  form  the  top  of  the  envelope.  These  intersect 
the  DLF  curve  for  tr  *  1.0  ms  because  it  was  the  minimum  rise  time  con¬ 
sidered. 

The  effect  of  shorter  decay  times  is  shown  in  Figure  V-4.  All  curves 
were  calculated  for  a  rise  time  of  1  ms.  The  curve  for  a  decay  time  of 
1000  ms  is  repeated  from  Figure  V-2  and  compared  with  decay  times  of  250  ms 
and  500  ms.  Maximum  amplification  occurs  for  the  longest  rise  times  as  ex¬ 
pected,  but  differences  are  small.  Thus,  by  using  the  envelope  of  Figure 
V-3  for  fundamental  periods  of  less  than  or  equal  to  1.2  ms  and  the  curves 
of  Figure  V-2  for  periods  greater  than  1.2  ms,  the  DLF  is  found  for  the 
worst  case  combination  of  rise  and  decay  times.  Of  course,  the  peak  pres¬ 
sure  is  obtained  from  data  presented  in  Section  III. 

The  above  results  are  strictly  correct  only  for  a  tank  which  is  36  m 
(118  ft)  long.  For  longer  tanks,  the  values  of  DLF  will  be  slightly  con¬ 
servative  because  the  rise  time  will  not  be  as  short  and  decay  times  greater 
t.han  1000  ms  will  ha”e  very  little  effect  on  the  DLF  for  structural  periods 
of  50  ms  and  below.  For  smaller  tanks,  the  results  could  be  unconservative 
because  the  rise  times  will  be  shorter. 

The  envelope  of  the  DLF,  corresponding  to  Figures  V-2  and  V-3,  can  be 
reproduced  for  any  tank  size  using  the  data  in  Appendix  C  if  the  approxima¬ 
tion  is  made  that  the  decay  time  of  the  sloshing  pressures  is  infinite. 
Because  the  maximum  decay  times  are  long  relative  to  the  structural  period 
of  interest  in  LNG  tankers,  this  assumption  is  very  reasonable  and  gives  a 
slightly  conservative  result. 

To  establish  the  envelope  of  DLF  for  a  tank  of  a  particular  size. 
Figure  C-3(b)  ,  Appendix  C,  is  used.  This  graph  gives  (DLF)^^  as  a  function 
of  the  rise  time,  tr,  and  the  structural  period,  T.  A  three-step  procedure 
is  required: 
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FIGURE  V-4.  EFFECT  OF  DIFFERENT  LOAD  DECAY  TIMES, 


(1)  Determine  the  minimum  rise  time,  ty ,  using  the  minimum 
value  of  the  nondimensional  time  T2  / g/i-  in  Table  V-l 
and  the  tank  length. 

(2)  Establish  a  range  of  values  of  T  depending  upon  the  range 
of  structural  periods  of  interest,  e.g.,  0.5  ms  <_  T  £  30  ms 
as  in  Figure  V-2. 

(3)  Determine  the  value  of  (DLF)may  for  each  value  of  T  by  enter¬ 
ing  Figure  C-2(b)  for  the  corresponding  values  of  tr/T  and 
taking  the  maximum  value  of  (DLF)nax  for  all  tr/T  greater 
than  or  equal  to  the  value  used  to  enter  the  figure. 

Table  V-2  Illustrates  this  procedure  for  the  36-m  (118-ft)  tank.  Note  that 
closely  spaced  periods,  T,  are  used  in  some  regions  to  better  define  the 
curve.  The  values  of  DLF  in  Table  V-2  arc  very  nearly  the  same  as  those 
given  by  Figures  V-2  and  V-3. 

It  is  not  necessary  to  establish  the  envelope  of  DLF  for  each  tank 
type.  For  each  structural  component  analyzed,  the  appropriate  value  of  DLF 
can  be  read  directly  from  Figure  C-2(b)  for  the  value  of  tr  calculated  for 
the  tank  size  and  T  calculated  for  the  structural  component.  As  for  Step  3 
in  the  determination  of  the  envelope,  the  maximum  value  of  DLF  is  selected 
for  all  tr/T  >_  the  calculated  value.  The  maximum  DLF  is  selected  for  all 
tr/T  >.  (the  value  used  to  enter  the  figure)  to  account  for  the  fact  that 
rise  times  can  be  greater,  but  hot  smaller,  than  the  minimum  value  used  to 
enter  the  figure. 

V.2.2  Frequency  Calculations 

From  the  review  of  LNG  tank  structural  details  in  Section  III, 3,  it 
is  clear  that  several  different  types  of  structural  elements  are  subjected 
to  sloshing  pressures.  These  include: 

o  beaus 
plates 

-  isotropic 

-  ortho tropic 

o  stiffened  panels 

uniaxial  stiff eneing 

-  biaxial  stiffening 

o  compression  members 

In  addition,  the  elements  are  often  submerged  in  liquid  on  one  side,  and 
added  mass  effects  of  the  liquid  will  alter  natural  frequencies.  Also,  non- 
structural  material,  such  as  insulation,  is  often  attached  to  the  wall  of  a 
tank,  and  this  mass  must  be  included  in  frequency  calculations. 


TABLE  Vt-2, 

,  CALCULATION  OF  ENVELOPEOF  DLF 

USING  FIGURE  C-2(b) 

.  * 
tr/T 

T  (ms) 

DLF  from  Figure  C-2(b) 

0.0625 

16.00 

1.990 

0.1250 

8.00 

1.970 

0.2500 

4.00 

1.900 

0.5000 

2.00 

1.640 

0.7500 

1.33 

1.300 

0.8100 

1.23 

1.215 

1.0000 

1.00 

1.215 

1.2500 

0.800 

1.215 

1.4000 

0.710 

1.215 

1.5000 

0.670 

1.210 

1.7500 

0.570 

1.130 

2.0000 

0.500 

1.130 

2.2500 

0.440 

1.130 

2.4500 

0.410 

3.130 

2.6500 

0.380 

1.100 

2.7000 

0.370 

1.090 

3.0000 

0.330 

1.090 

3.2500 

0.310 

1.090 

3.5000 

0.290 

1.090 

* 

Minimum 

<‘r> 

for  36  m 

(118-ft)  tank  is 

A 


T„  -  t_  -  0.0005  “V 


0.000958  5  0.001  sec 


Frequency  formulas  for  uniform  beams  in  bending  and  for  axial  members 
are  common  in  textbooks  on  structural  vibrations  and  will  not.  be  included 
here.  Formulas  for  orthotropic  plates  and  stiffened  panels  are  less  readily 
available,  and  methods  for  computing  frequencies  of  these  elements  are  given 
in  the  following  paragraphs .  The  equat ion  for  or ttiotroplc  plates  reduces  to 
that  for  isotropic  plates  when  isotropic  material  properties  are  used.  In 
addition  to  the  plate-panel  formulas,  a  formula  for  beam-.,  in  which’  shear  is 
significant  is  included.  Plates  with  stiffeners  can  often  be  analyzed  as  a 

beam  which  includes  one  of  the  stiffeners  and  part  of  the  stiffened'  plpte. 

*****  !  , 

The  formulas  presented  for  stiffened  panels  are  for  uniform  spacing 
of  uniform  stiffeners.  For  non-uniform  panels  and  girders,  more-general 
numerical  procedures  are  needed.  Finite -element  methods  are  very  well 
adapted  to  the  analysis  of  complicated  geometries.  Many  other  methods  have 
also  been  developed  specifically  for  stiffened  panels  and  grillages.  See, 
for  example,  the  work  of  Smith, (87)  Madsen, (88)  Chang, (89)  Chang  and  Michel- 
son,(90)  and  Chang  and  Pilkey.(91) 

V.2.2.1  Orthotropic  Plates  (Reference  92) 

Frequency  equation: 


where  p  ■  plate  material  density 

a  ■  long  dimension  of  the  plate  (x-direction) 
b  ■  short  dimension  of  the  plate  (y-direction) 


D 


x 


Exh‘ 


12(1  -  v  v  ) 
x  y 


(V-2a) 


D 


y 


12(1  -  v  v  ) 
x  y 


(V-2b) 


D 

xy 


D  v 
x  y 


(V-2c) 


where 
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h  *  plate  thickness 
v  *  Poisson’s  ratio 

E  ,E  *  elastic  moduli  parallel  to  x  and  y  axes,  respectively 
x  y 

6  *  shear  modulus 

The  constants  A,  B,  and  C  are  given  in  Table  V-3  for  different  boundary 
conditions . 

V.2.2.2  Panels  Stiffened  in  One  Direction  (Reference  93) 

Frequency  equation: 


Panel  geometry  and  definition  of  terms: 


(V-3) 


5  ■  material  elastic  modulus 
p  ■  material  density 
h  3  plate  thickness 

I  3  second  moment  of  area  of  stiffener  and  associated  plate 

S  3  sectional  area  of  stiffener  and  associated  plate  3 

A  +  p  h 
s  ^ 

A  3  sectional  area  of  stiffener 
s 

n  3  number  of  stiffeners 
s 
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TABLE  V-3.  FREQUENCY  COEFFICIENTS  FOR  THE  FUNDAMENTAL 
MODE  (EQUATION  V-l) 


CASE  NO. 


BOUNDARY  CONDITIONS 


^///////. 


rrrrr 


2 


4.730 


4.730 


151.300 


CJ 


4.730 


3.927 


141.379 


4.730 


121.396 


4m .  . . 

4 

4 

4  "rrrr  rr 


i 


3.927 


3.927 


3.927 


132.118 


113.444 


97.40? 


FIXED  BOUNDARY 

_  SIMPLY  SUPPORTED  BOUNDARY 
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The  effect  of  submergence  is  described  in  Section  V.2.2.5.  To  account  for 
uniformly  distributed  added  mass,  such  as  insulation  bonded  to  the  panel, 
adjust  the  density  p  as  follows. 


p'  -  P 


WA  +  WP 
W„ 


(V-4) 


where  W_  »  total  mass  or  weight  of  the  panel  * 

*  i 

W.  ■  total  mass  or  weight,  in  excess  of  the  panel  structure, 
but  which  moves  with  the  panel  during  vibration 

V.2.2.3  Panels  Stiffened  in  Both  Directions  (Reference  93) 

Frequency  equation: 


f 


1.5764 


Hz 


(V-5) 


Geometry  and  definition  of  terms: 
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n  *  number  of  stiffeners  (parallel  to  the  x-direction) 
s 

m  =  number  of  girders  (parallel  to  the  y-direction) 
s 

h  *  plate  thickness 


I  ,1  *  second  moment  of  area  of  stiffeners  and  girders, 

8  8  respectively,  plus  the  associated  plating 

A  ,A  *  cross-sectional  area  of  stiffeners  and  girders, 

8  respectively 


S 

s 


A  +  p  h  + 
s  *x 


A  +  p  h  + 

g  y 


X 

ms 

(V-6a) 

p  »h 

n_ 

(V-6b) 

Corrections  for  added  mass  effects  associated  with  submersion  in  a  fluid 
are  given  in  Section  V.2.2.5.  Account  for  a  uniformly  distributed  mass, 
such  8,a  Insulation  which  is  attached  to  the  pane] ,  by  adjusting  the  panel 
thickness  e  in  the  equation  for  S8  and  Sg. 

V.2.2.4  Beam  Formula  for  Plate  Stiffener  (Reference  87) 

Frequency  formula  for  the  fundamental  mode  (assumes  simple  support  of 
the  stiffener) : 


f 


HZ 


(V— 7) 


where 


V.2.2.5 


El  *  flexural  rigidity  of  the  plate-stiffener  combination 
GA  »  shear  rigidity  of  the  plate-stiffener  combination 
L  *  beam  length 

m  =  mass  per  unit  length  of  the  stiffener  and  associated 
plating 

Effect  of  Submergence  (Reference  93) 


To  obtain  the  frequency,  f^,  of  a  plate  with  one  side  exposed  to  air 
and  the  other  side  exposed  to  a  liquid,  modify  the  frequency  calculated  in 
air,  fa,  by  the  following  formula: 


(V-8) 


'  (V-9) 


(V-10) 

»  density  of  the  liquid 

Pp  *  density  of  the  plate 

a  -  plate  length 

b  ■  plate  breadth 

h  53  plate  thickness 

With  both  sides  of  the  plate  submerged,  becomes 


£z  =  V  * 


where 


K 


K  + 


K  ■  trh 


*  a2  b2 


V.2.3  Response  Calculations 


In  Section  V.2.2,  methods  are  presented  for  calculating  the  fundamen¬ 
tal  frequencies  of  typical  structural  components  in  the  walls  of  LNG  tanks. 
These  components  include  beams,  plates,  and  stiffened  panels.  Based  upon 
the  fundamental  frequency  and  the  minimum  rise  time  of  the  sloshing  forces, 
a  method  is  given  in  Section  V.2.1  for  determining  the  maximum  dynamic  load 
factor  (DLF)  of  the  structural  components.  The  DLF  plus  the  maximum  value 
of  the  sloshing  forces  (sloshing  pressure  multiplied  by  the  appropriate 
area  over  which  it  acts)  gives  an  equivalent  static  load  for  which  the 
structural  component  must  be  designed  to  withstand  the  worst  case  sloshing 
forces.  These  analyses  will  be  performed  for  structural  elements  such  as 
beams,  plates,  and  stiffened  panels. 

No  particular  guidance  is  necessary  for  the  analysis  of  these  compo¬ 
nents;  however,  coupling  between  structural  components  can  and  does  occur 
in  complex  structures  such  as  the  wall  of  an  LNG  tank.  For  these  cases*  the 
analyst  must  decide  whether  or  not  the  various  components  of  the  tank  wall 
can  be  analyzed  individually  or  should  be  analyzed  as  a  multi-degree  of 
freedom  system.  If  each  component  in  the  tank  structure  is  to  be  analyzed 
individually,  then  it  is  usually  conservative  to  assume  rigid  support  for 
each  element  (neglecting  flexibility  in  the  supporting  structure)  and  that 
the  sloshing  forces  are  transferred,  undistorted,  from  one  member  to  the 
next.  Often  the  response  is  attenuated  by  coupling,  but  it  can  be  in¬ 
creased  as  well. 


109 


aBfflBdssssBassMSsaaaoMBBBm 


To  resolve  such  questions  which  arise  when  applying  simple  analytical 
procedures,  a  multi-degree  of  freedom  solution  is  required.  These  can  be 
obtained  by  finite-element  method  or  by  similar  procedures.  If  a  transient 
response  calculation  is  made,  then  the  actual  F(t)  produced  by  the  slosh¬ 
ing  pressure  is  applied  to  the  structure.  To  be  certain  that  the  maximum 
response  is  obtained,  the  calculation  must  be  performed  with  the  minimum 
rise  time  and  with  several  rise  times  of  longer  duration.  Such  an  analysis 
procedure  can  be  quite  costly  to  perform,  if  many  dynamic  dof’s  are  in¬ 
cluded  in  the  calculation. 

An  alternate  approach,  which  avoids  a  dynamic  transient  solution  is  to 
perform  an  approximate  modal  analysis  using  the  dynamic  load  factor  deter¬ 
mined  for  a  one-dof  system  (as  already  described).  This  method  is  possible 
because  in  a  modal  superposition  analysis,  each  normal  mode  responds  to  the 
loading  independently  as  a  single  degree  of  freedom  system.  The  general 
procedure  is  as  follows: 

(1)  Calculate  the  mode  shapes  and  frequencies  of  the  multi¬ 
degree  of  freedom  model  of  the  structure. 

(2)  Calculate  the  modal  displacements  and  stresses  for  the 
peak  sloshing  force  (statically  applied). 

(3)  Find  the  dynamic  load  factors  for  each  mode  (each  modal 
frequency)  for  the  minimum  rise  time.  Figure  C-2(b)  is 
used  and  the  values  of  DLF  are  read  at  the  corresponding 
values  of  tr/T^  for  i  ■  1  through  N  modes.  Note  that  these 
are  not  the  envelope  values  described  in  Section  V.2.1. 

(4)  Multiply  the  modal  displacement  and  stresses  by  their 
respective  DLF’s. 

(5)  To  find  the  absolute  maximum  displacement,  sum  the  modal 
contributions.  To  find  the  most  probable  response,  take 
the  square  root  of  the  sum  of  the  squares  of  the  modal 
contributions. 

This  procedure  should  be  repeated  for  several  longer  rise  times  to  be 
certain  that  the  worst  case  loading  is  used.  Alternately,  a  conservative 
value  can  be  obtained  if  in  Step  (3)  values  of  DLF  for  each  mode  are  read 
from  an  envelope  of  DLF  established  for  the  tank. 

This  approximate  procedure  for  calculating  the  response  of  a  struc¬ 
ture  to  dynamic  loading  without  performing  a  transient  solution  is  commonly 
used  in  the  nuclear  industry(94)  for  earthquake  excitation.  A  more  detailed 
explanation  of  the  method  is  given  by  Biggs. (?9) 
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VI.  TASK  4  -  PRESENTATION  OF  RESULTS  - 
DESIGN  METHODOLOGY 


This  chapter  presents  the  methodology  for  designing  LNG  ship  tanks  to 
withstand  slosh  loads.  The  procedures  are  based  on  the  experimental  and 
analytical  studies  described  in  Chapters  III  -  V. 

VI. 1  Current  IMCO  Requirements  and  Proposed  Changes 

VI. 1.1  General 


Unified  Rules  for  Gas  Tankers:  Cargo  Containment,  as  published  by 
IMCO,  uses  the  following  classification  for  LNG  tanks: 

•  Integral 

.  Membrane  (see  Figure  II-2(a)) 

.  Semi-membrane 

•  Independent 

Type  A 

Type  B 

Type  C  (pressure  tanks)  (see  Figure  II-2(b)) 

Appendix  E  includes  definitions  for  all  of  the  above  mentioned  task  class¬ 
ifications.  With  the  exception  of  Type  C  independent  tanks,  all  tanks  must  be 
analyzed  for  local  sloshing  pressures  as  defined  in  Section  VI. 2. 2. 2.  For  in¬ 
tegral,  membrane  and  semi-membrane  tanks,  the  local  sloshing  forces  must  be  ac¬ 
counted  for  in  the  design  of  the  insulation,  and  the  ship  supporting  struc¬ 
ture  must  be  checked  for  adequate  strength.  Local  sloshing  pressures  are 
produced  on  the  longitudinal  bulkheads,  transverse  bulkheads,  and  on  the 
tank  top.  Areas  of  the  tank  which  must  be  checked  for  sloshing  pressures  will 
depend  upon  permissible  filling  depths.  In  addition  to  the  localized  slosh¬ 
ing  pressures,  the  support  attachments  for  all  independent  tanks  must  be 
checked  for  the  resultant  of  the  sloshing  forces.  Dynamic  effects  of  these 
forces  must  be  accounted  for. 

Where  tanks  or  the  ship's  hull  has  been  designed  in  accordance  with 
the  rules  and  formulas  of  the  classification  societies,  then  the  resulting 
scantlings  must  be  checked  for  the  slosh-induced  forces  and  increased  as 
necessary.  This  is  true  for  integral,  membrane,  semi-membrane,  and  Type  A 
independent  tanks.  In  Type  B  independent  tanks,  which  are  analyzed  by  re¬ 
fined  analytical  methods,  the  slosh-induced  forces  must' be  included  in  these 
analyses,  as  described  in  Section  V. 

Inadequate  information  is  available  to  determine  the  phasing  of  the 
slosh-induced  force  with  other  dynamic  forces.  Thus,  stresses  produced  by 
sloshing  are  assumed  to  be  independent  of  other  dynamic  loads  and  so  are 
combined  with  them,  taking  the  square  root  of  the  sum  of  the  squares.  It 


should  be  noted  that  stresses  produced  by  sloshing  forces  add  directly  to 
stresses  produced  by  vapor  pressure,  stationary  thermal  loads,  and  other 
static  loads  which  affect  the  tank,  but  they  are  not  combined  with  them. 

VI. 1.2  Specific  Requirements 

Specific  requirements  for  the  design  of  LNG  tanks  to  resist  sloshing 
forces  are  offered  which  fit  within  the  framework  of  the  IMCO  rules.  The 
following  revisions  and/or  additions  to  the  rules  are  suggested. 

IMCO  Para.  4.3.5  Sloshing  Loads:  Revise  Paragraph  (a)  to  read: 

The  risk  of  significant  loads  due  to  sloshing  induced  by  any  of  the  ship 
motions  mentioned  in  2.4.6  must  be  considered  for  all  filling  levels  less 
than  100%  full. 

IMCO  Para.  4.4.1  Integral  Tanks:  Add  Paragraph  (b) :  If  significant 
sloshing  induced  forces  are  found  to  occur  in  the  tanks  according  to  4.3.5, 
then  ship  hull  scantlings  in  the  transverse  and  longitudinal  bulkheads  and 
in  the  tank  top  must  be  checked  for  the  slosh-induced  forces  and  the  scant¬ 
lings  increased  as  necessary.  Dynamic  effects  of  the  sloshing  forces  are 
to  be  considered. 

IMCO  Para.  4.4.2  Membrane  Tanks:  Revise  (d)  to  read;  Special  atten¬ 
tion  is  to  be  paid  to  the  possible  collapse  of  the  membrane  and  the  insula¬ 
tion  due  to  an  over-pressure  in  the  interbarrier  space,  to  a  possible  vacuum 
in  the  cargo  tank,  to  the  sloshing  effects,  and  to  hull  vibration  effects. 

Add  Paragraph  (f):  Special  consideration  must  be  given  to  the  effect 
of  sloshing  forces  on  the  ship  hull  as  in  4.4.1. 

IMCO  Para.  4.4.4  Independent  Tanks  Type  A:  Add  Paragraph  (c):  If 
significant  slosh-induced  forces  are  found  to  occur  in  the  tank  according 
to  4.3.5,  then  ship  hull  scantlings  in  the  tank  walls,  roof,  and  inter¬ 
mediate  bulkheads  must  be  checked  for  the  slosh-induced  forces  and  the 
scantlings  increased  as  necessary.  Dynamic  effects  of  the  sloshing  forces 
are  to  be  considered. 

IMCO  Para.  4.6  Supports:  Add  Paragraph  4.6.8:  Supports  for  indepen¬ 
dent  tanks  must  be  analyzed  for  the  resultant  forces  produced  by  sloshing 
motions  in  the  tank.  Stresses  produced  by  the  dynamic  sloshing  forces  should 
be  combined  with  the  wave-induced  loads  of  4.6.3  according  to  4.5.1  for  inde¬ 
pendent  dynamic  stresses. 

The  IMCO  paragraphs  referenced  above  ar*  from  the  IMCO  rules  dated 

1976. 


A  design  methodology  is  presented  in  flow-chart  format  in  Table  VI. 1. 
The  design  procedure  is  shown  to  vary  with  tank  type,  i.e.,  independent 
pressure  tanks  (IMCO  type  C),  independent  gravity  tanks  (IMCO  types  A  Si  B), 
and  membranes  and  semi-membrane  tanks.  For  all  tanks,  the  design  sequence 
proceeds  from  (1)  comparing  ship  periods  to  predicted  resonant  slosh  periods 
to  (2)  establishing  loads  based  on  sloshing  pressure  and  force  coefficients 
to  (3)  designing  structures  by  specific  delineated  methods  which  vary  with 
tank  type . 


Each  important  block  in  the  design  sequence  is  coded  on  Table  VI. 1. 
A  description  of  the  procedures  for  each  coded  block  is  given  below. 


VI. 2.1 


Loads  Prediction  Methodology 

Resonant  Slosh  Periods  -  TR 


Resonant  sloshing  periods* which  are  a  function  of  tank  geometry  and 
fill  depth. are  determined  by  the  appropriate  equation  given  in  Table  III-2.  I 

If  partial  filling  is  allowed,  then  Tr  should  be  determined  for  fill  depths 
between  10%  and  98%.  Slosh  periods  considering  both  the  tank  length  (£)  \ 

and  breadth  (b)  must  be  determined  for  this  range  of  fill  depths.  If  par-  j 

tial  filling  is  prohibited,  resonant  periods  should  be  determined  for  fill  I 

depths  from  90  to  98%.  It  is  important  that  these  fill  depths  be  consid-  j 

ered  even  when  partial  filling  is  prohibited  since  boil-off  results  in  LNG  ‘ 

typically  being  carried  in  this  filling  range. 


Ship  Periods  -  T8  ! 

Using  standard  practices,  ship  periods  in  roil,  pitch,  and  surge  should  j 

be  determined  for  comparison  with  the  resonant  sloshing  periods.  In  the 
event  of  overlap  between  ship  periods  and  resonant  sloshing  periods,  either 
partial  filling  (less  than  about  98%  full)  must  be  prohibited  or  the  tank  f 

must  be  designed  to  withstand  slosh  loads.  j 


Design  Pressure  Coefficients  -  KP 

Several  methods  exist  for  establishing  design  pressure  coefficients.  j 
Based  on  the  results  given  in  this  report,  the  maximum  pressure  coefficient  , 
is  *  *25.  The  average  pressure  coefficient  is  KpAVG  *  25.  A  conser-  1 

vative  approach  would  be  to  assume  and  base  all  subsequent  design  | 

analysis  on  this  worst  case  pressure.  An  alternate  approach  would  be  to  j 
adopt  the  techniques  presented  in  Reference  44.  In  this  reference,  a  proba-  1 
bilistic  procedure,  involving  the  response  of  the  liquid  in  the  tank  at  dif-  j 
ferent  conditions  and  the  probability  of  these  operating  conditions  occur-  jj 
ring,  is  undertaken.  Although  there  is  insufficient  detail  in  Reference  44  j 
to  include  a  complete  example  of  how  peak  sloshing  pressures  are  computed  1 
using  the  probabilistic  procedure,  the  major  steps  can  be  identified.  l 


The  procedure  utilizes  a  statistical  distribution  of  pressures  ob¬ 
tained  from  scale  model  experiments  and  couples  this  statistical  definition 


I 
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TABLE  VI-1  DESIGN  METHODOLOGY  FLOW  CHART 
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PREDICT: 

pTm  MAX.  SLOSHING  PRESS.,  Pmax 
MIN.  RISE  TIME,  T2 
MAX.  DECAY  TIME,  T3 _ 


DETERMINE  PRESSURE 
PD-  PROFILE  ON  THE  TANK 
_ WALL _ 

TYPE  B  -  OPTION  1  4r 


COMPUTE  DYNAMIC  RESPONSE  OF 
TANK  STRUCTURE  TO  TRANSIENT, 
Bl-  DISTRIBUTED  PRESSURE 
P(x  or  y,z,t)  USING 
MULTIPLE-DOF  MODEL 


r®  (see  Page  116) 


ESTABLISH  DYNAMIC  LOAD 
FACTORS,  DLF,  FOR 
SLOSHING  PRESSURES 


COMPUTE  EQUIVALENT  STATIC 
DESIGN  LOADS  FOR  THE 
TANK  WALL 


TYPE  B  -  OPTION  2 


TYPE  A 


MEMBRANE  &  SEMI-MEMBRANE 
TANKS 


PERFORM  MODAL  ANALYSIS  AND 
COMPUTE  MODAL  STRESSES  AND 
DISPLACEMENTS  FOR  EQUIVALENT 
STATIC  LOADS 


RESIZE  TANK  SCANTLINGS  AS 
NECESSARY  BASED  UPON 
CLASSIFICATION  SOCIETY 
FORMULAS  AND  EQUIVALENT 
STATIC  DESIGN  LOADS 


CHECK  STRENGTH  OF  INTERNAL 
INSULATION  BY  TESTING 
OR  BY  ANALYSIS 


COMBINE  MODAl  STRESSES  AS 
INDEPENDENT  RESPONSES  TO 
OBTAIN  MOST  PROBABLE 
MAXIMUM  STRESS 


RESIZE  SHIP  SCANTLINGS  AS 
NECESSARY  BASED  UPON 
CLASSIFICATION  SOCIETY 
FORMULAS  AND  EQUIVALENT 
STATIC  DESIGN  LOADS 


© 

(see  Page  116) 


TABLE  VI-1  DESIGN  METHODOLOGY  FLOW  CHART  (CONTD.) 
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of  pressure  distributions  with  information  on  ship  motion  spectral  density 
in  the  frequency  ranges  around  the  resonant  slosh  frequencies.  Using  statis¬ 
tical  techniques,  this  procedure  then  determines  the  worst  case  sloshing 
loads  that  are  likely  to  occur  in  a  20-year  operational  period.  Pressure 
coefficients  greater  or  less  than  Kp  *  125  may  be  predicted  statistically 
depending  on  the  actual  conditions  encountered.  Based  on  the  results  re¬ 
ported  herein,  it  is  recommended  that  tank  designers  use  a  Kp  equal  to  125 
or  the  statistical  approach  of  Reference  44.  A  Kpj^  *  125  is  believed  to 

be  a  highly  conservative  number,  and  its  use  would  preclude  the  use  of  a  more 
complicated  procedure  which  is  based  on  probabilistic  conditions. 

Ship  Motion  Amplitude  -  A1 

To  determine  pressure  magnitudes,  the  amplitudes  in  oach  of  the  impor¬ 
tant  degrees  of  freedom  must  be  established.  Maximum  amplitudes  will  result 
in  setting  maximum  anticipated  pressures  for  design  purposes.  Maximum  amp¬ 
litudes  should  be  considered  for  roll,  pitch,  and  surge. 

Sloshing  Pressures  -  PT 

Sloshing  pressures  are  determined  from  the  pressure  coefficient  and 
amplitude  information  established  above.  Pressures  should  be  predicted  for 
each  of  the  appropriate  tank  walls,  whereby  the  appropriate  length  dimen¬ 
sion,  l  or  b,  is  utilized.  The  following  values  of  the  pressure  coeffic¬ 
ients  should  be  used  depending  on  the  amplitude. 

K£  «  Kp  x/fc  oi  *  <  0.1  (VI-1) 

KjJ  ■  Kp  ^0.075  +  °~9-25)  x/i  or  <J>  >  0.1  (VI-2) 

At  low  amplitude  sloshing,  the  pressures  are  shown  to  be  linear  with 
amplitude,  while  at  amplitudes  greater  than  0.1  the  pressures  do  not  in¬ 
crease  linearly.  (See  Figure  III-5,)  Equation  VI-2  uses  currently  avail¬ 
able  data  which  give  the  effect  of  amplitude  on  pressure  coefficients  and 
provide  a  more  realistic  definition  of  Kp  at  high  amplitudes. 

To  determine  a  characteristic  range  of  rise  and  decay  times  of  the 
dynamic  sloshing  pressures,  the  following  coefficients  should  be  utilized. 
These  values  are  taken  from  Table  V-l. 

Minimum  pressure  rise  time:  K^  *  0.0005 

Maximum  pressure  spike  duration:  =  0.9489 

Thus,  both  maximum  and  minimum  times  are  established  for  subsequent  struc¬ 
tural  response  calculations.  Both  length  and  breadth  dimensions  of  the  tank 
must  be  utilized  for  the  characteristic  length  dimension  in  these  equations. 


Tank  Wall  Pressure  Distribution  -  PD 


The  pressure  distributions  over  the  tank  walls  are  important  in  deter¬ 
mining  the  walls’  responses  to  slosh-induced  loads.  The  following  idealized 
pressure  distributions  are  for  (1)  vertical  tank  walls  and  (2)  the  tank  top. 
The  following  pressure  distributions  have  been  compiled  using  the  currently 
available  data.  Because  the  distributions  are  affected  by  a  great  many 
variables,  these  idealized  distributions  should  be  considered  as  conserva¬ 
tive  estimates. 

For  the  vertical  tank  walls,  the  pressure  may  be  considered  constant 
in  the  horizontal  direction.  In  the  vertical  direction,  the  pressure  vari¬ 
ation  is  described  by  the  following  cosine  function.  The  maximum  pressure 
occurs  at  the  static  liquid  filling  level  and  decays  to  zero  at  levels  20% 
above  and  below  the  static  liquid  level. 

KPZ  ■  I  KfMAX  [X  +  “*  (*4^)]  (VI-3) 

where  Kp^  *  pressure  coefficient  at  level  z 
KPmax  "  maximun  pressure  coefficient 
z  *  distance  from  the  tank  bottom 
h  »  liquid  filling  height 
H  •  tank  height 
and  with  limits: 

h  -  ^r£z_<h+  — 

For  localized  areas  of  4  m^  (2  m  x  2  m)  or  less,  the  pressure  can  be  as¬ 
sumed  to  act  uniformly  over  the  area. 

For  the  tank  top,  the  pressure  may  be  assumed  to  be  constant  (the 
maximum  value)  on  the  corner  areas.  These  areas  should  extend  in  from  the 
port  and  starboard  walls  one-third  of  the  width  of  the  tank  and  from  the 
fore  and  aft  walls  one-third  of  the  length  of  the  tank.  The  remaining  top 
areas  should  be  designed  to  withstand  sloshing  pressure  based  on  Kp^,. 


Total  Sloshing  Forces  -  FR 


The  total  forces  exerted  on  the  tank  bulkheads  by  the  sloshing  liquid 
can  be  determined  from  Figure  II1-6  for  prismatic  tanks  and  from  Figures 
III-10  -  III-17  for  spherical  tanks.  It  is  noted  that  for  prismatic  tanks 
these  forces  represent  the  inertia  forces  that  occur  on  the  tank  wall  due 


to  resonant  liquid  sloshing.  They  do  not  include  the  forces  exerted  on  the 
structure  by  the  inertia  of  tank  structural  components.  For  spherical 
tanks,  both  the  resultant  forces  (I  gures  III-10  and  III-ll)  and  the  slosh¬ 
ing  components  (Figures  III-12  through  III-17)  are  given.  The  resultant 
forces  represent  the  combination  of  static  weight  and  inertial  sloshing 
loads. 

VI . 2 . 2  Tank  Design  Methodology 

Tank  Support  Stress  -  C2 

The  sloshing  reaction  force  is  the  integral  of  the  sloshing  pressures 
over  the  tank  wall  area.  It  is  produced  by  the  motion  of  the  liquid  in  the 
tanks  and  does  not  include  the  inertia  of  the  tank  structure  or  the  weight 
of  the  liquid.  As  noted  in  block  FR,  some  results  are  given  for  spherical 
tanks  which  include  the  weight  of  the  liquid.  Rise  times  and  durations  are 
not  given,  but  it  is  safe  to  assume  that  the  reaction  is  quasi-static  for 
the  response  of  the  tank  on  its  supports,  and  so  the  supports  should  be  ana¬ 
lyzed  for  the  peak  sloshing  reaction  force,  statically  applied.  This  reac¬ 
tion  produces  stress  both  in  the  tank  and  in  the  supporting  structure.  The 
structural  arrangements  for  tank  support  can  be  very  complicated, and  com¬ 
plex  multi-dof  analysis  procedures  are  often  required  to  accurately  compute 
the  resulting  stresses. 

Dynamic  Tank  Structure  Response  -  B1 


Independent  tanks,  Type  B,  are  analyzed  using  refined  analytical  tech- 
nqiues  such  as  the  finite-element  method.  One  approach  (Option  1)  for  cal¬ 
culating  the  stresses  in  the  tank  produced  by  sloshing  pressures  is  to  per¬ 
form  a  dynamic  transient  solution  for  the  actual  pressure-time  or  force-time 
history.  This  calculation  procedure  is  discussed  in  Section  V.2.3.  An 
analysis  of  this  type  should  provide  the  most  accurate  solution  to  the  be¬ 
havior  of  the  tank  for  sloshing  pressures,  but  it  can  be  very  expensive  in 
terms  of  computation  time.  Multiple  analyses,  with  different  loading  con¬ 
ditions  (different  rise  times  for  the  sloshing  pressure)  are  usually  re¬ 
quired  to  determine  the  worst  case  loading,  making  this  approach  even  more 
costly.  Thus,  a  more  practical  approach  is  the  approximate  modal  analysis 
suggested  as  Option  2.  (see  Section  V.2.3). 

Dynamic  Load  Factors  -  LI 


Once  the  minimvn  rise  time  and  maximum  decay  time  of  the  sloshing 
pressures  have  been  determined, the  dynamic  load  factors  for  a  range  of  fre¬ 
quencies  of  the  tank  structure  can  be  established.  The  dynamic  load  fac¬ 
tors  depend  only  upon  the  rise  and  decay  times  of  the  loading  and  the  struc¬ 
tural  vibration  periods.  They  are  independent  of  the  pressure  magnitude. 
This  procedure  is  described  in  detail  in  Section  V.2.1. 

Equivalent  Static  Design  Loads  -  L2 

From  the  dynamic  load  factors  (DLF’s)  and  the  peak  value  of  the  slosh¬ 
ing  pressures,  equivalent  static  design  loads  can  be  determined  for  the  tank 
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structure.  Design  static  pressures  are  simply  the  peak  pressure  multiplied 
by  DLF.  Design  forces  are  the  design  static  pressures  integrated  over  the 
appropriate  wall  area  of  the  tank. 

Modal  Stress  Computation  -  B2 

An  alternate  approach  to  a  dynamic  transient  solution  for  the  computa¬ 
tion  of  stresses  in  the  Type  B  independent  tank  is  to  perform  an  approximate 
modal  solution.  This  approach,  Option  2,  is  based  upon  a  multi-dof  solution 
for  modal  frequencies  and  modal  stresses  produced  by  the  peak  value  of  the 
sloshing  pressure,  applied  statically  (see  Section  V.2.3).  Frequencies, 
mode  shapes,  and  static  stress  are  less  costly  to  obtain  than  a  dynamic 
transient  solution  described  as  Option  1  and,  once  obtained,  the  solutions 
for  stresses  under  different  loading  conditions  (different  rise  times)  are 
easy  to  obtain. 

Combined  Modal  Stresses  -  B3 


Once  modal  frequencies  have  been  calculated,  the  dynamic  load  factors, 
DLF's,  for  each  mode  can  be  determined.  Static  modal  stresses  multiplied 
by  the  modal  DLF  gives  the  modal  stress  for  the  dynamic  loading.  The  most 
probable  maximum  stress  in  the  structure  is  then  obtained  by  combining  modal 
stresses  by  the  "square  root"  method.  This  method  of  combining  stresses 
assumes  that  modal  stresses  are  independent.  An  upper  limit  on  stress  is 
obtained  by  simply  adding  mod^x  values. 

Set  Tank  Scantlings  (Type  A)  -  A1 

Scantlings  for  Type  A  independent  tanks  are  sized  according  to  formu¬ 
las  of  the  classification  society.  Once  the  equivalent  static  sloshing 
pressures  have  been  determined  for  the  tank,  these  same  formulas  can  be  ' 
used  to  size  or  resize  the  scantling  for  the  sloshing  forces.  Alternate 
methods,  such  as  those  suggested  for  Type  B  independent  tanks  or  for  the 
membrane  insulation,  are  also  applicable  to  Type  A  tanks.  Approval  to  use 
alternate  methods  must  be  granted  by  the  classification  society. 

Internal  Insulation  Strength  -  Ml 

Strength  of  the  insulation  which  supports  the  membrane  must  be  adequate 
for  the  sloshing  pressures.  As  demonstrated  in  Appendix  D,  strength  of  in¬ 
sulation  system  can  be  well  approximated  by  fairly  simple  classical  methods 
and/or  laboratory  tests  at  the  equivalent  static  design  pressure.  A  combi¬ 
nation  of  testing  and  analysis  is  recommended  because  of  uncertainties  in 
the  material  properties  of  some  insulation  systems.  Of  course,  more  refined 
analytical  methods  can  be  used,  but  these  are  recommended  only  if  material 
properties  of  construction  materials  are  well  defined. 

Set  Tank  Scantlings  (Membranes)  -  M2 

As  for  Type  A  independent  tanks,  the  scantlings  of  the  ship's  inner 
hull  are  determined  using  formulas  of  the  ship  classification  society.  Once 


the  equivalent  static  sloshing  pressures  have  been  determined  for  the  tanks, 
these  same  formulas  can  be  used  to  size  or  resize  scantlings  for  the  slosh¬ 
ing  forces.  Alternate  methods,  such  as  those  suggested  for  Type  B  indepen¬ 
dent  tanks  or  for  the  membrane  insulation  are  also  applicable  for  Type  A 
tanks.  Approval  must  be  granted  by  the  society  to  use  these  alternate  meth¬ 
ods. 

Combined  Stresses  from  Sloshing  and  Other  Dynamic  Loads  -  SI 


Insufficient  information  exists  to  determine  the  phasing  of  the  slosh 
loads  with  other  dynamic  loads  which  act  on  the  tank;  thus,  the  resulting 
stresses  are  treated  as  independent  with  respect  to  stresses  produced  by 
other  dynamic  loads.  Combination  of  these  independent  stresses  to  produce 
the  most  probable  maximum  dynamic  stress  is  by  the  square- root  method,  which 
is  the  square  root  of  the  sum  of  the  square  of  the  independent  values.  When 
combining  stress  produced  by  localized  sloshing  pressure  with  other  dynamic 
stresses,  stress  produced  by  the  dynamic  liquid  head  is  omitted.  Measured 
sloshing  pressures  include  the  dynamic  liquid  head.  Also,  because  sloshing 
pressures  occur  very  near  the  static  liquid  level,  or  above  as  in  the  case 
of  sloshing  pressure  on  the  tank  top,  the  liquid  head  is  very  small  where 
sloshing  pressures  are  highest. 

Combined  Static  and  Dynamic  Stress  -  S2 

The  resulting  dynamic  stress  is  added  directly  to  stresses  produced 
by  static  loads  such  as  thermal  loads,  still-water  loads,  vapor  pressure 
loads,  etc.,  to  find  the  maximum  stress  for  the  condition  being  investigated. 

VI. 3  Example  Problem  Utilizing  Tank  Design  Methodology 

The  following  example  problem  illustrates  the  use  of  the  tank  design 
methodology  presented  in  Table  VI-1.  The  example  analyzes  a  prismatic  mem¬ 
brane  tank  with  the  Gaz  Transport  membrane  system  (plywood  insulation  boxes) . 
Tank  dimensions  are  arbitrarily  set  at 

l  ■  36  m  (118  ft) 

b  -  36  m  (118  ft) 

H  *  24.4  m  (80  ft) 

and  tank  motion  amplitudes  (Al)  are  assumed  to  be: 

±5°  (  0.087  rad.) 

±10°  (0.175  rad.) 

±0.1 


^PitchMAX 

^RoHmax 

x 

S‘  SurgenAx 
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Tank  filling  levels  during  operation  are  placed  at  90  to  98%,  Current 
regulations  limit  most  LNG  tanks  to  maximum  operational  filling  levels  of 
97-98%,  with  partial  filling  prohibited.  From  Table  VI-1,  the  initial  de¬ 
sign  steps  are  to  calculate  (TR)  the  resonant  periods  of  the  liquid  in  the 
tank,  to  establish  (TS)  ship  motion  periods,  and  to  determine  if  the  tank 
resonant  periods  coincide  with  ship  motion  periods.  For  purposes  of  this 
example,  the  ship  motion  periods  are  assumed  to  be  in  the  range  of  6  to  12 
seconds. 

The  resonant  periods  of  the  liquid  in  the  tank  are  calculated  using 
Equation  III-6  for  prismatic  tanks.  For  filling  levels  between  90  and  98%, 
the  resonant  pitch  and  roll  sloshing  periods  are: 

6.89  sec  _<  T^  <_  6.94  sec 

0.9  £  h/H  £  0.98 

Thus,  the  tank  excitation  periods  overlap  the  resonant  periods  of  the  tank 
liquid  so  there  is  the  potential  for  resonant  sloshing. 

At  this  point,  the  designer  can  decide  to  modify  the  tank  size  or 
geometry  to  change  the  resonant  periods  of  the  liquid  in  the  tank.  In  this 
example,  the  tank  will  remain  unchanged,  so  it  must  be  designed  to  withstand 
resonant  sloshing  conditions. 

For  design  purposes,  the  maximum  nondiraensional  pressure  coefficient 
(KP)  is  125.  Equations  VI -1  and  VI -2  are  then  used  to  correct  this  coeffi¬ 
cient  for  tank  excitation  amplitude  effects.  As  previously  mentioned,  the 
amplitudes  for  this  example  are  assumed  values.  The  corrected  coefficients 
for  pitching,  rolling,  and  surging  motions  (PT)  are: 

for  pitching  (<}»  »  ±5°  or  0.087  rad.):  Kp^  *  125.0 

for  rolling  ($  ■  ±10°  or  0.175  rad,):  Kpj^  *  75.4 

for  surging  (  j  -  ±0.1):  Kp'^  »  125.0 

At  this  point,  it  is  necessary  to  choose  the  pressure  rise  and  dura¬ 
tion  time  values  (K?2  and  K13)  associated  with  the  pressure-time  history  of 
the  dynamic  impact  pressure.  Using  Table  V-l,  the  minimum  nondimen sional 
rise  time  (l^)  Is  0.0005.  The  maximum  spike  duration  (KT3)  is  0.9489. 

Now  that  the  nondimensional  design  pressures  and  times  have  been  es¬ 
tablished,  full  scale  pressures  and  times  can  be  calculated.  Pressures  can 
be  determined  using  Equation  1 1 1. 4.  For  this  example,  the  liquid  density  is 
493.19  kg/m3  (30.79  lbra/ft^).  The  design  pressures  are  found  to  be: 

for  pitching:  Pj^  =  1.892  N/(mm)^  (274.4  psi) 

for  rolling:  Pj^  =  2.295  N/(mm)^  (332.9  psi) 
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for  surging:  =  2.175  N/(mm)2  (315.4  psi) 

Full-scale  times  are  determined  using  Equation  II-2.  The  minimum  pressure 
rise  time  is  1.0  millisecond,  and  the  maximum  pressure  spike  duration  is 
1.82  seconds. 

To  compute  design  loads  for  structural  components  of  the  tank  (block 
LI  of  Table  VI-1),  the  dynamic  load  factors,  corresponding  to  the  rise  and 
decay  times  of  the  sloshing  pressures,  are  determined  for  a  range  of  struc¬ 
tural  frequencies.  This  was  done  for  a  tank  which  is  36  m  (118  ft)  long 
(or  wide)  in  Section  V.2,  and  the  results,  expressed  as  the  envelope  of  the 
dynamic  load  factors  (DLF) ,  are  plotted  in  Figures  V-2  and  V-3. 

To  determine  the  equivalent  static  design  load  for  structural  compo¬ 
nents  of  the  tank  (block  L2) ,  the  fundamental  vibration  period  of  the  com¬ 
ponent  must  be  calculated.  With  the  fundamental  period  known,  the  DLF  is 
read  directly  from  Figures  V-2  or  V-3.  These  steps  were  performed  for  the 
Gaz  Transport  insulation  boxes  and  are  described  in  detail  in  Appendix  D, 
part  D.l.  Two  components  of  the  box  were  analyzed  in  the  example:  the  box 
cover  and  the  longitudinal  stiffener  of  a  primary  box.  Results  for  these 
steps  are: 


Fundamental  Period, 

T 

DLF 

Box  Cover  (12-cell  box) : 

0.93  ms 

1.22 

Longitudinal  Stiffeners: 

1.96  ms 

1.63 

Equivalent  static  design  pressures  for  the  cover  and  stiffener  are  obtained 
by  multiplying  the  peak  sloshing  pressures  from  block  PT  by  the  DLF.  Using 
the  maximum  pressure,  P^ax,  which  occurs  in  roll  for  this  example,  the  de¬ 
sign  static  pressures  are: 

Cover:  P,  =  1.22  (2.295)  =  2.80  N/inm2  (406  psi) 

des  9 

Stiffener:  P^eg  =  1.63  (2.295)  =  3.74  N/nrai  (543  psi) 

Strength  of  the  insulation  boxes  (block  Ml)  was  established  in  Appen¬ 
dix  D  and  expressed  in  terms  of  the  peak  allowable  sloshing  pressure,  Pa. 
This  allowable  pressure  is  the  static  pressure  which  the  box  component  will 
withstand  divided  by  DLF.  The  allowable  pressure,  Pa,  would  be  compared 
with  the  maximum  sloshing  pressure,  Pmax»  rather  than  the  design  pressure, 
pdes*  To  obtain  a  pressure  for  comparison  with  Pjgg,  Pa  ^rom  Appendix  D  is 
simply  multiplied  by  DLF.  From  Appendix  D  we  have 

Cover:  P  =  0.853  N/mra2  (124  psi) 

a 

P®es  =1.22  (0.853)  =  1.04  N/mm2  (151  psi) 
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Stiffener:  P  «  0.778  -  0.937  N /ram  (113  -  136  psi) 
a 

=  1.63  (0.770  -  0.937)  =  1.268  -  1.527  N/mm2 
des  (184  -  222  psi) 

s 

where  P,jes  denotes  the  design  static  strength  of  the  component.  Comparing 
these  values  to  P^gg  computed  in  block  L2,  we  find 

Cover:  P®jeg  -  1.04  N/mra2  (151  psi)  <  P^eg  ■  2.80  N/ram2  (407  psi) 

Stiffener:  P®jes  *  1.268  -  1.527  N/ram2  (184  -  222  psi)  <  P^eg  ■ 

3.74  N/mm2  (543  psi) » 

Thus,  for  this  example,  the  design  strength  of  the  insulation  is  inadequate. 

In  this  example,  only  the  primary  insulation  box  was  analyzed  for  the  slosh¬ 
ing  pressures.  Of  course,  in  practice  the  designer  would  need  to  check  the 
strength  of  the  secondary  box  as  well. 

The  final  step  (block  M2)  in  evaluating  membrane  tanks  for  sloshing 
pressures  is  to  check  the  strength  of  the  ship's  inner  hull  which  supports  the 
membrane  system.  Scantlings  of  the  hull  are  checked  using  formulas  of  the 
classification  societies  and  the  equivalent  static  design  load.  In  this 
example  we  did  not  apply  classification  society  formulas  to  check  hull 
strength;  however,  alternate  methods  of  analysis  are  applicable,  but  appro¬ 
val  of  the  classification  society  must  be  obtained  before  alternate  proced¬ 
ures  can  be  applied.  As  an  example  of  alternate  methods,  see  the  analysis 
of  the  CONCH  independent  tank  in  part  D.3  of  Appendix  D. 

VI. 4  Summary 


The  design  methodology  presented  in  this  chapter  is  based  on  a  thor¬ 
ough  review  of  the  literature  supplemented  by  additional  experiments  and 
analytical  studies.  Since  information  required  for  many  of  the  steps  in 
this  procedure  is  limited,  the  procedures  developed  herein  have  been  based 
on  utilizing  conservative  approaches.  For  example,  pressure  coefficients 
are  based  on  the  worst  case  pressure  coefficients  determined  from  a  wide 
range  of  data  available  from  the  literature.  Likewise,  analytical  methods 
which  are  based  on. simplified  one  degree  of  freedom  analysis  have  been  de¬ 
veloped  with  conservatism  imposed.  Thus,  utilizing  these  techniques,  a 
conservative  design  could  occur.  It  is  believed  that  the  major  steps  in 
the  proposed  design  methodology  are  appropriate,  and  it  is  left  to  the  in¬ 
dividual  designer  to  substitute  more  realistic  loads  when  he  believes  it  is 
appropriate.  As  an  example,  sloshing  loads  could  be  more  realistically  es¬ 
tablished  from  scale  model  tests  on  replica  models  of  interest.  Scaled 
ship  motions  could  be  imposed  on  the  models  at  fill  depths  of  interest  and 
the  resulting  loads  and  load  distributions  determined.  This  would  provide 
more  accurate  data  relative  to  a  specific  design.  However,  if  resonant 
sloshing  for  all  filling  levels  were  investigated,  it  is  likely  that  a 

Kp„.v,  nearing  the  125  value  established  herein,  would  be  approached. 

4  MAX 
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VII  CONCLUSIONS  AND  RECOMMENDATIONS 


As  a  result  of  the  study  reported  herein  the  following  conclusions 
and  recommendations  are  presented. 

VII. 1  Summary  and  Conclusions 

o  A  significant  amount  of  scale  model  sloshing  data  for  both 
prismatic  and  spherical  tanks  and  covering  a  wide  range  of 
fill  depths,  excitation  frequencies  and  amplitudes,  tank  wail 
pressure  measurement  locations  and  total  force  measurements  is 
available  in  the  literature. 

o  A  compilation  of  this  data  and  presentation  in  a  common  format 
has  been  accomplished  to  provide  design  pressure  and  force 
coefficients. 

o  Because  the  complexity  of  resonant  large  amplitude  sloshing 
causes  a  wide  variation  in  the  impact  pressure  magnitude  on 
each  successive  sloshing  cycle,  both  worst  case  and  average 
pressure  coefficients  have  been  defined  for  design  purposes. 

These  pressures  are  assumed  to  potentially  occur  at  any  loca¬ 
tion  on  the  tank  wall  if  partial  filling  is  allowed. 

o  Pressure-time  histories  have  been  established  for  sloshing  dy¬ 
namic  impact  pressures.  No  correlation  was  found  to  exist 
between  impact  pressure  magnitude  and  duration.  A  representa¬ 
tive  range  of  impact  pressure  rise  times  and  durations  were 
established . 

o  Dynamic  loads  test,  where  predicted  full-scale  pressure/time 
histories  were  reproduced  on  representative  segments  of  a  LNG 
membrane  tank  structure,  established  structural  response  and 
failure  characteristics.  It  is  noted  that  failure  occurred 
at  simulated  loads  that  could  occur  during  "worst  case" 
resonant  sloshing. 

o  Combined -degree -of -freedom  resonant  sloshing  test  (pitch  and 
heave,  surge  and  heave)  produced  impact  pressures  no  greater 
than  the  corresponding  pressures  measured  in  single-degree-of- 
freedom  resonant  sloshing  (at  the  same  amplitudes). 

o  Materials  properties  measured  on  plywood  used  in  some  LNG 

membrane  structures  showed  considerable  variation.  This  data 
is  valuable  in  the  design  of  tank  wall  structures  using  such 
materials. 

o  Analytical  methods  were  developed  to  predict  structural  response 
to  dynamic  slosh  loads.  Localized  areas  of  the  tank  structure 
are  treated  as  equivalent  one  degree- of- freedom  systems  and  a 
dynamic  load  factor  (DLF)  is  established  based  on  the  fundamental 
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period  of  the  structure  and  a  worst  case  pressure-time  history. 
Examples  for  typical  structural  elements  are  given  and  show 
that  some  typical  structures  are  not  designed  to  withstand 
worst  case  impact  loads. 

o  Design  methods  are  presented  for  independent  pressure  and 

gravity  tanks  and  membrane  and  semi-membrane  tanks.  The  design 
sequence  proceeds  from  1)  comparing  ship  periods  to  predicted 
resonant  slosh  periods  to  2)  establishing  loads  based  on  loads 
coefficients  to  3)  designing  structures  by  specific  delineated 
methods  which  vary  with  tank  type. 

VII. 2  Recommendations 


o  Additional  literature  review  and  experiments  are  required  to 
evaluate  loads  on  liquid  bulk  carriers  with  significant  tank 
internal  structures  (i.e.,  typical  crude  and  chemical  carriers). 
Emphasis  should  be  placed  on  localized  loads  on  web  frames, 
pump  columns  and  swash  bulkheads.  Total  forces  should  also 
be  established. 

o  Full  scale  LNG  sloshing  data  currently  being  obtained  should 

be  analyzed  in  light  of  the  results  reported  herein.  Simulated 
full  scale  ship  motion  conditions  on  scale  models  should  be 
undertaken  with  model  sloshing  loads  compared  to  full  scale. 

o  Additional  model  experiments  should  be  conducted  to  further 

establish  effects  of  excitation  amplitude  on  measured  pressure 
coefficients.  Also  end  wall  pressure  distributions  for  a 
complete  range  of  prismatic  tank  fill  levels  should  be  estab¬ 
lished. 

o  Analytical  methods  and  design  procedures  presented  herein 

should  be  further  refined  after  the  above  recommended  work  is 
accomplished  and  feedback  from  industry  is  obtained. 
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APPENDIX  A 


SLOSHING  FACILITIES  FOR  ANGULAR  MOTION  AND  SIMULTANEOUS  HORIZONTAL  AND  VERTICAL  MOTION 
A.l  Angular  Motion  Slosh  Rig 

The  angular  motion  slosh  test  facility  (see  Figure  A.l)  consists  of  a 
trunnion-mounted  "picture  frame"  weldment  made  from  8-inch  square  steel 
tubing.  The  trunnions  are  pivoted  in  pillow  blocks  mounted  on  top  of  struc¬ 
tural  steel  pedestals  attached  to  a  base  frame,  which  is  anchored  to  a 
concrete  floor.  An  adjustable  crank  attachment  is  fastened  to  one  of  the 
vertical  legs  of  the  "picture  frame".  A  hydraulic  actuator  piston  rod 
is  connected  to  this  crank,  while  the  cylinder  is  anchored  to  the  pedestal, 
to  provide  the  excitation. 

Dynamometers  attached  to  the  upper  and  lower  horizontal  members  of  the 
"picture  frame"  carry  a  vertical  support  near  each  side.  The  model  support 
frame  is  attached  to  these  vertical  supports  at  a  height  determined  by 
the  desired  model  vertical  location  relative  to  the  pitch  (roll)  axis. 

Thus,  all  loads  are  transmitted  into  the  "picture  frame"  through  the  dyna¬ 
mometers. 

An  overall  view  of  the  facility  is  shown  in  the  photograph  of 
Figure  A-l  and  clearance  dimensions  within  the  support  system  and  picture 
frame  are  shown  on  Figure  A. 2.  The  maximum  possible  angular  travel  is 
about  ±15°,  provided  the  actuator  force  capacity  is  not  exceeded. 

The  maximum  model  tank  width  and  height  are  fixed  by  the  clearances 
shown  in  Figure  A. 2.  Models  which  do  not  span  the  width  of  the  support 
frame  must  be  mounted  on  rails  attached  to  either  the  top  or  bottom  of  the 
support  frame.  The  actual  size  of  the  model  tank,  however,  must  be  such 
that  the  dynamometers  and  actuators  are  not  overloaded.  The  determining 
factors  in  selection  of  model  size  are;  (1)  The  total  mass  of  the  tank 
plus  liquid;  (2)  The  combination  of  excitation  frequency,  excitation  ampli¬ 
tude,  vertical  offset  of  the  model  (tank  plus  liquid),  center  of  gravity 
from  the  pitch  (roll)  axis,  and  the  horizontal  offset  of  model  from  the 
pitch  axis  which  results  in  the  greatest  inertia  moment  about  the  pitch  axis. 

These  variations  determiut  the  greatest  force  reactions  on  the  dynamometers. 

(3)  The  combination  of  the  above  parameters  which  result  in  the  highest 
reaction  force  on  the  actuator.  This  •  oobination  may  not  be  the  same  as 
that  which  results  in  the  greatest  dynamometer  reactions. 

A  maximum  model  size  cannot,  therefore,  be  determined  until  the 
entire  range  of  test  parameters  has  been  defined.  Reactions  on  the  dyna¬ 
mometers  and  actuator  must  be  investigated  for  every  combination  of  test 
parameters. 

Excitation  of  the  system  is  provided  by  a  linear  hydraulic  actuator 
connected  to  a  crank  on  the  "picture  frame".  The  actuator  has  a  maximum 
force  output  of  about  10,000  pounds  at  its  highest  velocity,  and  up  to 
14,000  pounds  at  very  low  velocities  (low  frequency  and  low  amplitude). 

A  self-contained  servo-valve  and  feedback  transducer  within  the  actuator  are 
used  in  conjunction  with  an  external  controller  to  control  the  actuator. 

The  controller  has  provisions  for  centering  the  actuator  at  any  point 
within  its  stroke  (at  the  center  of  the  stroke  for  normal  testing),  and 
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CROSS-SECTIONAL  DIMENSIONS  FIGURE  A  3  DRTVF  cvctcm  d,  ^ 

H.J.  DRIVE  system  block  diagram 


for  setting  the  excitation  amplitude.  Excitation  waveform  and  frequency 
are  determined  externally  to  the  controller  input.  A  function  generator 
connected  to  the  input  can  excite  sinusoidal,  square,  triangular,  or  sawtooth 
motions  of  the  actuator  at  any  discrete  frequency. 

A  tape  input  can  be  used  to  excite  random  motions  with  various  spectral 
contents,  or  any  other  specified  motion  (such  as  typical  ship  motions) 
within  the  limitations  of  the  actuator. 


The  actuator  is  limited  to  total  useful  stroke  of  about  2.2  inches 
peak-to-peak  and  a  maximum  linear  velocity  (restricted  by  the  servo-valve) 
of  about  8  inches  per  second.  The  limits  of  the  actuator  are  determined 
by  the  maximum  stroke  for  an  excitation  frequency  up  to  1.2  Hz,  then  by  the 
maximum  velocity  up  to  the  frequency  where  the  reaction  force  on  the 
actuator  rod  exceeds  the  maximum  actuator  force.  A  block  diagram  of  the 
drive  system  is  shown  on  Figure  A. 3. 

A. 2  Translational  Motion  Slosh  Rig 

The  translational  motion  slosh  rig  operates  on  a  system  very  similar* 
to  the  angular  motion  rig.  The  translational  machine  consists  of  a  5  foot 
by  5  foot  mounting  table  that  can  be  excited  with  simultaneous  vertical 
and  horizontal  motion  that  is  independent  along  each  axis.  Maximum  table 
payload  capacity  is  6000  pounds  dead  weight.  As  with  the  angular  motion 
unit,  the  translational  table  drive  mechanisms  are  servo-controlled  elec- 
trohydraulic  actuators,  and  have  the  following  capabilities: 


Frequency  Range 
Force  Capacity 
Maximum  Stroke 
Maximum  Velocity 
Maximum  Acceleration* 


Horizontal  Vertical 


0-200  Hz 
10,000  lb 
8.0  in. 
45  in. /sec 

16  g 


0-100  Hz 
20,000  lb 
6.6  in. 
22  In./sec 
10  g 


*(At  zero  payload) 


The  excitation  signals  are  provided  by  function  generators  or  analog 
tape  signals.  Table  displacement  is  accurately  controlled  by  automatic 
feedback  to  respond  to  an  arbitrary  voltage  signal. 
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APPENDIX  B 

PRESSURE-TIME  HISTORY  DATA  FOR  TRANSDUCER  LOCATIONS  2-13 

TABLE  B-l.  NOND IMENS IONAL  PRESSURE-TIME  HISTORY  VALUES 
AT  TRANSDUCER  LOCATION  2  FOR  A  25%  FULL  TANK 

% 

TRANSDUCER  LOCATION  2 
FLUID: WATER 

EXCITATION: F‘ITCH<AMP.  =+/-  4.800  DEG.) 

NUMBER  OF  CYCLES  PER  TEST:  £00 

PERCENT  FILLING:  £5 


PT.2 

PT.3 

PT.4 

PT.5 

AVG.KP 

6.6493 

1.0960 

1 . 0303 

0.0435 

MAX.KP 

15.6515 

£.1194 

1 . 88?3 

0. 94?4 

MIN. KP 

3.1??? 

8.235? 

0.3548 

0.002? 

AVG. TIME 

0. 1804 

0.3225 

1.4294 

1.9184 

MAX. TIME 

0.3??? 

0.5694 

1.6614 

2.2314 

MIN. TIME 

0.0398 

0. 1889 

1. 1393 

1.661? 

PT.l  IS  ALWAYS  ZERO. 

AIL  VALUES  PRINTED  IN  THE  TABLE  ARE  NOND I MENS IONAL 


NOTE:  See  Figure  IV-5  for  definitions. 
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NOTE:  See  Figure  IV-5  for  definitions.  NOTE:  See  Figure  IV-5  for  definitions 
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NOTE:  See  Figure  IV-5  for  definitions.  'NOTE:  See  Figure  IV-5  for  definitions. 


TABLE  B-12.  NONDIMENSIONAL  PRESSURE-TIME  HISTORY 
VALUES  AT  TRANSDUCER  LOCATION  13  FOR  A  75%  FULL  TANK 


TRANSDUCER  LOCATION  13 
FLU! D: WATER 

EXCITATION: PITCHCAMP. =+/-  4.000  DEG. > 
NUMBER  OF  CYCLES  PER  TEST:  200 
PERCENT  FILLING:  TO 


PT.2 

pr.3 

PT.4 

PT.5 

AVG.K’P 

5.6172 

0.0000 

0.0000 

0.0000 

MAX. KP 

15.3374 

0.0013 

0.0000 

0.0000 

MIN. KP 

2.5381 

0.0000 

0.0000 

0.0000 

AVG. TIME 

0. 1676 

0.4O66 

0.3276 

1.4816 

MAX. TIME 

O. 3303 

0.6182 

1.1873 

1.8528 

MIN. TIME 

O. 035S 

0.2842 

0.7128 

1.1405 

PT.l  IS  ALWAYS  ZERO, 

ALL  VALUES  PRINTED  IN  THE  TABLE  ARE  NONDIMENSIONAL. 


NOTE:  See  Figure  IV-5  for  definitions. 
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FIGURE  B-l.  NONDIMENSIONAL  FRESSURE- 
TIME  HISTORY  VALUES  FOR  200 
RESONANT  SLOSHING  CYCLES  AT  TRANS¬ 
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TANK 
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FIGURE  B-2.  INTEGRATED  NONDIMENSIONAL 
PRESSURE  VALUES  FOR  200  RESONANT 
SLOSHING  CYCLES  AT  TRANSDUCER 
LOCATION  4  FOR  A  25%  FULL 
TANK 
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FIGURE  B-3.  NONDIMENS IONAL 
PRESSURE  VS  IMPULSE  RISE  TIME 
FOR  200  RESONANT  SLOSHING 
CYCLES  AT  TRANSDUCER  LOCATION  4 
FOR  A  25*. FULL. -TANK 
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FIGURE  B-4.  NONDIMENSIONAL 
PRESSURE  VS  IMPULSE  DURATION 
FOR  200  RESONANT  SLOSHING 
CYCLES  AT  TRANSDUCER  LOCATION 
4  FC?  A  25*  FULL  TANK 
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FIGURE  B-5.  NONDIMENSIONAL 
IMPULSE  DURATION  VS  IMPULSE 
RISE  TIME  FOR  200  RESONANT 
SLOSHING  CYCLES  AT  TRANS¬ 
DUCER  LOCATION  4  FOR  A  25* 
FULL  TANK 


FIGURE  B-6.  NONDIMENSIONAL 
PRESSURE-TIME  HISTORY  VALUES 
FOR  200  RESONANT  SLOSHING 
CYCLES  AT  TRANSDUCER  LOCATION 
7  FOR  A  75*  FULL  TANK 
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FIGURE  B-7.  INTEGRATED  NUN¬ 
DIMENSIONAL  pRr^p;  VALUES 
FOR  200  RESONAfJ  SLASHING 
CYCLES  AT  TRANSC-CLR  LOCATION 
7  FOR  A  75%  FULL  TANK 
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FIGUkL  :  9.  NONDIMENS IONAL 
PRESSURE  VS  IMPULSE  DURATION 
FOR  200  RESONANT  SLOSHING 
CYCLES  AT  TRANSDUCER  LOCATION 
7  FOR  A  75%  FULL  TANK 
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FIGURE  B-8.  NONDIMENSIONAL 
PRESSURE  VS  IMPULSE  RISE  TIME 
FOR  200  RESONANT  SLOSHING 
CYCLES  AT  TRANSDUCER  LOCATION 
7  FOR  A  75%  FULL  TANK 
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FIGURE  B-10.  NONDIMENSIONAL 
IMPULSE  DURATION  VS  IMPULSE 
RISE  TIME  FOR  200  RESONANT 
SLOSHING  CYCLES  AT  TRANSDUCER 
LOCATION  7  FOR  A  75*  FULL  TANK 
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FIGURE  B-ll.  NONDIMENSIONAL 
PRESSURE-TIME  HISTORY  VALUES 
FOR  200  RESONANT  SLOSHING 
CYCLES  AT  TRANSDUCER 
LOCATION  11  FOR  A  75*  FULL 
TANK 
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FIGURE  B-12.  INTEGRATED  NON- 
DIMENSIONAL  PRESSURE  VALUES 
FOR  200  RESONANT  SLOSHING 
CYCLES  AT  TRANSDUCER  LOCATION 
11  FOR  A  75%  FULL  TANK 
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FIGURE  B-13.  NONDIMENSIONAL 
PRESSURE  VS  IMPULSE  RISE 
TIME  FOR  200  RESONANT 
SLOSHING  CYCLES  AT  TRANS¬ 
DUCER  LOCATION  11  FOR  A 
75%  FULL  TANK 


FIGURE  B-14.  NONDIMENSIONAL 
PRESSURE  VS  IMPULSE  DURATION 
FOR  200  RESONANT  SLOSHING 
CYCLES  AT  TRANSDUCER 
LOCATION  11  FOR  A  75%  FULL 
TANK 
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FIGURE  B-15.  NONDIMENSIONAL 
IMPULSE  DURATION  VS  IMPULSE 
RISE  TIME  FOR  200  RESONANT 
SLOSHING  CYCLES  AT  TRANS¬ 
DUCER  LOCATION  11  FOR  A  75 % 
FULL  TANK 
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FIGURE  B-17.  INTEGRATED  NON- 
DIMENSIONAL  PRESSURE  VALUES 
FOR  200  RESONANT  SLOSHING 
CYCLES  AT  TRANSDUCER  LOCATION 
13  FOR  A  75%  FULL  TANK 
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FIGURE. B-16.  NONDIMENSIONAL 
PRESSURE-TIME  HISTORY  VALUES 
FOR  200  RESONANT  SLOSHING 
CYCLES  AT  TRANSDUCER  LOCATION 
13  FOR  A  75%  FULL  TANK 
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FIGURE  B-18.  NONDIMENSIONAL 
PRESSURE  VS  IMPULSE  RISE  TIME 
FOR  200  RESONANT  SLOSHING 
CYCLES  AT  TRANSDUCER  LOCATION 
13  FOR  A  75%  FULL  TANK 
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FIGURE  B-19.  NONDIMENSIONAL 
PRESSURE  VS  IMPULSE 
DURATION  FOR  200  RESONANT 
SLOSHING  CYCLES  AT  TRANS¬ 
DUCER  LOCATION  13  FOR  A 
75*  FULL  TANK 
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FIGURE  B-20.  NONDIMENSIONAL 
IMPULSE  DURATION  VS 
IMPULSE  RISE  TIME  FOR  200 
RESONANT  SLOSHING  CYCLES 
AT  TRANSDUCER  LOCATION  13 
FOR  A  75*  FULL  TANK 
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APPENDIX  C 


ONE-DEGREE- OF-  FREEDOM  EQUIVALENT  SYSTEMS 

To  derive  one-degree~of-freedom  (one  dof)  "equivalent"  systems  for 
what  are  normally  regarded  as  multi-dof  structural  elements  requires  that 
some  assumptions  be  made  about  the  behavior  of  the  structure.  The  basic 
assumption  is  a  deformation  pattern  which  defines  displacements  of  the 
total  element  in  terms  of  a  single  displacement  at  some  point  in  the 
structure. 

As  an  example,  consider  the  simply  supported  beam  of  Figure  C-l. 

If  the  displacement  along  the  beam  can  be  described  in  terms  of  a  single 
variable,  for  example,  the  center  displacement,  wQ,  then  the  motions  are 
reduced  to  a  single  dof.  This  is  accomplished  by  assuming  a  deformation 
pattern  for  the  beam  which  is  some  function  of  w0.  Various  choices  can  be 
made,  e.g.,  the  fundamental  mode  shape,  the  static  deformed  shape  for  the 
load  distribution  of  the  dynamic  loading,  or  simply  some  approximate  shape 
which  resembles  the  fundamental  mode  or  static  deformed  shape  and  matches 
the  appropriate  boundary  conditions.  For  elastic  behavior,  Biggs (79) 
chooses  the  static  deformed  shape  of  the  structure.  For  the  simply 
supported  beam,  uniformly  loaded  by  the  distributed  load  N*g(t)»  [g(t)J  is 
the  time  function  for  the  load  ;  the  static  deformed  shape  for  elastic 
behavior  is  given  in  Figure  C-l. 


N-g(t) 


{III  H" 
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FIGURE  C-l.  DEFORMATION  PATTERN  FOR  SIMPLE  BEAM 


Reducing  displacements  to  a  single  dof  does  not  guarantee  similarity 
between  the  computed  response  of  a  one-dof  system  and  the  structural 
element.  One  way  that  similarity  can  be  obtained  is  to  equate  the  energies 
of  the  distributed  and  single-dof  systems.  At  any  time » the  internal  strain 
energy  and  the  system  kinetic  energy  must  equal  the  external  work  as  shown 
by  Equation  C-l. 


W  =  U  +  KE 


(C-l) 


C-l 


W  =  external  work 


U  *  internal  strain  energy 
KE  =  kinetic  energy 

If  kinematic  similarity  (equal  displacements  and  velocities  in  this  case) 
is  maintained  between  the  distributed  and  single-dof  systems,  then  equating 
energies  of  the  two  systems  will  assure  similarity  of  the  computed  behavior. 
This  is  true  for  elastic  or  elastic-plastic  behavior,  but  in  ship  tanks 
only  elastic  behavior  will  be  considered. 

As  a  first  example,  we  chose  the  fundamental  mode  of  the  simply 
supported  beam  as  the  deformed  shape.  This  choice  will  permit  a  check  on 
frequency  with  the  exact  solution  to  demonstrate  kinematic  equivalence  of 
the  distributed  and  equivalent  systems.  The  deformed  shape  for  the  funda¬ 
mental  mode  is 


w(x)  -  wrt  sin  r— 

O  L 


(C-2) 


Evaluating  W,  U,  and  KE  for  this  deformation  pattern  gives 


L 

W  ■  /  N*g(t)  w(x)  dx  «  -  N«g(t)  w  L 

*  71  O 

0 


where 

E  *  material  modulus 
I  =  beam  section  moment  of  inertia 


(C-3a) 


(C-3b) 


(C-3c) 


A  =  beam  cross-sectional  area 


p  =  material  density 

wq  =  the  beam  center  displacement  and  is  some  function  of  time. 

These  same  energies,  evaluated  for  the  single-dof  system  of  Figure  C-2, are 


FIGURE  C-2.  EQUIVALENT  ONE-DOF  SYSTEM 

Equating  equations  C-3  and  C-4  and  recognizing  that 

pAL  =  m,  total  beam  mass 
NL  =  F,  total  beam  static  load 


the  following  relationships  are  obtained 


F 


kl  = 

~  *  0.6366 

(C-5a) 

k  = 
e 

irAEI  n  , 

-  *  0.634  k 

(C-5b) 

2L 

m 

K  « 

m 

—  =  0.500 
in 

(C-5c) 

where 

k  ■  spring  constant  for  the  simply  supported  beam 

(N*L  divided  by  static  center  deflection) 

^  ■  load  factor 

K  ■  mass  factor 
m 

Notice  that  the  equivalent  spring  constant  is  very  nearly  equal  to  the 
static  spring  constant  multiplied  by  the  load  factor. 

Because  these  relationships  were  developed  on  the  basis  of  kinematic 
equivalance  of  the  displacement  v0,  it  is  instructive  to  compare  the  fre¬ 
quency  given  by  the  equivalent  mass  and  spring  constants  with  the  exact 
solution  for  a  simply  supported  beam.  The  exact  solution  for  the  funda¬ 
mental  mode  is  given  by  Den  Hartog(95)  a8 


2  El 

“  *  —3  (C— 6) 

>  mL 


where 

m  *  total  beam  mass  =  p  AL 
For  the  equivalent  system, 


1  » 


As  expected,  the  assumption  of  the  fundamental  mode  for  the  deformed 
shape  gives  a  frequency  which  is  exactly  equal  to  the  beam's  fundamental 
frequency. 

Now,  K,,  1^,  and  ke  will  be  determined  for  the  static  deformed  shape 
of  the  beam  for  comparison  with  the  results  of  Biggs  0$).  The  static 
deformed  shape  is  given  in  Figure  C-l,  and  for  this  case  the  energies  of 
the  distributed  system  are: 


W  =  N’g(t)  - - — 

(C-8a) 

3072  El  w  2 

U  -  0 

U  3 

1251/ 

(C-8b) 

1984  .  2 

^  "  7875  pAL  Wo 

(C-8c) 

Equating  to  equations 

(C-4),  as  before,  we  obtain 

F 

K L  -  f  '  »■« 

(C-9a) 

,  6144  El  „  ,,  . 

k  «* - r —  -  0.64  k 

1251/ 

(C-9b) 

m 

K  -  —  -  0.5026  -  0.50 
m  m 

(C-9c) 

These  values  of  stiffness  and  mass  of  the  equivalent 

of 

system  give  a  frequency 

w  =  9.889  /— s* 

7»l3 

(C-10) 

which  is  within  0.2%  of  the  fundamental  mode  of  the  beam  [Equation  (C— 6) ] . 
Because  the  static  deformed  shaps  is  made  up  of  many  normal  modes,  in  addi¬ 
tion  to  the  first,  some  difference  in  frequency  is  to  be  expected.  If  the 
beam  truly  deforms  dynamically  with  the  static  deformed  shape,  then  it  would 
respond  with  the  equivalent  system  frequency. 
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Several  observations  can  be  made  as  a  result  of  these  calculations 
for  elastic  behavior: 

1.  A  rational  method  is  available  for  deriving  an  equivalent 
one-dof  system  for  a  distributed  structure. 

2.  Kinematic  equivalency  is  maintained  if  the  assumed  deformation 
pattern  exactly  matches  the  behavior  of  the  distributed  system. 

3.  The  values  of  K.,  Kjj,  and  ke  match  those  given  by  Biggs^) 
and  others  for  the  static  deformed  shape.  Table  C-l  gives 
these  values  for  simply  supported  beams. 

4.  For  the  static-deformed  shape,  the  equivalent  stiffness,  ke, 
can  be  computed  from  the  static  stiffness  of  the  distributed 
system  and  the  load  factor. 

5.  The  fundamental  mode  shape  and  the  static  deformed  shape  give 
very  similar  results  for  a  uniform  loading.  Using  the  static 
deformed  shape  should  give  the  best  results  for  other  than 
uniform  loading  or  uniform  structures  because  the  deflected 
shape  will  not  match  the  fundamental  mode  so  closely. 

Transformation  factors  as  in  Table  C-l  have  been  developed  by  Biggs (79) 
and  others  for  many  different  structural  elements  and  different  loading 
conditions  for  each  element.  For  example,  Biggs  gives  values  for 

o  simply-supported  beams 

O  clamped  beams 

o  propped  cantilevers 

o  one-  and  two-way  slabs  with  edge  supports 

o  two-way  slabs  with  interior  supports 

C.l  Displacements  of  the  Equivalent  Systems 

Displacements  of  these  equivalent  one -'degree -of -freedom  systems  can 
be  computed  using  simple  numerical  methods  or  closed  form  solutions.  A 
convenient  way  to  present  the  solution  for  elastic  behavior  is  in  terms  of 
a  dynamic  load  factor  (DLF).  The  dynamic  load  factor  is  simply  the  maximum 
displacement  produced  by  the  dynamic  load,  divided  by  the  displacement 
produced  by  the  peak  Jocid  applied  statically.  Thus,  the  DLF  can  be  either 
smaller  or  greater  than  one  depending  upon  the  nature  of  the  loading  and 
vibration  period  of  the  responding  structure. 

Plots  of  dynamic  load  factors  for  different  types  of  loading  functions 
are  given  in  Figure  C-3.  Note  from  these  figures  that  the  maximum  DLF  is 
2.0,  whereas  the  minimum  approaches  zero  for  certain  types  of  loads  of 
very  short  duration. 


As  Biggs (79)  states, 

"It  is  important  to  recognize  that  the  dynamic  reactions  of 
the  real  structural  element  have  no  direct  counterpart  in  the 
equivalent  one-degree  system,  i.e.,  the  spring  force,  is  not  the 
same  as  the  real  reaction.  This  is  true  because  the  simplified 
system  was  deliberately  selected  so  as  to  have  the  same  dynamic 
deflection  as  the  real  element,  rather  than  the  same  force  or 
stress  characteristics." 

These  comraei^s  are  particularly  true  for  elastic-plastic  behavior,  where 
the  maximum  resistance  is  limited  by  permitting  plastic  deformation  to 
occur  in  the  structure.  For  this  case  the  maximum  resistance  is  not  a 
good  indication1  of  the  maximum  shear  which  has  occurred  in  the  structural 
element.  When  only  elastic  behavior  occurs,  the  maximum  resistance  of  the 
member,  the  applied  load  multiplied  by  DLF,  gives  a  good  indication  of  • 
the  peak  shear  in  the  member.  Here,  only  elastic  behavior  will  be  treated 
in  keeping  with  the  philosophy  that  only  elastic  behavior  will  be  tolerated 
in  ship  tanks. 

To  determine  the  dynamic  reaction,  V,  of  Table  C-l,  a  free-body 
diagram  of  the  element,  including  inertia  effects,  must  be  used.  For  the 
uniformly-loaded,  simply-supported  beam  of  Figure  C-l,  such  a  free-body 
diagram  is  given  in  Figure  C-4.  Inertia  forces  are  distributed  according 
to  the  assumed  deformation  pattern  and  are  opposed  to  the  motion.  Only  ' 
the  period  of  time  up  to  the  maximum  response,  t,,,,  is  considered. 

From  the  free-body  diagram  of  one-half  of  the  beam,  an  expression  for 
V(t)  is  obtained  by  taking  moments  at  x  -  a,  the  line  of  action  of  the 
resulting  inertia  forces. 

V(a)  -  \mh/2  -  F(t)  a  -  |  -  0  (C-ll) 

For  elastic  behavior,  locations  of  the  resultant  of  the  inertia  forces  and 
the  mid-span  moment  are  determined  from  the  static  deformed  shape.  They  are 


FIGURE  C-4.  SIMPLY  SUPPORTED  BEAM  IN  EQUILIBRIUM  UNDER  DYNAMIC 
LOADING 


Substituting  into  Equation  (C-ll),  the  dynamic  reaction  at  the  beam  end 
for  elastic  behavior  is 

v  B  .30 6 El  w  +  o.!066  F(t)  (C-12) 

L3  ° 

where  F(t)  is  the  total  dynamic  load  on  the  beam  [N*g(t)*L]  and  wQ  is 
understood  to  be  a  function  of  time.  Thus,  the  shear  reaction  is  a  function 
of  the  loading  as  well  as  the  elastic  resistance  in  the  member.  Notice 
in  Equation  C-12  that  the  maximum  reaction  may  occur  when  wD  is  a  maximum, 
when  F(t)  is  a  maximum  or  somewhere  in  between. 

For  elastic  behavior  only,  it  is  worthwhile  to  compare  the  shear 
reactions  predicted  by  Equation  C-12  or  the  equation  in  Table  C-l  with  that 
obtained  from  the  equivalent  static  loading.  This  can  be  easily  done  for 
a  simply  supported  beam,  uniformly  loaded  by  a  triangular  force  with  zero 
rise  time.  This  loading  most  closely  approximates  the  worst  case  sloshing 
load  because  the  rise  time  is  only  about  1.0  ms  and  is  very  short  relative 
to  the  decay  time. 

The  shear  equation  from  Table  C-l  is 

V  =  ,39R  +  .llF(t)  (C-13) 
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For  elastic  behavior  the  resistance  at  the  time  of  maximum  response,  t^, 
will  equal  the  dynamic  load  factor,  DLF,  multiplied  by  the  peak  applied 
load  Fj.  This  gives 


V  -  .39*DLF*F1  +  .11  F(t) 


(014) 


Also,  for  a  triangular  pulse  of  duration  t^,  F(t)  can  be  expressed  in  terms 


of  the  peak  load  F^  as 


<td  -  t) 

F(t)  -  F,  - -  for  o  <  t  <  t. 

!  td  d 


(C-15) 


t,-t 

V  -  .39*DLF*F.  +  ,11  F  -~- 
1  1  cd 


Fx  [.39DLF  +  .ll(td-t)/tdJ 


(C-16) 


At  t  ■  0,the  second  term  in  brackets  is  maximum  but  the  first  is  zero  (the 
system  has  not  displaced).  The  first  term  in  the  brackets  corresponds  to 
the  maximum  resistance  and  occurs  at  t  ■  t  .  Because  sloshing  pressures 
are  of  long  duration,  the  shear  will  most  Tikely  peak  at  t  ■  t  .  Thus, 


-  .39‘DLF  4  ,ll(t .  -  tj/t. 

a  m  a 


(C-17) 


This  expression  is  compared  to 


Vf  -  O.Stf^DLF) 


(018) 


which  is  the  maximum  shear  in  a  dimply  supported  beam  acted  on  by  a  uniform 
static  load  F^'DLF.  In  non-dimensional  form 


0.5'DLF 


(019) 


Equating  equation  (019)  to  equation  (0-17)  and  setting  the  second  term  in 
equation  (017)  to  its  maximum  value  of  0.11  we  find 


.5DLF  «  .39* DLF  +  .11 


(020) 


The  sides  are  only  equal  when  DLF  is  1.0  (static  loading).  When  DLF  is 
larger  than  1.0,  the  left  side  is  greater  because  the  term  0.11  is  no*-  multi¬ 
plied  by  DLF.  For  DLF  <  1.0  the  right  side  can  be  greater,  depending  upon 
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the  contribution  of  the  second  term  for  which  0.11  is  the  maximum  value. 
Because  sloshing  pressures  have  short  rise  times  and  long  duration  (rela¬ 
tive  to  the  structural  periods  of  interest),  the  DLF  will  always  be  greater 
than  1.0  for  LNG  sloshing  forces  in  ship  tanks. 


C-12 


APPENDIX  P 


EXAMPLE  CALCULATIONS  FOR  MEMBRANE  AND  PRISMATIC  TANKS 

Calculations  were  performed  for  typical  LNG  tank  structures  to  illus¬ 
trate  the  methods  of  computing  structural  frequencies  and  the  dynamic 
load  factors.  In  addition,  calculations  of  structural  strength  have  been 
made  and  are  expressed  in  terms  of  allowable  sloshing  pressures  which  the 
structural  components  can  withstand.  Allowable  sloshing  pressures  for  the 
Gaz  Transport  insulation  boxes  are  compared  with  results  from  laboratory 
tests  conducted  on  the  boxes.  Results  of  the  laboratory  tests,  in  which 
both  static  and  dynamic  pressure  loads  were  applied  to  the  cover  of  the 
boxes,  are  summarized  in  Section  IV. 5.  It  should  be  noted  that  the  structural 
strength  calculations  are  only  examples  of  how  such  analyses  may  be 
performed.  Certainly,  more  rigorous  methods  are  available,  as  discussed 
in  Section  V.2.3,  and  might  be  preferred;  however,  for  these  cases  the 
one  dof  analysis  is  a  good  approximation  to  actual  structural  behavior  and 
should  give  conservative  results. 

D.l  Plywood  Insulation  Boxes 

In  this  example, the  strength  of  the  cover  and  longitudinal  stiffened 
in  the  12-cell  box  will  be  determined.  These  are  the  primary  load-carrying 
members  in  the  primary  boxes  and  failures  were  produced  in  these  members 
during  laboratory  testing.  Dynamic  load  factors  for  the  boxes  are  based  on 
a  tank  which  is  36m  (118  ft)  long.  Tank  length  sets  the  time  constants 
for  the  sloshing  pressures  and  determines  the  dynamic  amplification  as 
shown  in  Section  V.2.1. 


Material  Properties 

The  insulation  boxes  are  constructed  of  7-ply  birch  plywood.  Flexural 
strength  and  stiffness  were  measured  in  3-point  bending  tests  both  at  room 
and  at  LN2  temperatures.  Exterior  plies  were  oriented  perpendicular  (cross-grain) 
to  the  direction  of  the  bending  stress;  this  orientation  corresponds  to  the 
principal  loading  direction  of  the  box  covers.  Results  from  these  tests 
are  summarized  in  Section  IV. 4. 

Material  properties  used  in  the  analysis  are  given  in  Table  D-l. 

Values  for  the  elastic  modulus  and  rupture  stress  in  the  perpendicular  (1) 
direction  are  taken  from  Section  IV. 6.  The  elastic  modulus  parallel  (||) 
to  the  grain  was  calculated  using  properties  for  white  birch  from  the 
Wood  Handbook^^)  and  procedures  outlined  in  the  Plywood  Design  Guide(°8)t 
Plywood  density  was  measured  on  box  material  and  Poisson’s  ratios  parallel 
and  perpendicular  to  the  grain  of  the  exterior  ply  were  estimated  using 
values  for  birch  from  Reference  97 . 

Cover  Strength 

The  corner  cells  are  the  weakest  part  of  the  box  covers  because,  for 
uniform  pressure  loadings,  the  boundaries  over  the  stiffeners  are  treated 
as  clamped  and  the  edges  along  the  side  and  end  of  the  box  are  treated  as 
simply  supported.  All  other  cells  have  three  clamped  edges.  Dimensions 
a»,d  boundary  conditions  of  the  corner  cell  are: 
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TABLE  D-l.  PLYWOOD  MATERIAL  PROPERTIES 


ar  N/mm2 (psi) 

Ejj  N/mm2(psi) 

V 

p  kg/nm3(ii-S|si) 
in. 

Ea  N/nm2(psi) 


New  Plywood  -  Average  Properties  -  RT 
Direction  With  Respect  to  Grain  of  Exterior  Ply 
Perpendicular  (y)  Parallel  (x) 


55.96  (8,117) 
3070  (445,304) 
.206 


108.4  (15,727) 
7149  (1,037,000) 
.268 


55.36  x  10"8  (5.18  x  10~5) 
464  (67,260)  618  (89,680) 


FIGURE  D-l 


Static  Pressure 


Timoshenko  and  Woinowsky-Krieger  (  99  )  give  a  solution  for  bending 
stresses  in  a  plate  with  these  boundary  conditions e  The  maximum  moment 
occurs  near  the  center  of  side  b  and  is  given  by 

My  -  0.1180  Pb2  for  a/b  «  2  (D-l) 

For  the  corner  cell 

a/b  ■  3.47, 

but  the  magnitude  .f  the  momenc  changes  very  little  for  b/a  >  2  so  the 
value  given  is  close  to  the  true  value.  The  stress  produced  by  this  moment 
is 

6M 

o  *  (D-2) 

hJ 


Substituting  Equation  (D-l)  for  My  and  the  rupture  stress  Ov  for  o,  an 
expression  for  the  allowable  static  pressure  pa,  can  be  obtained. 


p  °rh 

p  ss  - 

ra  2 

,708b^ 


(D-3) 


For  or  =  55.96  N/tnm2  (8117  psi) 
h  =  9  mm  (.354  in.) 
b  =  88  mm  (3.46  in.) 


Equation  (D-3)  gives 
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P  =  0.827  N/mm2  (120  psi) 

d 

2 

Thus,  0.827  N/mm  (120  psi)  is  the  allowable  static  pressure  of  new  boxes  at 
room  temperature  based  upon  average  values.  This  gives  a  50%  probability  of 
failure  if  the  properties  have  a  normal  distribution. 

Test  results  on  three  new  12-cell  boxes  produced  failures  at  static 
pressures  of  120  psi,  130  psi  and  150  psi.  Table  IV-22  in  Section  IV. 5 
shows  the  damage  and  describes  the  failure.  In  three  tests  failures 
occurred  in  the  corner  cells.  In  one  test»the  outer  cell  in  the  center 
of  the  long  cell  also  failed.  This  is  the  second  weakest  cell  in  the  box 
cover.  The  damage  is  cited  as  shear  failure  because  the  cover  was  completely 
severed  at  the  stiffener.  Shear  properties  were  not  measured  for  the  covers 
and  these  failures  were  not  predicted.  Rough  estimates  indicate  that  the 
cover  has  about  equal  shear  and  bending  strength,  and  the  resulting  failure 
is  some  combination  of  shear,  bending,  and  local  compression.  Compression 
tends  to  cut  through  the  inner  ply;  tension  tends  to  crack  the  outer  ply. 

These  together  weaken  the  plywood  in  shear. 

Dynamic  Pressure 

To  determine  the  permissible  dynamic  or  sloshing  pressure,  rise  and 
decay  times  of  the  loading  rauot  be  known.  These  times  are  a  function  of 
tank  length  and  were  determined  in  Section  V.2.1  for  a  tank  36  m  (118  ft) 
long.  From  these  times, the  envelope  of  the  dynamic  load  factor  (DLF)  was 
also  determined  and  is  giver,  in  Figures  V-2  and  V-3.  Thus,  these  calculations 
will  be  for  a  tank  which  is  36m  long  so  that  the  DLF  already  calculated 
can  be  used. 

To  read  the  value  of  DLF  from  Figure  V-2  and  V-3, the  fundamental 
frequency  of  the  structural  component  must  be  computed.  Equations  (V-l 
and  (V-2)  of  Section  V.2.2.1  will  be  used.  Substituting  plywood  properties 
from  Table  D-l  and  the  plywood  thickness  from  Figure  D-l,we  have 

y>3 

Dx  "  12(l-'v"v  )"  458’426  N,mra  (4058  lb  ln-> 

'  "xy 

E  h3 

Dy  -  12(1.V  v  )m  196»837  N,tn®  <1742  lb  ln‘> 

x  y 

rh3  *  • 

D  *  D  V  +  -r— —  ~  D  /4.8  *  95,505  N-mm  (845  lb  in.) 
xy  x  y  6  x 


The  frequency  equation  is 


From  Table  V.l  A  =  B  =  3.927 
C  =  132.118 


*  Measured  by  Herman  (100)  on  five-ply  maple  plywood. 
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Substituting  the  known  parameters,  the  fundamental  frequency  and  period  of 
the  corner  cell  in  the  cover  are  found  to  be 

f  -  2058  Hz;  T  -  -  -  0.485  ms 


These  values  are  for  the  box  cover  vibrating  in  air,  whereas  in  the  tank  it 
is  exposed  to  LNG  on  one  side  (even  though  separated  by  the  membrane).  The 
effect  of  the  membrane  on  the  added  mass  of  the  liquid  is  uncertain; 
however,  including  the  effect  of  the  liquid  on  the  frequency  will  increase 
the  period  and  produce  a  higher  DLF  (see  Figure  V-3),  a  conservative 
effect.  Thus*the  added  mass  will  be  included.  Section  V.2.2.5  gives 
approximate  formulas  which  account  for  submergence.  The  frequency  in  air, 
fa,  is  multiplied  by  a  parameter  ^  to  find  the  frequency,  f^,  in  the  liquid 
ip  is  given  by 


r 


K 


♦-  «75 


x  £  *  'p 


where 


For  LNG 


Pi  "  *48pH,0  "  *48  X  10 


r6  K  ! 4.63  x  I<f5 


nun 


in. 


Substituting  for  p.  and  the  panel  dimensions  from  Figure  D-l  we  have: 


K  -  .3346  f  -  \pf  =  1076  Hz 

e  a 

il>  =  .523  T  «  1/f  =  0.928  ms 

e 


From  Figure  V-3  the  dynamic  load  factor  is  found  to  be 


DLF  =  1.22 


To  predict  the  allowable  sloshing  pressure  for  the  box  cover,  the  allowable 
static  pressure  is  divided  by  this  factor;  however,  to  be  representative 
of  in-service  conditions,  the  allowable  static  pressure  of  .827  N/mm^ 

(120  psi) ,  previously  calculated,  should  be  adjusted  to  account  for  LKG 


temperatures  and  in-service  effects.  The  average  strength  of  the  plywood 
at  LNG  temperatures  is  increased  as  shown  in  Table  IV- 26  from  55.96  N/mm2 
(8117  psi)  to  89.05  N/mm2  (12,916  psi).  Section  IV-4  also  gives  the 
degradation  of  strength  for  in-service  conditions  as  79%.  Accounting 
for  the  three  effects  of  dynamic  loading,  LNG  temperatures  and  in-service 
conditions,  the  allowable  sloshing  pressure  for  the  cover  is 


N 


.827 

1.22 


mm 


(ffritX0*79) "  *853  ps1> 


Average  rather  than  minimum  properties  have  been  used  in  these  calculations, 
so  there  is  a  50%  probability  of  a  bending  failure  at  this  pressure  for  a 
normal  distribution  of  plywood  properties. 

Strength  of  Longitudinal  Stiffener 

From  the  analysis  of  the  corner  cell  in  the  box  cover,  it  is  apparent 
that  the  majority  of  the  pressure  loading  on  the  cover  is  transferred  to 
■the  long  edges  of  the  cell.  These  edges  are  supported  by  the  sides  of  the 
box  or  by  the  longitudinal  stiffeners.  The  longitudinal  stiffener  picks  up 
load  from  two  adjacent  cells  and  so  is  more  highly  loaded  than  the  box  sides. 
Figure  III-  19  shows  the  arrangement  of  the  stiffeners  in  the  box  and  their 
support  points  on  the  transverse  stiffeners  of  the  secondary  boxes. 

* 

Figure  D-2  shows  the  longitudinal  stiffener  at  one  end  of  the  primary 
box.  A  uniform  loading  N  is  applied  to  the  top  of  the  stiffener  from  the 


305  nun 
(12  in.) 


ms 


FREE  EDGE 
(SLOT  FOR 
TRANSVERSE 
STIFFENER) 


EDGES  SIMPLY  SUPPORTED 
AT  TOP, END  AND  BOTTOM 
OF  BOX 


182  mm 
(7.17  in.) 


-SUPPORT  FROM 
TRANSVERSE  STIFFENER 
IN  SECONDARY  BOX 


FIGURE  D-2 


cover  and  the  bottom  is  assumed  to  be  supported  only  at  the  transverse 
stiffener  in  the  secondary  box.  Simple  solutions  do  not  exist  for  this 
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combination  of  geometry,  support  conditions  and  loading.  A  rather  domplex 
multidegree-of-freedom  analysis,  such  as  with  finite  elements,  would  be 
required  for  a  good  representation  of  these  conditions;  however,  some 
approximations  can  be  made  to  bound  the  buckling  strength  of  the  stiffener 
and  this  approach  will  be  followed  here. 

Static  Pressure 


Two  approximations  to  the  actual  conditions  were  checked.  They  are 
shown  schematically  in  Figure  D-3. 


N 

N  ♦  b 

rrm 

1 

b 

b 

a 

a 

♦  ♦  *  ♦ 

t 

N 

N  •  b 

(a)  UPPER  BOUND 

(b)  LOWER  BOUND 

FIGURE  D-3 


Solutions  for  these  cases  are  given  by  Roark  and  Young  (101).  Noting 
that  the  load  N  is  equal  to  the  pressure  applied  to  the  box  cover 
multiplied  by  the  width  of  a  cell  (N  ■  88mm*P)  we  can  express  the  allowable 
buckling  .load  directly  in  terras  of  the  pressure  on  the  box  cover.  Material 
properties  for  the  load  perpendicular  to  the  grain  of  the  exterior  ply 
(the  weakest  direction)  have  been  used. 

2 

Case  (a)  N  ■  (88mm)P  ■  4.27  —  -r  f •  h  (D-4) 

3  i-v2'b' 


For  the  load  perpendicular  to  grain  of  exterior  ply 

2 

-(-2L) 

>  i  / 

(.206  / 


pa  -  »'-2 


Case  (b)  N*b  »  (88mm)  Pa*b  =*  j  — E-  b 


(1-V2)b 


(D-5) 


Again  for  loading  perpendicular  to  the  grain  of  the  exterior  ply  the 


D-7 


allowable  pressure  Is 

P  *  .30  N/ram^  (43  psi) 

These  two  cases  should  give  reasonable  bounds  for  the  buckling 
strength  of  the  longitudinal  stiffeners,  expressed  in  terms  of  a  uniform 
static  pressure  applied  to  the  box  cover.  Of  course,  average  properties 
have  been  used  and  the  grain  of  the  stiffener  is  assumed  to  be  parallel 
to  the  box  axis  and  perpendicular  to  the  load  direction. 

In  laboratory  tests  on  the  boxes  (refer  to  Table  IV- 22),  buckling 
of  the  longitudinal  stiffener  in  the  end  cell  was  observed  at  a  static 
pressure  of  1.03-1.24  N/mnr  (150-180  psi).  The  exterior  grain  of  the  stif¬ 
fener  was  parallel  to  the  box  axis,  corresponding  to  the  direction  in  the 
predictions  above.  Thus,  the  upper  bound  estimate.  Case  (a)  appears  to 
offer  a  reasonable  approximation  to  the  behavior  of  the  stiffener.  It  is 
interesting  to  note  in  Table  IV-  22  that  on  Test  No.  7  the  static  pressure 
reached  1.65  N/mm^)  (240  psi)  without  stiffener  failure.  Stiffeners  for 
the  box  used  in  this  test  were  oriented  with  the  exterior  grain  perpendicular 
to  the  box  axis  (the  strong  direction).  The  predicted  buckling  pressure 
for  this  orientation  would  Increase  by  the  ratio  of  the  elastic  modulus 
in  the  two  directions  or  would  be  approximately  doubled. 

Dynamic  Pressure 

Under  dynamic  conditions , the  response  of  the  cover  was  increased  by  a 
factor  of  1.22.  Depending  upon  the  frequency  of  the  stiffener  in  compression, 
similar  amplification  may  occur.  To  estimate  the  frequency  of  the  stiffener 
for  compressive  loads , use  a  spring  mass  model  as  in  the  sketch  below.  Den 
Hartog  (  95  )  gives  the  frequency  equation  for  this  case  as 


M  «  M  +  M 
c  e 


M  *  mass  of  cover 
c 

Mg  *  added  mass  of  the  liquid 
m  ■  mass  of  stiffener 

K  =  axial  spring  constant  of  the  stiffener 
For  a  unit  width  of  the  stiffener  and  cover  the  values  are: 


M 


c 


(88mm) (1) (9mm) (55.36  x  10  ®)  =  4.385  x  10  ^  kg 

(2.498  x  10“6  1*>~se--2) 

In. 


B 


iKt* 
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M  *  -  I  •  (1) (49.48  x  10~8)  -  1.51  x  10~3kg 

*  2  ^  2 
/ ^  t  rt"* ^  lb-sec  \ 

(8*572  x  10  — jjj - ) 

m  -  (1)  (182)  (9)  (55.36  x  10"8)  «  9.068  x  loAg 
(S.17  x  l<f6 

in. 

K  -  AE/L  -  (9) (1) (464)/182  -  22.94  N/n a  (131  lb/in.) 


Substituting  into  the  frequency  equation  gives 


f  -  509  Hz;  T  ■  1/f  15  1.96  ms 

This  period  is  approximately  twice  that  of  the  cover.  Generally,  for  this 
condition  (a  large  spread  in  the  stiffnesses  of  a  coupled  system)  the  load 
transferred  to  the  supporting  system  will  be  attenuated;  however,  for  this 
approximation  it  was  assumed  that  the  load  is  unaltered.  From  Figure  V-2 
the  dynamic  load  factor  for  this  frequency  is 


DLF  -  1.63 


As  for  the  cover,  to  predict  the  allowable  buckling  load  produced 
by  sloshing  pressures,  the  increase  in  stiffness  which  occurs  at  reduced 
temperatures  and  the  reduction  in  strength  for  in-service  conditions  should 
also  be  accounted  for.  From  Section  IV. 6,  Table  IV-26,  we  find 


4785 

3070 


and  the  in-service  reduction  is  79%  as  before.  Combining  these  factors,  the 
allowable  sloshing  pressure  for  dynamic  buckling  (using  measured  buckling 
load  from  Test  10  in  Table  IV-22)  is 

P  ■  (1.03  -  1.24  N/mm2)  (^)(0,79) 
a  1.63 

»  .778  -  .937  N/mm2 (113  -  136  psi) 


*M£  is  taken  as  the  mass  per  unit  length  of  a  liquid  which  is  contained 
within  a  half  cylinder  which  lies  above  one  cell  of  the  box  cover. 


This  allowable  sloshing  pressure  for  buckling  of  the  longitudinal  stiffener 
is  based  upon  one  measured  value  of  the  static  buckling  load;  however,  the 
experiment  did  correlate  well  with  the  upper  bound  on  buckling  strength  com¬ 
puted  for  the  stiffener  using  average  material  properties.  As  for  the  cover, 
the  material  is  assumed  to  be  oriented  in  its  weakest  direction,  and  this 
condition  does  exist  in  the  boxes  tested. 


Technigaz  Membrane  Support  Svstem 


The  Technigaz  membrane  system  is  described  in  Section  III. 3.1  and 
shown  in  Figure  III. 3. 4.  Structural  support  for  the  membrane  is  provided 
by  a  composite  panel  formed  from  plywood  sheets  with  a  balsa  core.  The 
panel  is  separated  from  the  membrane  by  an  "end  grain"  balsa  pad  and 
supported  from  the  inner  hull  by  grounds.  Sloshing  pressures  on  the  membrane 
produce  bending  and  shear  stresses  in  the  composite  panel  between  the 
grounds.  The  allowable  sloshing  pressures  will  be  based  upon  these  stresses. 


A  section  of  the  composite  panel  between  the  secondary  grounds  is 
shown  in  Figure  D-4. 


SECONDARY  GROUND 
(TYPICAL) 


FIGURE  0-4 


Material  properties  for  balsa,  sweet  birch  and  sweet  sugar  maple  were  taken 
from  Reference  97.  Pertinent  properties  are  given  in  Table  D-2.  Details 
cf  the  plywood  construction  were  unknown  so  the  following  assumptions  were 
.Bade : 

(1)  Birch  plywood  is  7  ply  with  grain  of  exterior  plies  parallel 
to  the  grounds  (weakest  orientation). 

(2)  Maple  plywood  is  5  ply  with  grain  of  exterior  plies  parallel 
to  the  grounds. 


TABLE  D-2.  MATERIAL  PROPERTIES  FOR  BALSA. 
SWEET  BIRCH,  AND  SUGAR:  MAPLE 


Specific  gravity 

•* 

Elastic  modulus, 
longitudinal  dir,  EL 

Modulus  of  rupture 

Compression,  parallel 
to  grain 

Elastic  modulus,  transverse 
(ET/EL) 

Shear  parallel  to  grain 


Sweet  Sugar  i: 

Mlsa  Birch  Maple 

•17  .65  .63 

3792  N/mm2  14,962  N/ram2  12,618  S/m2 
(.55x10®  psi)  (2.17x10®  psi)  (1.83  x  10®  psi) 

117  N/mm2 
(16,900  psi) 

54  N/mm2 
(7,830  psi) 

015  .050  .050  (est.) 


2.07  I?/mm^ 
(300  psi) 


Shear  modulus,  (Gjjj/E) 
Longitudinal-radial  plane 


.054 


Based  upon  the  above  assumptions  and  material  properties  in  Table  D-2*the 
bending  stiffness  of  the  composite  panel  can  be  calculated.  For  a  25.4  mm 
wide  strip  between  the  grounds  the  axial  stiffnesses  of  the  plywood  faces 
were  found  to  be: 

birch:  (EA)b  -  (14962  N/mm2)(12  mm) (25.4  mm) (3/7  +  .050*4/7) 

-  2.085  x  106  N  (.469  x  106  lb) 


maple:  (EA)m  -  12,618  (3.4) (25.4) (2/5  +  .050*3/5) 
-  .4686  x  106  N  (.105  x  106  lb) 


Neglecting  the  contribution  of  the  balsa  core  to  the  panel  stiffness,  the 
locations  of  the  neutral  axis  (N.A.)  and  the  bending  stiffness  (El)  of  the 
panel  are: 


Y 

NA  EEA 


35.8  mm  , 
(1.41  in.) 


El  -  EEA(y-yNA)2 

•  1.008  x  1010  Nmm2  (3.513  x  106  lb-in.2) 


Mass  per  unic  length  of  the  25.4  mm  (1  in.)  strip  of  panel  between  the 
grounds  must  include  the  mass  of  the  "end  grain"  balsa  pad  and  the  membrane. 
Measured  specific  gravity  of  birch  plywood  is  Y  *  0.72  and  this  value  was 
used  in  place  of  the  values  given  in  Table  D-2  for  birch  and  maple.  For 
a  stainless  steel  membrane  with  a  thickness  of  1.3  ram  (.050  in.),  the  mass 
per  unit  length  is: 


Panel:  m  »  25.4  mm  [(.72  x  10-^  ■^■j)(12  mm  +  3.4  mm)  + 

**  mmJ 

.17  x  10"6  ^-,(155  mm))  -  9.509  x  10~4  ^(1.379  x  10~4lb-sec2/in.2) 
j  mm 

mm 

Pad:  mpD  *  25.4  mm(55  mm) (.17  x  10  ^ 

mm 

=  2,375  x  10~4  &  (3.445  x  10’5  lb-sec2/in.2) 
mm 

Membrane:  m  =  25.4  mm(1.3  mm) (7.84  x  10  ^  ^v>) 
m  J 

A  1^  IQd  n  ty 

=  2.589  x  10  &  (3.755  x  10~5  lb-sec  Via.  ) 
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To  account  for  the  effect  of  submersion  in  LNG  Equations  (V-8)  through 
(V-iO)  are  used  with  a  *  ».  Compute  the  panel  density,  pp,  from  the  mass 
per  unit  length  m^. 
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Density  of  the  LNG 


is 


- — - j  "  *335  x  10’8 

(25.4)(170)mm  ran 


<3.136  x  10 


-7  lb  sec 


Pe  -  .48  x  10'8  (4.49  x  10 

mm 


-7  lb  sec 


in. 

2 


-) 


in. 


Substituting 


K  ■  irh 


'\/~T  +  "T  "  ^0.70  mm)*/0  ■t 
Va^  d  * 


(487.5  mm)' 


-  1.096 

♦■-y/r? 


pl/pp 


.658 


T  -  1.20/. 658  -  1.82  ms 


For  a  tank  which  is  36  m  (118  ft)  long* the  dynamic  load  factor  is  read 
directly  from  Figure  V-2  for  the  fundamental  period  of  1.82  tas. 


DLF  =  1.57 


The  maximum  stress  produced  by  bending  will  occur  in  the  maple  plywood 
and  is 


o 


Me 

El 


where 

c  *  h-Y... 

M  EA 

EM  *  •  T&hnxJ)- 5426  N/-"2  <78-696  ps1) 

M  =  maximum  moment  for  a  clamped  beam 
^2  (b’.Pa.DLF)L2 

"  lF  =  12 
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b*  =  beam  width  =  25.4  mm  (1  in.) 

Solving  for  the  allowable  sloshing  pressure,  P  ,we  have 

3 


p  a  - 

a  b’.DLF.L^h-Yj^Xy 

2 

For  a  permissible  compressive  stress  of  54  N/mm 
allowable  sloshing  pressure  is 


(D-6) 


in  the  maple  plywood, the 


P  *  0.957  N/mra2  (137  psi) 

3 

An  allowable  sloshing  pressure  for  shear  strength  of  the  panel 
will  be  based  upon  a  constant  shear  stress  through  the  core.  The  shear 
reaction  at  each  end  of  the  25.4  mm  (1  in.)  wide  beam  is 


V  -  l/2-(487.5)  (25.4)  (P  ) -1.57 

a 


Using  the  allowable  stress  for  shear  parallel  to  the  grain  (a  conservative 
value),  the  allowable  sloshing  pressure  can  be  found. 


a 


6 


2.07  N/mm2 


v  1/2(487. 5)(25.4)Pa(l. 57) 
A  "  155(25.4) 


P  -  0.84  N/mm2 (122  psi) 

3 


Note  that  the  shear  and  bending  strength  of  the  panel  are  in  good  balance. 

The  numbers  should  be  somewhat  conservative  because  the  shear  strength 
of  the  face  plywood  was  neglected,  the  bending  strength  of  the  balsa  core 
was  neglected  and  the  width  between  the  centerline  of  the  grounds  was 
used  for  the  panel  width. 

D.3  CONCH  Independent  Tank 

The  structural  arrangement  of  the  CONCH  independent  tank  is  described 
in  Section  III .3.  One  of  the  weakest  parts  of  the  tank  is  the  swash 
bulkhead.  In  this  analysis* we  will  predict  the  allowable  sloshing  pressure 
for  a  section  of  the  bulkhead  between  Frames  6  and  7  (refer  to  Figures 
III -23  and  III -24  ). 

A  schematic  of  the  section  between  Frame  6  and  7  is  given  in  Figure 
D-5.  Primary  dimensions  are  in  inches  as  given  in  the  engineering  drawings. 
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FIGURE  D-5 


Stiffened  Panel 


The  fundamental  frequency  of  the  section  of  the  bulkhead  between  the 
frames  and  between  the  brackets  will  be  estitua‘  <u.  using  Equation  vV-3). 
Input  parameters  for  the  equation  are 

h  =  13  mm  (0.50  in.) 

=  71,000  N/mm2  (10.3  x  106  psi) 


E 


p  =  2.65  x  10  u  ^  (2.49  x  10  lb  sec  /in.  ) 
mm 

I  *  61.086  x  106mm4(146.76  in.4) 

S  -  10587  mm2  (16.41  in.2) 


p  ■  733  mm  (28.875  in.) 

L  -  3454  mm  (136  in.) 

o 

Boundary  conditions  for  the  panel  can  be  approximated  by  clamped  boundaries 
at  the  frame  and  simply  supported  along  the  stiffener  to  which  the  brackets 
are  attached.  For  these  boundary  conditions  the  constants  from  Table  V-2 
are 


A'  -  1.5764 
B*  -  4.85 
C'  -  .0833 


Now  substituting  into  Equation  (V-3)  we  find 

X1 


i-"\4 


v!  + 

La2  p  s  V 


T  *  y  *  ms 


The  dynamic  load  factor  for  the  bulkhead  will  be  determined  using  Figure 
C-2(b)  and  the  minimum  nondimen sional  rise  time  T2  [tr  in  Figure  C— 2 ib) ] 
from  Table  V-l.  For  a  nondimensional  time  of  0.0005  and  a  tank  length, 
L,  of  15.608  pi  (51.208  ft)  to  the  swash  bulkhead,  the  rise  time  tr  is 


T^g/L  *  0.0005 


tr  =  T2  =  0.63  ms 


The  ratio  of  the  rise  time  to  the  fundamental  frequency. 


‘r/T  ’  05  ‘  -°76 


is  used  to  enter  Figure  C-2(a).  The  dynamic  load  factor  is 


DLF  ~  2.0 


Strength  of  the  panel  can  be  approximated  by  considering  the  bending  of 
one  stiffener  and  the  associated  plating  between  the  two  frames  for  a 
uniform  pressure  loading  acting  over  the  width  p  (see  Section  V.2.2.2). 
The  loading  per  unit  length,  N,  is 


N  -  p-DLF»P 

a 


where  Pa  is  the  allowable  sloshing  pressure.  Using  the  standard  equation 
for  a  uniformly  loaded  beam  with  clamped  ends  and  taking  the  allowable 
yield  stress,  a  ,  for  5083  aluminum  at  -196°C  as  158.6  N/mm2  (23,000  psi), 
the  allowable  sloshing  pressure  is 

121  °y  ,  2 

P  ■  - "  *034  N/nan  (5  psi) 

pDLF  Ls  C 

where  c,  the  distance  from  the  neutral  axis  to  the  extreme  fiber,  is 
193  mm  (7.585  in.).  This  allowable  sloshing  pressure  is  quite  low  and 
simply  indicates  that  the  tops  of  LNC  tanks  are  designed  for  relatively 
low  static  vapor  pressure  and  small  liquid  loads  and  not  for  dynamic 
sloshing  pressures. 
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APPENDIX  E 


IMCO  TANK  TYPE  DEFINITIONS 


CHAPTER  IV  -  CARGO  CONTAINMENT 


4.1  General 

Administrations  should  take  appropriate  steps  to  ensure  uniformity  in 
the  implementation  and  application  of  the  provisions  of  this  chapter.* 

4.2  Definitions 

In  addition  to  those  in  1.4,  the  following  definitions  applv  through¬ 
out  the  Code. 

4.2.1  Integral  tanks 

(a)  Integral  tanks  form  a  structural  part  of  the  ship’s  hull  and 

are  influenced  in  the  same  manner  and  by  the  same  loads  which 

stress  the  adjacent  hull  structure. 

(b)  The  "design  vapour  pressure"  Pe  as  defined  in  4.2.5  should  not 
normally  exceed  0.25  kp/cm^.  If,  however,  the  hull  scantlings 
are  increased  accordingly,  P0  may  be  increased  to  a  higher 
value  but  less  than  0.7  kp/cm  . 

(c)  Integral  tanks  may  be  used  for  the  products  provided  that  the 

lowest  temperature  in  any  part  of  the  hull  structure  under  no 

circumstances  will  fall  below  -10°C.  A  lower  temperature  may 
be  accepted  by  the  Administration  subject  to  special  consider¬ 
ation. 

4.2.2  Membrane  tanks 

(a)  Membrane  tanks  are  non-self-supporting  tanks  which  consist  of  a 
thin  layer  (membrane)  supported  through  insulation  by  the  adja¬ 
cent  hull  structure.  The  membrane  is  designed  in  such  a  way 
that  thermal  and  other  expansion  or  contraction  is  compensated 
for  without  undue  stressing  of  the  membrane. 

(b)  The  design  vapour  pressure  P0  should  not  normally  exceed  0.25 
kp/cm^.  If,  however,  the  hull  scantlings  are  increased  accord¬ 
ingly,  and  consideration  is  given,  where  appropriate,  to  the 
strength  of  the  supporting  insulation,  PQ  may  be  increased  to 

a  higher  value  but  less  than  0.7  kp/cm^. 

(c)  The  definition  of  membrane  tanks  does  not  exclude  designs  such 
as  those  in  which  non-metallic  membranes  are  used  or  in  which 
membranes  are  included  or  incorporated  in  insulation.  Such 
designs  require,  however,  special  consideration  by  the  Adminis¬ 
tration. 

*  Reference  is  made  to  the  published  Rules  of  members  and  associate  members  of 
the  International  Association  of  Classification  Societies  and  in  particular 
to  IACS  Unified  Requirements  Nos.  92  and  93. 
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4.2.3 


Semi-membrane  tanks 


(a)  Semi-membrane  tanks  are  non-self-supporting  tanks  in  the  loaded 
condition  and  consist  of  a  layer,  parts  of  which  are  supported 
through  insulation  by  the  adjacent  hull  structure,  whereas  the 
rounded  parts  of  this  layer  connecting  the  above-mentioned  sup¬ 
ported  parts  are  designed  also  to  accommodate  the  thermal  and 
other  expansion  or  contraction. 

(b)  The  design  vapour  pressure  P0  should  not  normally  exceed  0.25 
kp/cm2.  If,  however,  the  hull  scantlings  are  increased  accord¬ 
ingly,  and  consideration  is  given,  where  appropriate,  to  the 
strength  of  the  supporting  insulation,  P0  may  be  increased  to 

a  higher  value  but  less  than  0.7  kp/cm2. 


4.2.4 


Independent  tanks 


Independent  tanks  are  self-supporting;  they  do  not  form  part  of  the 
ship's  hull  and  are  not  essential  to  the  hull  strength.  The  three  categories 
of  independent  tanks  are: 


(a)  Independent  tanks  type  A  which  are  designed  primarily  using 
Recognized  Standards  of  classical  ship-structural  analysis 
procedures.  Where  such  tanks  are  primarily  constructed  of 
plane  surfaces  (gravity  tanks),  the  design  vapour  pressure  PD 
should  be  less  than  0.7  kp/cm2. 


(b)  Independent  tanks  type  B  which  are  designed  using  model  tests, 

refined  analytical  tools  and  analysis  methods  to  determine  stress 
levels,  fatigue  life  and  crack  propagation  characteristics. 

Where  such  tanks  are  primarily  constructed  of  plane  surfaces 
(gravity  tanks)  the  design  vapour  pressure  PQ  should  be  less 
than  0.7  kp/cm2. 


(c)  Independent  tanks  type  C  (also  referred  to  as  pressure  vessels) 
are  tanks  meeting  pressure  vessel  criteria  and  having  a  design 
vapour  pressure  not  less  than: 

P0  -  2  +  AC(p)  £ kp/cm2 J 

where 

A  ■  °-0185  (b|) 

with 


Om  =  design  primary  membrane  stress 

AaA  *  allowable  dynamic  membrane  stress  (double  amplitude 
at  probability  level  Q  =  10“®) 

5.5  kp/mm2  for  ferritic/martensitic  steel 

2.5  kp/mm2  for  aluminum  allow  (5083-0) 


c 


»  a  characteristic  tank  dimension  to  be  taken  as  the 
greatest  of  the  following: 

h;  0.75b;  or  0.45& 

with 

h  *  height  of  tank  (dimension  in  ship’s  vertical 
direction)  (m) 

b  *  width  of  tank  (dimension  in  ship's  transverse 
direction)  (m) 

i  m  length  of  tank  (dimension  in  ship's  longitudinal 
direction)  (m) 

p  *  the  relative  density  of  the  cargo  (p  ■  1  for  fresh  water) 
at  the  design  temperature. 

However,  the  Administration  may  allocate  a  tank  complying  with 
the  criterion  of  this  sub-paragraph  to  type  A  or  type  B,  depen¬ 
dent  on  the  configuration  of  the  tank  and  the  arrangement  of 
it 8  supports  and  attachments. 
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